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1.  BACKGROUND 


Bonded  repair  technology  provides  enhanced  structural  efficiency  to  extend  the 
fatigue  life  of  aircraft  structures.  It  provides  industry  with  an  efficient  and  cost 
effective  solution  to  aircraft  structural  defects.  This  has  prompted  many 
organisations  to  use  bonded  repair  technology  as  a  preferred  repair  method  over 
mechanically  fastened  repairs. 

The  USAF  has  utilised  bonded  repair  technology  to  prolong  the  structural  lives  of 
aircraft  such  as  the  F-16  and  the  C130.  But  the  verification  process  of  such  repairs  can 
be  expensive  and  time  consuming,  since  detailed  Finite  Element  Analysis  and  testing 
may  be  required  to  ensure  an  efficient  repair.  The  use  of  validated  analytical 
procedures  for  deriving  repairs  is  therefore  desirable  to  reduce  verification  and 
testing  time.  The  US  Airforce  academy  has  developed  a  software  package  called 
"CalcuRep"  for  this  purpose. 

Robert  Fredell  and  Walter  VanBarneveld  first  developed  CalcuRep  at  Delft 
University  in  the  early  1990's.  Based  on  work  done  by  Robert  Fredell  [1]  at  the  Delft 
University  of  Technology,  CalcuRep  was  designed  to  repair  cracks  specifically  on  an 
Aluminium  aircraft  fuselage.  This  initial  version  of  the  software  was  developed  in 
DOS  and  provided  the  user  with  repair  design  tools  with  a  basic  user  interface.  This 
prompted  the  development  of  CalcuRep  for  Windows  (Beta)  at  tjre  US  Air  Force 
Academy  in  1994,  which  consisted  of  a  user-friendly  menu  system  and  a  graphical 
interface.  CalcuRep  for  Windows  Version  1.0  [2],  released  in  1997,  incorporated  user 
requirements  and  suggestions. 

The  RAAF  Engineering  Standard  C5033  [3,  4,  5]  was  released  in  1994  to  provide 
bonded  repair  design  tools  to  assist  in  the  repair  of  structural  defects  in  aircraft  used 
by  the  RAAF.  Aimed  more  towards  providing  engineers  with  a  set  of  design 
standards  rather  than  a  software  package,  C5033  provides  a  set  of  constraints  to 
ensure  quality  of  repair  design,  specifications  and  application  processes. 

In  early  1999,  the  US  Air  Force  Academy  awarded  a  contract  through  the  DSTO 
Centre  of  Expertise  at  Monash  University  to  verify  and  extend  the  capability  of 
CalcuRep  version  1.0  for  Windows.  The  work  is  noted  as  the  CalcuRep  review  task. 
The  verification  aspect  of  the  task  (know  as  part  1  of  the  work)  was  performed  by 
comparing  the  CalcuRep  approach  with  the  RAAF  Engineering  Standard  C5033 
procedures.  A  significant  example  case  of  a  RAAF  F-111C  wing  repair  was  used  to 
compare  the  CalcuRep  and  C5033  repair  design  techniques.  The  results  are  detailed 
in  this  report. 

Part  2  of  the  task  was  to  extend  the  capability  of  CalcuRep  to  include  unsupported 
single-sided  repairs  and  repairs  to  stiffened  panels.  The  results  for  part  2  are  detailed 
in  Chapter  2. 
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The  work  compiled  in  reference  [8]  establishes  a  new  design  approach  from  the  most 
recent  works  developed  at  the  DSTO  in  the  field  of  bonded  repair  design.  Applying 
these  recent  findings  in  such  a  manner  to  design  repairs  will  prove  to  benefit  both  the 
DSTO  and  the  USAF.  The  content  of  reference  [8]  may  also  be  considered  by  the 
RAAF  in  expanding  the  current  RAAF  Engineering  Standard  C5033  to  deal  with 
single-sided  repairs  to  unsupported  structure.  But  overall,  this  task  proves  to  be 
beneficial  to  both  the  USAF  and  the  DSTO  by  attempting  to  reduce  any  differences 
between  the  US  and  Australian  repair  design  methodologies. 


2.  INTRODUCTION 

Verification  of  procedures  and  mathematical  tools  used  for  bonded  repair  design  is 
essential  in  providing  an  efficient  and  feasible  repair.  The  RAAF  Engineering 
standard  C5033  uses  such  tools  and  was  introduced  to  provide  a  conservative  design 
approach  to  bonded  repairs.  Using  this  standard  in  conjunction  with  a  valid  design 
case  provides  for  a  baseline  repair  design  that  can  be  used  to  validate  the  CalcuRep 
software  package.  The  particular  case  study  used  for  this  purpose  is  the  1994  F-111C 
Lower  Wing  Skin  Boron  Repair. 

Part  1  of  the  CalcuRep  review  consists  of  two  steps  as  follows: 

a  A  detailed  review  of  the  CalcuRep  approach  to  repair  design  and' analysis,  and 

b  A  case  study  comparison  using  CalcuRep  to  analyse  a  bonded  composite  repair  to 
the  lower  wing  skin  of  a  RAAF  F-111C  aircraft,  originally  designed  according  to  the 
RAAF  Engineering  Standard  C5033  [3]. 

Several  errors  and  limitations  were  identified  and  corrected.  A  revised  version  of  the 
CalcuRep  code  was  then  implemented  in  a  Mathcad  module  [9]  and  used  to  analyse 
the  F-lll  repair  case.  The  results  are  detailed  in  this  report. 


3.  EXISTING  TOOLS  AND  TECHNIQUES 

3.1  RAAF  Engineering  Standard  C5033 

3.1.1  Overview  of  the  C5033  Approach 

Designing  an  effective  repair  is  essential  in  extending  aircraft  structural  life.  Assuring 
the  quality  of  repair  design,  specifications  and  processes  is  fundamental  in  achieving 
this  goal.  The  RAAF  engineering  standard  C5033  [3]  in  conjunction  with  the 
following  two  documents,  AAP  7021.016-1  [4]  and  AAP  7021.016-2  [5]  has  been 
developed  to  manage  these  issues.  Mr  Max  Davis,  a  Senior  Research  Scientist  for  the 
RAAF,  conducts  two  tuition  courses  on  the  Standard.  The  first  is  a  familiarisation 
course  centred  towards  providing  the  basic  concepts  behind  repairs  and  processes. 
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while  the  second  details  repair  design,  specification  and  process  techniques.  The 
courses  are  held  twice  annually  at  the  RAAF  Base  Amberley,  Queensland. 

C5033  provides  design  methods  for  the  following  repair  cases,  assuming  the  patch 
and  parent  structure  to  be  both  isotropic: 

•  Double  Sided  Bonded  Repairs 

•  Single  Sided  Supported  Bonded  Repairs 

•  Single  Sided  Supported  Bonded  Repairs  to  a  Part  Through  Thickness  Defect 

•  Single  Sided  Supported  Bonded  Doubler 

•  Double  Sided  Bonded  Doubler 

•  Single  Sided  Supported  Bonded  Repairs  for  Stress  Corrosion  Cracking 

•  Double  Sided  Bonded  Repairs  for  Stress  Corrosion  Cracking 

•  Scarf  Repairs  to  Composite  Structures 

C5033  does  not  equip  the  user  with  a  software  package.  Instead  it  provides  strict 
guidelines  on  repair  design,  specification  and  processes.  Therefore  the  user  must 
perform  repair  design  calculations  either  manually,  or  by  programming  them  into  an 
appropriate  software  package.  MathCad  version  8.0  [9]  was  used  in  this  report  for 
this  purpose,  since  it  provides  visual  inspection  of  the  equations  used  and  quick 
accessibility  of  user  input  data.  ' 

An  overview  of  C5033's  design  methodology  for  single-sided  supported  repairs  to  a 
cracked  structure  is  presented  below: 

3.1.1.1  Determination  of  Patch  Thickness 

Restoration  of  stability  and  compatibility  characteristics  of  a  cracked  or  cut-out 
structure  is  a  result  of  restoring  the  original  stiffness  of  the  structure.  C5033  restores 
the  original  stiffness  of  the  structure  by  providing  the  following  limitation  on  patch 
design  thickness. 


Eoto 


Patch 


Eih 


Figure  1  Parent  Material  or  Structure 


Eit;  ^  1.2 Eiti 

<tn< - 

Eo  E0 


(1) 


Flere  Ei  Elastic  Modulus  of  the  Parent  Material 

E0  Elastic  Modulus  of  the  Patch 
ti  Thickness  of  the  Parent  Material 
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to 


Thickness  of  the  Patch 


3.1.1.2  Evaluate  Rapid  Repairability  Criterion  (RRC) 

The  rapid  repairability  criterion  is  a  management  tool  to  ensure  the  adhesive  is  never 
the  critical  element  in  a  repaired  joint. 


Strength 


Figure  2 


Provided  the  adhesive  operates  within  region  'A'  shown  above,  the  parent  material 
will  fail  before  the  adhesive  does.  On  the  other  hand,  adhesives  operating  in  section 
B  above  will  fail  before  the  parent  material  does.  An  adhesive  will  never  be  the 
critical  element  (therefore  operates  within  region  A  above)  in  a  repaired  joint 
provided  its  load  capacity  can  meet  the  following  criteria. 

P>1.2xauuriXti  (2) 


where  'F  is  the  lesser  of 


P  =  2Ef)[a0  - aieff\RT  - Tcure)+ {a0  -  a^Toper  - J?r)]+ 


~re+rP 


1  + 

V  E0t0  J 


or 


P  =  2 E0t0\aieff  -  a0\RT  -  Tcu n)+(a,  -  -  i?r)]+  ^2rjrf 

_afi  +  vI) 


f\ 

v2 


-Ye+Yp 


>  (  E  t  ^ 

Et  1  +  ^ 

l  Efij 


aieff  ~ 


(3) 

(4) 

(5) 


Here  P 

do 

CXi 

Otieff 

RT 

Tcure 


Load  capacity  of  the  adhesive 

Coefficient  of  Thermal  Expansion  of  the  Patch 

Coefficient  of  Thermal  Expansion  of  the  Parent  material 

Effective  Coefficient  of  Thermal  Expansion  of  the  Parent  material 

Room  Temperature 

Cure  Temperature 


\ 


5 


Toper  Operating  Temperature 
r|  Adhesive  Thickness 

xp  Adhesive  Plastic  Shear  Stress  Limit 

7e  Adhesive  Elastic  Shear  Strain  Limit 

yp  Adhesive  Plastic  Shear  Strain  Limit 

cruLTi  Ultimate  Strength  of  the  Parent  Material 

Equations  (3)  and  (4)  where  developed  by  Hart-Smith  [10]  and  include  a  bi-linear 
thermal  residual  stress  modification  developed  by  Mr  Max  Davis  to  allow  for  both 
the  cure  and  aircraft  operating  cycles,  and  will  be  discussed  in  more  detail  in  Section 
3.1.3. 


3.1.1.3  Assess  Structural  Integrity  of  the  Repaired  Structure 
Step  1  Calculate  the  Applied  Stress  (a*) 


Figure  3 

The  applied  stress  is  determined  as  follows: 

•  Where  the  design  ultimate  load  is  known,  it  must  be  used. 

•  Where  the  design  limit  load  is  known,  the  applied  stress  will  be  1.5  times  the 
design  limit  load. 

•  Where  the  design  ultimate  and/or  the  design  limit  load  are  not  known,  the 
applied  stress  is  taken  as  the  lower  value  of  material  ultimate  strength  or  1.5 
times  the  yield  strength.  Load  path,  stress  concentrations,  net  section  stresses  and 
buckling  limitations  must  be  taken  into  account. 

Step  2  Assess  Structural  Integrity  at  the  end  of  the  Patch 


Figure  4 

The  stress  in  the  structure  at  the  end  of  the  patch  can  be  expressed  by  the  following: 
a"  =  £2,0-'  +  £,((a„ -a,€\RT-Tm)+(a0 -Rt))  (6) 


where,  Ql  is  the  load  attraction  factor 

Equation  6  consists  of  two  components,  an  applied  stress  component  as  per  [11],  and 
a  thermal  residual  stress  component  as  per  [12]. 
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The  Load  Attraction  Factor 


Load  is  attracted  into  the  repaired  region  due  to  the  inclusion  of  extra  stiffness 
provided  by  the  patch  as  shown  below: 


Figure  5 


In  a  uniaxial  load  case,  the  following  stresses  are  attracted  into  the  repaired  region, 
which  is  denoted  as  section  A- A. 

1st  Principle  Stress  Direction 


Figure  6 

Section  A-A  is  shown  below 


1st  Principal  Stress  Direction 
a*  +  p 


■  2nd  Principal  Stress  Direction 

q 


Figure  8 


Here  'p'  and  'q'  are  the  resultant  stress  components  attracted  into  the  repaired 
region.  Therefore  the  load  attraction  factor  becomes  [11]: 
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(7) 


* 
a 

C5033  does  not  calculate  this  value,  but  instead  uses  a  typically  conservative  value 
for  the  load  attraction  as  follows: 


Ql=1.2 

The  required  margin  of  safety  for  a**  is  as  follows: 

CfULTi/  o”  -1>0.2 

Step  3  Determine  the  Stress  Under  the  Repair 

The  stress  in  the  structure  at  the  edge  of  the  patch  is  distributed  between  the  parent 
material  and  the  patch  inside  the  repaired  region.  The  corresponding  stress 
distributed  into  the  parent  structure  [11]  is  as  follows: 


cr  = 


E,t 


E,t  +  E,t; 


(8) 


B  A 


B  A 
Figure  9 


Here  section  A-B  is  as  shown  below 


Figure  10 


where  o0  is  the  stress  under  the  repair  in  the  absence  of  a  flaw 

or  is  the  stress  in  the  repair  in  the  absence  of  a  flaw 


Step  4  Evaluate  the  Adhesive  Shear  Strain 
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The  adhesive  shear  strain  calculation  relies  on  an  elastic  perfectly  plastic  shear 
stress/ strain  curve  assumption.  This  assumption  is  reasonable  provided 
conservation  of  strain  energy  is  achieved  between  the  true  shear  stress/  strain  curve 
and  the  equivalent  elastic  perfectly  plastic  shear  stress/  strain  curve,  see  Figure  11. 


Equal  Areas 


Figure  11 

where  yma*  =  ye  +  yp 

The  maximum  elastic  shear  strain  [13]  is  given  as  follows: 


Jo  = 


(9) 


(10) 


where,  X  Elastic  Exponent 

G  Adhesive  Shear  Modulus 


The  elastic  Exponent  is  given  by  [13]: 


r| 


1  1 

+  ■ 


E0t0  j 


(11) 


If  yo^YE,  the  adhesive  is  plastic  and  equation  10  is  replaced  by  the  following: 


Yo  = 


2  G 


f 

(  \ 

2^ 

1  + 

g  oti^ 

V  XP  J 

V 

J 

(12) 


It  is  required  that  —  +  —  - 1  >  0.2 

Vo 


Step  5  Determine  the  Reduced  Stress  Intensity  Factor  in  the  Structure  Post  Repair. 
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Reference  [11]  demonstrates  that  the  stress  intensity  factor  for  an  adhesive  bonded 
repair  never  exceeds  a  limiting  value. 


where  Kr  Stress  Intensity  Factor  (SIF)  after  repair 
Ko  Stress  Intensity  Factor  before  repair 
Kco  Upper  Bound  or  Limiting  Stress  Intensity  Factor  . 

•  Characteristic  Crack  Length 

The  characteristic  crack  length  is  a  physical  parameter  independent  of  actual  crack 
length  and  applied  stress  and  provides  the  upper  bound  stress  intensity  factor 
through  the  following  relationship. 

Koo=a0(7tA)1/2  (13) 

Depending  on  whether  the  adhesive  is  elastic  or  plastic,  the  following  two  equations 
are  commonly  used  to  predict  this  upper  bound  stress  intensity  factor. 

Elastic  Behaviour  in  the  adhesive  [11] 


Plastic  Behaviour  in  the  adhesive  [13]  &  [14] 


K  = 


<7Jn 

f 

1  + 

"cr.V 

< 

( 

1  +  2 

'  °oXti 

3  ^ 

G 

o  p 

V 

\  XP  J 

\ 

) 

3  M, 

l  TP  J 

) 

(15) 
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It  is  required  that  —  - 1  >  0.2 . 

Where  Kc  Parent  Material  Fracture  Toughness. 

Step  6  Calculate  Peel  Stress 

Peel  stresses  occur  due  to  load  path  eccentricity  inherent  in  single-sided  repairs. 
The  maximum  peel  stress  likely  to  occur  in  a  single-sided  repair  is  [10]: 


where  oc  max  Maximum  Peel  Stress 

Ec  Effective  Transverse  Tensile  Modulus  for  the  Adhesive 

The  effective  transverse  tensile  modulus  of  an  adhesive  for  a  joint  comprising  of 
dissimilar  materials  can  be  approximated  by 

1_  A_ 

'  E'c  Ec  +  E,  +  E0 


(17) 


where  Ec  Elastic  Modulus  of  the  Adhesive  System. 

Reliable  peel  stress  data  is  not  readily  available.  There  is  also  no  current  test 
procedure  to  determine  a  valid  peel  stress  allowable,  therefore  Hart-Smith  [15] 
devised  the  following  general  guild-line  values: 

crc  max  ^  6000  psi  (Brittle  Adhesive  Systems) 
ac  max  <  10000  psi  (Ductile  Adhesive  Systems) 

Step  7  Assess  the  Structural  Integrity  of  the  Patch 

The  structural  integrity  of  a  patch  is  found  assuming  the  entire  load  is  transferred 
into  the  patch.  This  assumption  is  valid  provided  the  actual  crack  length  is 
sufficiently  larger  than  the  characteristic  crack  length  [11]: 
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3.1.2  Evaluate  Repair  and  Structure  for  Fatigue  Susceptibility 


This  section  is  used  to  eliminate  fatigue  critical  repairs.  Note,  it  does  not  provide  a 
crack  growth  calculation.  To  deem  a  repair  fatigue  uncritical,  the  repair  must 
demonstrated  fatigue  resistance  when  subject  to  a  constant  amplitude  fatigue  loading 
with  a  maximum  cyclic  load  (of)  equal  to  the  lower  of: 

•  60%  of  DLL  (where  known)  or 

•  40%  of  cr* 

at  a  stress  ratio  of  R=0. 


Fatigue  susceptibility  is  determined  as  follows. 

•  Step  1  Evaluate  the  Fatigue  Susceptibility  of  the  Structure  at  the  end  of  the 
Patch. 

where  a##  is  the  stress  at  the  edge  of  the  repair  due  to  the  fatigue  loading. 

Check  if  reference  data  contains  anticipated  damage  at  a  stress  level  of  a##.  If  so, 
check  whether  this  level  of  damage  is  acceptable  for  the  corresponding  repair  case. 

•  Step  2  Stress  Under  the  Patch 


a  = 


El, 


.## 


EJ.  +  EJ, 


Where,  a#  Stress  under  the  repair  due  to  the  fatigue  loading. 


(20) 


•  Step  3  Calculating  the  Adhesive  Shear  Strain  at  the  Defect 


(21) 


If  y0f>YE,  the  adhesive  has  become  plastic  and  the  adhesive  shear  strain  is  given  by: 


f 

cr#i/ 0 

2\ 

1  + 

l 

T 

V 

\  P  ) 

J 

(22) 


Where,  yt  Maximum  Shear  Strain  in  the  adhesive  due  to  the  fatigue  loading. 
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The  adhesive  is  said  to  be  fatigue  resistant  if  the  following  condition  is  met. 

Yf  <  2yE 


•  Step  4  Calculate  the  Stress  Intensity  Factor  in  the  Repair  Parent  Material 

Elastic  Behaviour  in  the  adhesive 


Plastic  Behaviour  in  the  adhesive 


(23) 


where  Koof,  is  the  upper  bound  SIF  due  to  the  fatigue  loading 

\ 

Using  crack  growth  data  for  the  parent  structure,  compare  the  reduction  in  the  crack 
growth  rates  between  the  SIF  produced  above  and  the  unrepaired  SIF.  An  extension 
to  allow  for  visplastic  adhesive  behaviour  is  given  in  Chapter  3. 

•  Step  5  Check  the  Fatigue  Resistance  of  the  Patch 


(25) 


where  crpf  Maximum  stress  in  the  patch  due  to  the  fatigue  loading. 

Check  if  reference  data  contains  anticipated  damage  at  a  stress  level  of  apf.  Determine 
whether  this  level  of  damage  is  acceptable  for  the  repair  case. 

3.1.1.4  Establish  Patch  Dimensions 

Patch  Length 

The  elastic  load  transfer  length  is  the  inverse  of  the  elastic  exponent  (Equation  11).  To 
provide  creep  resistance,  a  factor  of  3  has  been  introduced  as  shown  in  equation  27. 
The  plastic  load  transfer  length  [18]  has  been  derived  to  be  sufficiently  large  enough 
to  carry  all  of  the  unnotched  parent  material  ultimate  stress.  The  total  length  of  the 
patch  must  include  any  pre-existing  bond-line  defects  such  as  fastener  heads. 
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LT  =le  +  Ip  (26) 

‘e  =  l  (27) 

,  .  WJL  (28) 


where  Lt  Total  transfer  length 

L  Overlap  Length 

le  Elastic  Load  Transfer  Length 

lp  Plastic  Load  Transfer  Length 

8  Defect  Size 


•  Patch  Width 

h-  Width  "H 


Figure  14 
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3.1.1.5  Determine  Inspection  requirements 

In  order  to  ensure  the  adhesive  is  never  the  critical  element,  the  total  length  of  defects 
in  a  given  line  must  be  less  than  or  equal  to  (j>max: 

*™=A+0-u  <29> 

where  lmm  =  l-5x  21  p  (30) 

If  the  defect  size  is  greater  than  this  value,  the  joint  will  not  be  able  to  transfer 
material  ultimate  stress  through  the  plastic  zone  at  all  times.  A  margin  of  safety  of 
50%  is  recommended. 

3.1.2  Limitations  of  C5033 

Two  major  limitations  of  the  RAAF  Engineering  Standard  C5033  are  the  inability  to 
deal  with  the  following  scenarios: 

•  Strong  Biaxial  Loading  Applications 

•  Unsupported  single  sided  repairs 

•  Limited  treatment  of  interlaminar  failure  mechanisms  in  the  patch 

•  Limited  treatment  of  visco-plastic  behaviour  in  the  adhesive.  * 

3.1.2.1  Strong  Biaxial  Loading  Applications 

If  the  parent  structure  and  repair  material  have  identical  Poisson's  ratios  and  are 
assumed  to  be  isotropic,  2nd  direction  principal  stresses,  i.e.  ay+q  (see  Figure  15),  do 
not  influence  the  stress  in  the  major  direction  under  the  patch.  However,  if  the 
Poison's  ratios  differ  in  a  significant  biaxial  stress  field,  the  stress  under  the  repair 
will  be  effected  significantly  by  both  principal  stresses. 


Cy+q  Major  Direction 

t  t  t  - ’ 

◄ 

G*+p  * 

i  i  1 


a*+p 


cjy+q 

Figure  15 

where,  ay  is  the  applied  stress  in  the  2nd  principal  stress  direction 
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Equation  (8)  produces  a  good  approximate  solution  for  uniaxial  load  cases  when  the 
parent  structure  and  repair  consist  of  different  Poisson's  ratios.  For  similar  Poisson's 
ratios,  equation  8  produces  the  exact  solution. 

3.1.2.2  Unsupported  Single-Sided  Repairs 

When  a  single  sided  repaired  structure  is  unsupported,  out  of  plane  bending  effects 
are  present.  These  effects  cause  the  stress  inside  the  repaired  region  to  become  non- 
uniform,  see  Figure  16. 


Resulting  Stress 


Figure  16 


Applied  Stress 


The  C5033  approach  does  not  take  these  effects  into  account  therefore  is  limited  to 
only  supported  repair  cases. 

3.1.3  Improvements  made  Since  Release  of  the  Original  Document 

t 

The  RAAF  standard  C5033  is  a  living  document  and  changes  when  new 
developments  are  made  which  can  benefit  the  repair  process.  Improvements  made  to 
previous  editions  of  C5033  are  as  follows: 


•  Load  Capacity  of  the  Adhesive 


1994  Draft  C5033  Document 


P  =  2 (a0  -  a^ATE^  +  jlt]T , 

and 


1 

tYe  +  Yp 
V2  j 


'  f  E,t;  ^ 


Eft 


1  4.  _i_L 
V  EOtO  ) 


(31) 


1  ^ 

f1  Et) 

tYe  +Yp 

Et 

O  0 

1+  00 

U  ) 

l  Ef,  J 

(32) 


Differences  between  equations  3,  4,  31  and  32  occur  in  the  thermal  residual  stress 
component.  Equations  31  and  32  contain  a  thermal  residual  stress  for  the  cure  and 
operational  cycles  as  one  cycle,  ie  2(a0-ai)ATEiti  from  (31),  where  AT=TCUre-RT.  While 
equations  3  and  4  deal  with  both  the  cure  and  operational  thermal  cycles 
independently,  i.e.  2Eiti((a0-aieff)(RT-Tcure)+(ao-ai)(Toper-RT)),  developed  by  Mr 
Max  Davis  [3].  ^  j 


Cure  Cycle  Operational  Cycle 
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During  the  cure  cycle,  the  repaired  region  is  heated.  Constraint  is  provided  from 
cooler  material  surrounding  the  heated  region,  which  effectively  reduces  the 
coefficient  of  thermal  expansion  of  the  parent  material.  Therefore  the  aieff  term  is 
included  in  the  cure  cycle  term  above.  On  the  other  hand,  when  the  aircraft  is  in 
operation,  the  repaired  structure  is  cooled  uniformly  at  altitude,  therefore  no 
constraint  against  expansion  is  provided.  Therefore  an  effective  coefficient  of  thermal 
expansion  for  the  parent  structure  is  not  necessary  as  depicted  in  the  operational 
component  of  equations  3  and  4. 

•  Stress  at  the  Edge  of  the  Patch 

1994  Draft  C5033  Document 

<t"  =  QLcr*  +  QrsEt  (at  -  a0  )A T  (33) 

where  Qt  is  the  thermal  inclusion  factor 

Similarly,  differences  occur  between  equation  6  and  33  due  to  the  inclusion  of  an 
independent  operational  thermal  stress  cycle.  Including  this  cycle  eliminates  the  need 
for  the  thermal  inclusion  factor  in  the  current  equation. 

3.2  CalcuRep  , 

3.2.1  Overview  of  the  CalcuRep  Approach 

CalcuRep  was  designed  to  provide  users  with  analysis  tools  and  a  user-friendly 
graphical  interface  to  assist  in  the  design  of  adhesive  bonded  repairs.  It  has  been 
designed  specifically  for  single-sided  supported  repairs  on  aircraft  fuselage,  although 
it  can  be  adapted  to  design  repairs  to  most  locations  on  aircraft.  CalcuRep  considers 
post  repair  stresses  in  a  more  detailed  approach  to  that  of  C5033  for  single-sided 
supported  repairs.  It  attempts  to  expand  the  repair  design  methodology  to  biaxial 
loading  systems  to  model  loads  present  in  aircraft  fuselage,  assuming  the  patch  to  be 
orthotropic  and  the  parent  structure  to  be  isotropic. 

When  the  software  is  run,  the  user  is  prompted  graphically  with  a  choice  of  four 
input  buttons.  The  first  three  buttons  allow  the  user  to  input  skin,  patch  and  adhesive 
data,  manually  or  by  the  use  of  CalcuRep' s  material  data  library.  The  fourth  button 
allows  the  user  to  input  repair  geometry  and  loading  parameters.  Once  the  data 
entry  is  complete,  the  user  can  then  choose  the  'Calculate'  button.  CalcuRep  then 
provides  a  brief  summary  of  total  stresses,  strains  and  stress  intensity  factors  within 
the  repaired  region  along  with  their  respective  margins  of  safety.  The  user  also  has 
the  option  of  viewing  a  detailed  four  page  report  on  the  calculated  applied,  thermal 
and  bending  stress  contributions  of  the  parameters  previously  discussed  and 
®^§&®®tmns  on  what  to  change  to  reducing  stresses  after  repair. 
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The  following  sections  provide  an  overview  of  CalcuRep's  analytical  methodology. 

3.2.1.1  Determining  the  Stress  at  the  Edge  and  Under  the  Repair 

CalcuRep  uses  results  derived  from  [11]  to  determine  the  stress  beneath  and  at  the 
edge  of  the  repair.  Reference  [11]  considers  a  single-sided  supported  repair 
containing  an  infinitely  large  orthotropic  plate  with  an  attached  elliptical  orthotropic 
repair,  see  Figure  17. 


Figure  17 


where,  2a  Patch  length 
2b  Patch  width 

m  =  patch  aspect  ratio  =  a/b 
P  applied  stress  in  the  1st  principle  direction 
Q  applied  stress  in  the  2nd  principal  direction 


The  first  step  is  to  consider  the  repaired  region  as  an  inclusion,  ie;  as  one  separate 
entity  with  the  following  load  applied. 

1st  Principal  direction  2nd  Principal  direction 


P+p 


Q+q 


Q+q 


Figure  18 
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where  'p'  and  'q'  are  stresses  attracted  into  the  repaired  region  (V  &  'y'  directions 
respectively). 

Determining  the  displacement  functions  for  an  infinite  plate  loaded  internally  by 
tractions,  tractions  being  functions  of  'p'  and  'q'  [19],  and  the  inclusion  loaded  as 
shown  above,  see  Figure  18,  matching  the  displacement  in  the  V  and  'y'  directions  at 
the  hole  boundary  produces  the  following; 

(For  an  isotropic  plate) 


<X  +  P+  1 
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1 
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where 


a  imaginary  part  of  the  first  root  of  X 

P  imaginary  part  of  the  second  root  of  X 


(E, 

u 


■2v, 
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+  1  =  0 


X  characteristic  equation  parameter 

Ex  Plate  Modulus  in  x-direction  (isotropic  plate) 

Exi  Inclusion  Modulus  in  x-direction 

Eyi  Inclusion  Modulus  in  y-direction 

Vxi  Inclusion  Poisson's  Ratio  in  xy-direction 

V12  Poisson's  Ratio  of  the  parent  material 


(34) 


(35) 


•  Inclusion  Material  Constants 


Using  the  elastic  properties  below  (Equation  36),  and  the  following  relations  within 
the  repaired  region,  (P+p)t=a0t+aRtR  and  ei=e0=SR,  holding  for  arbitrary  strains;  the 
inclusion  properties  are  derived  as  follows: 


A,  = 


E, 


&  Axr  = 


E. 


xR 


11  XiX  i 

-v12  1 


V\2RV2\R 


A  =  x 
Ay  ■,  2 

1- V, 


&  AyR  = 


JyR 
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1  V\2RV2\R 


where  Ax,  Ay,  Axr,  AyR  Elastic  properties  for  the  plate  and  repair 
Exr  Elastic  modulus  of  repair  in  x-direction 

-  EyR  Elastic  modulus  of  repair  in  y-direction 


(36) 
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V12R 

V21R 

G  &  Gr 


Repair  Poisson  ratios  in  xy-direction 
Repair  Poisson  ratios  in  yx-direction 
Shear  Modulus  of  plate  and  repair,  respectively 


The  inclusion  has  the  following  elastic  properties 


+  AjR  &  A  1  _  Ayt  +  AyRtR 
v  _  V12 Ayt  +  V\U(AyRtR 


'xl 


y  _  V12^/  +  V21  rAxr[r  ^29) 


V 


+  AxRtR 


(37) 


=  4/0  -  lw)&£,/  =  4/0  -  >w)  (40> 


where,  Axi  &  Ayi  Elastic  properties  of  the  inclusion  in  the  'x'  &  'y' 

directions 

Vyi  inclusion  Poisson's  ratio  in  the  yx-direction 

Completing  the  analysis,  see  [11]  for  details,  the  stress  at  the  edge  of  the  repair  and 
under  the  repair  are  as  follows: 

\ 

•  Stress  at  the  Edge  of  the  Repair 

&IXX  =  P  +  P  (41) 

°yy=Q  +  <l  (42) 

stress  in  the  parent  structure  at  the  edge  of  the  repair  in  the 
x-  direction. 

stress  in  the  parent  structure  at  the  edge  of  the  repair  in  the 
y-direction. 

•  Stress  Under  the  Repair 


where  CTxx1 

Gyy1 
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(43) 


where  a0  is  the  stress  under  the  repair  in  the  x-direction 

3.2.1.2  Thermal  Residual  Stresses 

CalcuRep  uses  an  approach  developed  by  Robert  S.  Fredell  in  reference  [1],  which 
centres  on  satisfying  various  equilibrium  and  compatibility  requirements  to 
determine  resulting  stresses  in  the  repaired  and  surrounding  regions  of  the  structure, 
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post  cure.  The  following  entails  the  process  used  by  CalcuRep.  (Note .a  separate  set  of 
notation  will  be  used  for  this  section  to  keep  it  consistent  with  the  CalcuRep  code 
and  is  summarised  in  the  appendix.) 

•  Geometry  as  Used  in  CalcuRep 


CalcuRep  manipulates  the  above  parameters  prior  to  equation  substitution  as 
follows: 

dx=dx/ 2,  dy=dy/2, ,  Rx=1.2Rxand  Ry=Ry/2 


where 


dx  Heated  region  length  in  the  x-direction 

dy  Heated  region  length  in  the  y-direction 

Rx  Half  frame  spacing 

Ry  Stringer  spacing 


•  Effective  Coefficient  of  Thermal  Expansion  of  the  Parent  Structure 

As  stated  previously,  during  the  cure  process,  cooler  material  surrounding  the 
heated  area  provides  constraint  against  expansion,  thus  effectively  reducing  the 
coefficient  of  thermal  expansion  for  the  parent  structure.  CalcuRep  uses  the 
following  equation  specifically  for  aluminium  structure  as  per  [1]. 
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where  ap(eff)x  &  ap(eff)y  Effective  CTE  of  the  parent  structure  in  the  V  and  'y' 

directions 

a  Parent  structure  coefficient  of  thermal  expansion  (CTE) 

Vxy  Poisson's  ratio  of  the  parent  structure 


(46) 

(47) 


Tx  &  Ty  Temperature  at  the  edge  of  the  parent  structure  at  Rx  &  Ry 

Tc  Cure  temperature 

Tamb  Ambient  temperature 

t  Thickness  of  parent  structure 


•  Loads  in  the  Repaired  Structure  after  Cure. 

As  the  repaired  region  is  heated,  the  parent  structure  experiences  initial  thermal 
stresses,  while  the  patch  is  free  to  expand.  Once  the  adhesive  is  cured,  a  rigid  bond  is 
assumed  to  exist  between  the  patch  and  plate.  The  patch  and  structure  is  then  cooled 
to  room  temperature  leaving  thermal  residual  stresses  in  the  patch  and  parent 
structure.  Fredell  [1]  solves  for  unknown  loads  per  unit  width  acting  on  the 
following  regions  within  the  repaired  structure  on  cooling  after  cure,  by  considering 
equilibrium  and  compatibility  conditions: 

Enforcing  the  Rigid  Bond  Assumption  Yields  (noted  as  Pi  loads) 

After  cure,  the  patch  is  assumed  to  be  rigidly  bonded  to  the  structure.  The  following 
loads  ensure  this  assumption  is  valid. 
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where 


PlxP=-PlxR-Plx 

PlyP=-PlyR=Ply 


_  V 


*+*t. 


4 

B3 — 2. 


and  P[x  =  P[  y  —  +  — 


(48) 


Pixp  &  Piyp  Load  applied  to  plate  to  satisfy  the  rigid  bond 

assumption  in  the  'x'  &  'y'  directions  respectively 
PixR  &  PiyR  Load  applied  to  repair  to  satisfy  the  rigid  bond 

assumption  in  the  V  &  'y'  directions  respectively 
Pix  &  Piy  Load  tractions  equal  in  magnitude  to  the  above 
loads  (V  &  'y'  directions  respectively) 

The  remaining  constants,  numbered  as  per  [1],  are: 


B-Xt  ExRtR 

Ba  =  (<xP(eff)x 


B2=  —  +  -^~  53=  — +  — — 
Bxt  BxRtR  Ext  EyRtR  . 

~  axR  B5  =  {aP(eg-)y  -  ayR  )a7 


(49) 


where 


Ex  Elastic  Modulus  of  plate  in  x-direction  (same  as  ‘y’  since 
Isotropic) 

Exr  Elastic  Modulus  of  repair  in  the  x-direction 

EyR  Elastic  Modulus  of  repair  in  the  y-direction 

VxyR  Poisson's  ratio  of  repair  in  the  xy-direction 

VyxR  Poisson's  ratio  of  repair  in  the  y x-direction 

vxy  Poisson's  ratio  of  parent  structure 

aXR  Coefficient  of  thermal  expansion  of  repair  in  the  x-direction 
ayR  Coefficient  of  thermal  expansion  of  repair  in  the  y-direction 


External  Loads  on  the  Inclusion  and  Surrounding  Matrix  (noted  as  P2  loads) 


Figure  21 
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Loads  in  the  Repaired  Structure  During  Operational  Altitude. 

As  discussed  previously,  at  operational  altitude  the  aircraft  structure  is  cooled 
uniformly,  therefore  no  constraint  is  provided  against  contraction  of  the  parent 
structure.  This  results  in  an  additional  thermal  cycle  which  effects  the  resulting 
residual  stresses  in  the  patch  and  the  parent  structure.  The  operational  thermal  stress 
component  is  found  by  using  the  same  procedure  listed  in  the  previous  section  with 
the  following  modifications: 


AT-Toper-Tamb  (61) 

ap(eff)x=a  (62) 


where,  Tambis  the  temperature  at  sea  level  during  operation. 

The  total  loads  due  to  both  cure  and  operational  thermal  cycles  are  determined  by 
simply  adding  the  respective  loads  found  for  both  cycles. 

Converting  Loads  into  Stresses  . 

Once  the  loads  are  determined  within  various  locations  of  the  repaired  structure,  see 
Figures  20  to  22  ,  CalcuRep  calculates  the  resulting  stresses  in  the  parent  structure 
under  the  repair,  at  the  edge  of  the  plate  and  in  the  repair  as  follows:' 


Stress  under  the  repair  in  the  parent  structure 

(Plx+P2xFl)/ 1 

(63) 

Stress  at  the  edge  of  the  repair 

P2x/t 

(64) 

Stress  in  the  repair 

(P  2xRI"P lx)/  tR 

(65) 

3.2.1.3  Out  of  Plane  Bending  Stresses 

Out  of  plane  bending  stresses  exist  due  to  load  path  eccentricity  predominantly 
found  in  single  sided  repairs.  As  a  result  of  this  phenomenon,  the  applied  load  on  the 
structure  tends  to  shift  towards  the  neutral  axis  as  shown  below: 
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Here  P  is  the  applied  stress  and  zx  defines  the  location  of  the  neutral  axis. 

Fredell  [1]  used  a  similar  approach  to  L.R.F  Rose  [13],  which  considers  the  patch  and 
structure  as  a  one-dimensional  infinitely  wide  strip  with  an  applied  load  as  shown 
above.  Using  a  simple  relationship  between  bending  moment  and  deflection  [1],  the 
maximum  stress  at  the  edge  of  the  patch  is  as  follows: 


where 


° Bend  max 


A  3>/Q^-  tanh(Q/g) 

2  <t>^  +  <E>/  tanh(0/a) 


where  OBendmax  Bending  stress  at  the  edge  of  the  repair 

Eli  Bending  stiffness  per  width  of  the  inclusion 

EIm  Bending  stiffness  per  width  of  the  plate  outside  the 

repaired  region 

®m  Root  of  the  characteristic  equation  for  the  plate 

®!  Root  of  the  characteristic  equation  for  the  inclusion 

CalcuRep  uses  the  following  equation  to  estimate  the  reduction  in  bending  stress  at 
■  the  edge  of  the  patch  due  to  patch  taper  ratio: 


Figure  24 
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where 


tRmin 

CTBendmin 


Minimum  patch  thickness 

Reduced  bending  stress  due  to  patch  taper 


3.2.1.4  Remaining  Parameters 

Once  the  stresses  have  been  determined  CalcuRep  then  calculates  the  resulting 
reduced  stress  intensity  factor  after  repair,  the  maximum  stress  in  the  patch  and  the 
maximum  adhesive  shear  strain. 

•  Reduced  Stress  Intensity  Factor  after  Repair 

Firstly,  CalcuRep  determines  the  total  stress  in  the  parent  structure  under  the  repair 
by  a  summation  of  the  following; 

Co  New=  (o0+  (Plx+P2xPl)  / 1)  (73) 

where  o0  New  is  the  total  Stress  under  the  repair  as  used  in  CalcuRep.  The  term  ct0 
New  replaces  a0  in  either  equation  14  or  15  depending  on  whether  the  adhesive  is 
elastic  or  plastic,  thus  producing  the  reduced  stress  intensity  factor. 

•  Maximum  Stress  in  the  Patch 

The  maximum  stress  in  the  patch  is  then  determined  as: 

„  I  t  .  1 1*  +  ^2 xPI  |  ^2 xRI  -  P\x 

G„  =  CJV  —  + - + - 


•  Maximum  Adhesive  Strain 

CalcuRep  uses  either  equation  10  or  12.  replacing  a0  with  a0  New  depending  on 
whether  the  adhesive  is  elastic  or  plastic,  to  calculate  the  maximum  adhesive  shear 
strain. 

3.2.2  Limitations  of  CalcuRep 

CalcuRep  has  a  number  of  limitations  on  its  use,  some  of  which  are  listed  below: 

•  Input  Parameters 

•  Material  Data  Library 

•  The  Effect  of  Adhesive  Visco-Plasticity,  see  Chapter  3 

•  Interlaminar  Failure  in  the  Patch,  see  Chapter  3 

•  Aircraft  Operational  Temperature 

•  Patch  Dimension  Design 

•  Unsupported  Single-Sided  Repairs 
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3.2.2.1  Input  Parameters 

CalcuRep  enforces  strict  restraints  on  various  user  input  data.  For  example, 
CalcuRep  only  excepts  parent  structure  thickness  in  the  following  range: 

0.02"  <  t  <  0.1333" 

Most  other  input  parameters  within  the  software  contain  similar  limits,  which  have 
been  provided  to  prevent  the  user  from  placing  an  unrealistic  parameter  value  into 
the  software.  This  unfortunately  can  restrict  rare,  but  realistic  repair  scenarios. 

3.2.2.2  Material  Data  Library 

CalcuRep' s  adhesive  material  data  library  contains  data  for  FM73,  FM-300-2  and  AF- 
163-2  for  one  temperature  case  only.  Since  adhesive  material  properties  are  sensitive 
to  temperature  change,  using  the  adhesive  material  data  library  can  seriously  limit 
repair  design  for  varying  temperature  environments.  In  this  case  the  user  should 
manually  input  adhesive  parameters  for  each  temperature  case  into  the  software  to 
determine  the  largest  resulting  stresses,  strains  and  stress  intensity  factors.  Then 
continue  the  design  process  to  develop  a  repair,  which  can  be  adapted  to  all 
conceivable  temperature  environments. 


3.2.2.3  Aircraft  Operational  Temperature 


As  shown  previously,  CalcuRep  determines  the  thermal  residual  stresses  by 
considering  the  cure  and  operational  thermal  cycles  of  the  repaired  structure.  To 
determine  the  aircraft  operational  thermal  cycle,  the  temperature  at  which  the 
aircraft  operates  is  necessary.  CalcuRep  uses  atmospheric  lapse  rate  equations  to 
determine  this  temperature  at  a  user  specified  altitude  as  follows. 

For  Hoper  <  11000m  then  Toper  =  288.15  -  0.0065 Hoper  (75) 

For  Hoper>11000m  then  Toper=216.65  K 
where  Hoper  is  the  operational  altitude 

The  above  equations  can  be  used  for  operating  temperatures  equal  to  or  below  the 
temperature  at  sea  level.  For  aircraft  that  operate  at  higher  temperatures  (not 
necessarily  related  to  altitude  temperatures)  the  above  equations  are  not  valid  to 
produce  the  desired  operational  temperature.  Producing  above  sea  level 
temperatures  from  equation  67  would  require  an  equivalent  negative  height. 
Unfortunately  CalcuRep  does  not  allow  a  negative  height  to  be  used  within  the 
software  therefore  limiting  its  analysis  to  only  a  select  range  of  aircraft. 
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3.2.2A  Patch  Dimension  Design 

CalcuRep  provides  no  procedure  to  determine  patch  dimensions.  Instead  it  provides 
the  user  with  a  set  of  suggestions  to  reduce  stresses,  strains  and  stress  intensity 
factors  by  changing  patch  material  and  dimensions.  This  procedure  is  adequate  if  the 
user  has  a  preliminary  size  approximation  to  begin  design  with.  On  the  other  hand, 
this  procedure  can  be  difficult  to  begin  when  the  user  has  no  concept  of  an  initial 
patch  size. 

3.2.2.5  Unsupported  Single-Sided  Repairs 

CalcuRep  attempts  to  deal  with  unsupported  single-sided  repair  effects  by 
considering  out  of  plane  bending  (section  3. 2.1 .3).  The  maximum  stress  contribution 
at  the  edge  of  the  patch  is  determined  by  summing  the  corresponding  applied, 
thermal  and  bending  stress  components.  This  calculation  simply  provides  a 
maximum  stress,  not  the  actual  stress  variation  through  the  thickness,  and  provides 
for  a  good  check  to  see  whether  the  material  will  yield. 

Determining  the  stress  under  the  patch,  the  reduced  stress  intensity  factor  and  the 
shear  strains  due  to  out  of  plane  bending  is  far  more  complex'  than  the  process 
explained  above.  Firstly,  the  stress  under  the  repair  in  the  structure  will  vary  through 
the  thickness.  Secondly,  the  reduced  stress  intensity  factor  from  equations  14  &  15  is 
specifically  for  a  uniform  stress  variation  through  the  thickness,  therefore  can't  be 
used  to  predict  results  with  a  stress  variation.  Lastly,  the  maximum  shear  strain  is  not 
given  by  equations  10  &  12  since  out  of  plane  bending  provides  a  larger  contribution 
to  the  no  bending  solution. 

CalcuRep  does  not  deal  will  these  issues  therefore  is  similar  in  process  to  C5033 
which  limits  its  analysis  to  single-sided  supported  repairs. 

3.2.3  Errors  Identified 

A  number  of  discrepancies  have  been  identified  in  CalcuRep's  bonded  repair 
analysis  tools.  These  errors  have  seriously  hindered  reliable  repair  design  by 
CalcuRep,  therefore  rectifying  these  discrepancies  is  fundamental  in  providing  a 
software  package  that  can  be  used  with  confidence.  Errors  in  analysis  tools  used  by 
the  software  have  been  identified  in  the  following  analysis  components: 

•  Out  of  plane  bending 

•  Thermal  residual  stresses  post  repair 

Solving  a  non-linear  differential  equation  derives  the  out  of  plane  bending  analysis 
tools.  The  moment  used  to  form  this  differential  equation  inside  the  repaired  region 
consists  of  a  deflection  contribution  and  a  neutral  axis  offset  contribution.  The 
current  analysis  in  [1]  only  considers  the  deflection  contribution  in  setting  up  this 
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governing  differential  equation,  thus  over  estimates  the  maximum  bending  stress  at 
the  edge  of  the  patch. 

Errors  in  the  thermal  residual  stress  analysis  tools  tend  to  produce  compressive 
stresses  at  the  edge  of  the  repair  in  the  parent  structure,  rather  than  tensile  stresses 
for  above  ambient  cure  temperatures  when  the  plate  CTE  is  higher  than  the  patch 
CTE.  Revisiting  the  derivations  of  these  formulas  as  listed  in  [1],  identified  where 
these  errors  occur.  The  following  sections  provide  the  correct  derivations. 

3.2.3.1  Out  of  Plane  Bending  Discrepancies 


See  Appendix  A1  for  the  correct  derivation.  Equations  58  to  63  should  be  replaced  by 
the  following  corrections.  (Note:  Notation  as  per  section  3.2.1.2) 


Bending  Stiffness 
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Ef 
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E  t 3 

n  _  _fV__ 
12(l -Vj22) 


E  t 3 

^xRLR 


EJr 


Y 2 


-  +  —  -zv 
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•  Neutral  Axis  Offset 


;  _  ) 

2(1  +  5) 

where 

“  V\2RV2\r) 


(77) 

(78) 


•  Maximum  Bending  Stress  at  the  Edge  of  the  Patch 


where 


a 


Bend  max 


Ej_  tanh(qQ7) 

2  <bM  +<b/tanh(aO/) 


(79) 


(80) 


3.2.3.2  Thermal  Residual  Stress  Discrepancies 

In  this  case  the  corrected  equations  are  presented  in  the  Appendix. 
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Enforcing  the  rigid  bond  assumption  produces  the  following  Pi  loads: 
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•  External  Loads  on  the  inclusion  and  surrounding  matrix  post  cure  produces  the 
following  P2  loads: 
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•  The  Inclusion  Experiences  Displacements  due  to  P2  Loads: 
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•  P2  Loads  are  Distributed  into  the  Parent  Structure  and  Patch  in  the  Repaired 
Region 
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Having  the  internal  and  external  loads  within  the  repaired  structure  on  cooling  after 
cure,  stresses  at  various  locations  are  determined  as  follows: 

•  Stress  Under  the  Repair 


P  +P 

(J  —  1*  2  xPI 


(90) 
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Stress  in  Parent  Structure  at  the  Edge  of  the  Patch 


P  +  P 

r\x  T  r2x 

t 


•  Maximum  Stress  in  the  Patch 


R  max 


Pu+P* 


(91) 
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4  Case  Study  of  the  F-111C  Lower  Wing  Skin  Repair 


In  February  of  1994,  a  fatigue  crack  was  discovered  on  an  Australian  F-111C  aircraft, 
tail  number  A8-145,  located  in  the  lower  wing  skin  at  Forward  Auxiliary  Spar  Station 
281.28.  Failure  to  repair  such  a  flaw  on  primary  structure  would  have  seen  the  wing 
discarded  from  service.  Rather  than  discarding  the  wing,  it  was  decided  that  the 
wing  would  be  repaired.  Two  repair  options  were  available  to  the  RAAF,  a 
mechanically  fastened  repair  or  an  adhesively  bonded  repair^  system.  Using  a 
mechanically  fastened  repair  would  have  required  additional  hole  locations  for 
fasteners.  Additional  holes  induce  more  stress  concentrations  thus  increase  the 
probability  of  initiating  fatigue  problems.  Also  mechanically  fastened  repair  transfer 
load  at  the  fastener  locations  only  and  some  fasteners  take  up  more  load  that  others 
which  degrades  its  efficiency  in  load  carrying  capability.  On  the  other  hand,  adhesive 
bonded  repairs  provide  efficient  load  transfer  since  the  adhesive  is  active  throughout 
the  repaired  region.  Therefore  the  RAAF  decided  to  use  an  adhesive  bonded  repair. 

The  repair  was  designed  using  the  RAAF  Engineering  Standard  C5033  [3,  4,  5].  The 
repair  was  designed  to  be  applicable  to  the  following  temperature  cases;  40,  -25  and 
167  degrees  Fahrenheit  [16].  Repair  specifications  and  processes  were  listed  as  per 
references  [4,  5]  to  ensure  the  high  quality  of  repair  installation.  NDI  inspection 
found  no  crack  growth  after  850  Airframe  Flight  Hours  prior  to  the  wing  being 
retired  from  service  as  a  fatigue  test  and  tear-down  item. 

A  comprehensive  repair  substantiation  program  was  conducted  at  the  DSTO  and  is 
•listed  in  references  [6]  and  [7],  Repair  efficiency  was  demonstrated  by  reference  [6] 
and  [7]  and  indicated  that  the  C5033  repair  methodology  produced  a  good  and 
feasible  repair.  The  particular  parameter  of  interest  to  this  report  from  reference  [6]  is 
the  reduced  stress  intensity  factor.  The  Finite  Element  results  are  summarised  in  the 
'results'  section  below. 

Having  a  real-life  case  study  from  the  preliminary  design  stage  to  the  actual 
implementation  of  the  repair  is  beneficial  in  providing  a  baseline  standard  to 
compare  analytical  results.  This  design  case  was  used  to  compare  the  current  RAAF 
engineering  standard  C5033  with  the  original  CalcuRep  and  error  free  CalcuRep 
codes  in  the  following  sections. 
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Due  to  the  limitations  identified  in  section  3.2.2,  a  Mathcad  module  containing  the 
exact  coding  of  the  CalcuRep  software  was  developed.  This  enabled  most  limitations 
to  be  lifted  and  provided  the  flexibility  of  material  data  entry  for  the  module. 
Likewise,  the  design  tools  used  in  C5033  [3]  where  entered  into  a  Mathcad  module. 
The  benefit  of  this  process  allows  the  user  to  view  the  mathematical  expressions  used 
in  each  module  and  allows  changes  to  be  made  where  deemed  appropriate.  This  also 
allows  programmers  to  view  the  values  of  intermediate  variables  in  order  to  verify 
software  against  the  modules.  Two  example  cases,  which  will  be  discussed  later, 
have  been  provided  in  appendix  A4  and  A5,  listing  the  current  RAAF  C5033 
standard  and  the  error  free  CalcuRep  software. 

4.1  Material  Input  Data 
Skin  Data  [15] 

Skin  Type:  2024-T851 

Skin  Thickness:  0.14in 

Yield  Strength:  59ksi 

Ultimate  Strength:  65ksi  ' 

Fracture  Toughness:  42ksiVin 

Elastic  Modulus:  10.5xl03ksi 

Poisson's  Ratio:  0.31 

Coefficient  of  Thermal  Expansion  12.6xl0-6in/in/degrees  F 

Adhesive  Data  [16] 

Adhesive  Type:  FM-73  (0.06  psf) 

•Adhesive  Thickness:  0.01  in 

Cure  Temperature:  180  Degrees  F 

Poisson's  Ratio:  0.33 

Material  property  variation  per  temperature  environment  as  follows: 
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Coefficient  of  Thermal  Expansion  (Longitudinal  Direction):  2.3x1c)-6  /F 
Coefficient  of  Thermal  Expansion  (Transverse  Direction):  10.6xl0'6  /F 

Major  Poisson's  Ratio:  0.21 


Following  properties  [6]  per  the  following  ply  lay-up  [0,0,45,-45,0, 0,0]s 


Elastic  Modulus  Longitudinal:  22.6xl03ksi 

Elastic  Modulus  Transverse:  4.319xl03ksi 

Coefficient  of  Thermal  Expansion  (Longitudinal  Direction) [20]:  3.486xl0‘6  /F 

Coefficient  of  Thermal  Expansion  (Transverse  Direction)  [20]:  9.41x1c)-6  /F 

Major  Poisson's  ratio:  0.575 

Patch  Thickness:  0.0728  in 

Loading 

Applied  axial  stress:  37.1ksi 

Patch  Dimensions 
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C5033  provides  a  method  for  determining  the  patch  dimensions  as  discussed 
previously.  Therefore  determining  the  required  patch  dimensions  in  C5033  and 
placing  them  into  CalcuRep  will  allow  a  feasible  comparison  to  be  made.  The 
dimensions  calculated  by  C5033  are  summarised  below. 


Temperature(F) 

-40 

-25 

75 

167 

Patch 

Length(in) 

11.34 

11.44 

12.53 

20.74 

Patch 

Width(in) 

12.7 

12.7 

12.7 

12.7" 

Table  2 


Establishing  a  CalcuRep's  Repair  Geometry  Similar  to  C5033 

CalcuRep  and  C5033  do  not  run  with  the  same  assumptions.  CalcuRep's  repair 
geometry  consists  of  a  finite  plate  with  a  finite  patch;  the  plate  being  fully  restrained 
at  its  boundary.  C5033  uses  an  infinite  plate  and  a  finite  patch  with  the  plate  being 
unrestrained  at  its  boundary.  C3055  attempts  to  restraint  the  parent  structure 
boundary  by  using  an  effective  coefficient  of  thermal  expansion,  equation  (5),  which 
is  derived  assuming  a  fixed  edge  condition,  see  [20]  for  details.  To  set  CalcuRep  to 
run  similarly  to  C5033,  its  plate  dimension  must  be  increased  to  infinity.  Therefore  Rx 
and  Ry,  being  the  outer  plate  radius  in  the  'x'  and  'y'  directions  respectively,  must  be 
set  to  infinity. 

Case  Study  Run  for  C5033  and  CalcuRep 

Having  the  appropriate  material  data,  the  next  step  is  to  decide  what  parameters  to 
compare.  Common  parameters  computed  by  C5033  and  CalcuRep  are  (the  following 
notation  will  be  used  in  Tables  3-6  in  the  results  section): 

•  cta  Stress  under  the  repair  in  the  parent  structure  due  to  the  applied  stress 

•  ctt  Stress  under  the  repair  in  the  parent  structure  due  to  thermal  residual  stresses 

•  Go  Total  stress  under  the  repair  in  the  parent  structure  (applied+thermal) 

•  ctia  Stress  in  the  parent  structure  at  the  edge  of  the  patch  due  to  applied  stresses 

•  on  Stress  in  the  parent  structure  at  the  edge  of  the  patch  due  to  thermal  residual 
stresses 

•  gi  Total  stress  in  the  parent  structure  at  the  edge  of  the  repair 

•  yMax  Maximum  shear  strain  in  the  adhesive 

•  Kao  Upper  bound  or  limiting  stress  intensity  factor  post  repair 

•  Op  Max  Maximum  patch  stress  assuming  the  presence  of  a  crack 
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Running  the  F-111C  Boron/Epoxy  Repair  Case 

Using  the  material  data  listed  previously,  the  parameters  above,  as  calculated  by 
C5033,  the  original  CalcuRep  and  the  error  free  CalcuRep,  can  be  determined  at  each 
operating  temperature  environment,  ie  -40F,  -25F,  75F  and  167F.  The  parameters  can 
then  be  compared  for  each  case  to  determine  whether  the  CalcuRep  approach  is  valid 
against  the  already  verified  C5033  [6,  7]  approach.  An  example  of  each  code  in 
MathCad  form  has  been  provided  in  as  part  of  this  report  and  consists  of  the  -25F 
case  for  C5033  and  the  error  free  CalcuRep  code.  The  remaining  runs  are  summarised 
in  the  following  section. 

4.2  Analytical  Comparison 


4.2.1  Casel:  -40  F 


Parameter 

C5033 

CalcuRep 

(Original) 

CalcuRep 

(Corrected) 

cta  (ksi) 

21.01 

20.58 

20.58 

gt  (ksi) 

7.67 

5.80 

13.56 

a0  (ksi) 

28.68 

26.38 

34.14 

oia  (ksi) 

44.52 

44.50 

44.50 

an  (ksi) 

16.26 

-4.09 

22.07 

cji  (ksi) 

60.78 

40.41 

>66.57 

YMax 

0.198 

0.174 

0.262 

Koo  (ksiVin) 

23.02 

20.97 

.  28.15 

Op  Max  (ksi) 

116.90 

72.36 

128.02 

i 

:able  3 

4.2.2  Casel:  -25  F 


Parameter 

C5033 

CalcuRep 

(Original) 

CalcuRep 

(Corrected) 

ga  (ksi) 

21.01 

20.61 

20.61 

gt  (ksi) 

7.00 

5.15 

12.08 

Go  (ksi) 

28.01 

25.76 

32.69 

44.52 

44.56 

44.56 

14.83 

-3.65 

19.75 

gi  (ksi) 

59.35 

40.91 

64.31 

YMax 

0.172 

0.25 

22.48 

20.48 

26.91 

Op  Max  (ksi) 

114.10 

73.92 

123.67 

i 

fable  4 
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4.2.3  Case3:  75  F 


Parameter 

C5033 

CalcuRep 

(Original) 

CalcuRep 

(Corrected) 

cta  (ksi) 

21.01 

20.92 

20.92 

cnr  (ksi) 

2.48 

2.23 

2.22 

a0  (ksi) 

23.49 

23.15 

23.14 

44.52 

45.19  • 

45.19 

■KSXtm 

5.26 

-1.64 

4.29 

gi  (ksi) 

49.78 

49.48 

YMax 

0.179 

0.175 

0.175 

ICo  (ksiVin) 

19.33 

19.00 

18.99 

Op  Max  (ksi) 

95.72 

81.70 

95.16 

Table  5 


4.2.4  Casel:  167  F 


Parameter 

C5033 

CalcuRep 

(Original) 

CalcuRep 

(Corrected) 

oa  (ksi) 

21.01 

22.91 

<  22.91 

GT  (ksi) 

-1.67 

6.18 

-6.85 

Go  (ksi) 

19.33 

29.10 

,16.06 

oia  (ksi) 

44.52 

49.38 

49.38 

oit  (ksi) 

-3.55 

-4.35 

-9.94 

gi  (ksi) 

40.97 

45.03 

39.44 

YMax 

0.597 

0.206 

Koo  (ksiVin) 

21.41 

36.59 

CTp  Max  (ksi) 

78.79 

80.88 

75.85 

Table  6 


4.2.5  FE  results  from  reference  [6] 


Temperature 

-40F 

-25F 

167F 

fvo(ksiVin) 

24.25 

24.25 

22.82 

1 

'able  7 

5  DISCUSSION 

The  RAAF  Engineering  Standard  C5033  and  CalcuRep  version  1.0  are  founded  on 
almost  identical  mathematical  analysis  tools,  with  exception  to  their  methods  used  to 
calculate  thermal  residual  stresses.  Other  differences  between  the  two  approaches 
have  been  discussed  in  detail  in  Sections  3.1  and  3.2  and  mainly  consist  of 
assumption  differences.  Regardless  of  these  differences,  the  results  produced  by  the 
two  approaches  are  expected  to  correlate  well. 
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Finite  Element  Analysis  and  testing  have  demonstrated  validation  of  the  RAAF 
standard  C5033  via  reference  [6]  and  [7],  Therefore  the  C5033  approach  has  been 
used  to  provide  the  baseline  comparison  data  listed  in  tables  3  to  6.  A  good 
correlation  between  C5033  and  CalcuRep  will  directly  indicate  the  validity  and 
accuracy  of  the  CalcuRep  approach. 

To  isolate  specifically  where  and  why  discrepancies  occur  between  the  three 
Mathcad  codes,  i.e  C5033,  CalcuRep  (original)  and  CalcuRep  (corrected),  a  discussion 
of  each  variable  listed  in  Tables  3  to  6  will  be  provided  as  follows; 

•  Stress  Under  the  Repair  in  the  Parent  Structure  due  to  Applied  Stresses  (a a) 

Reference  [11]  derives  the  expression,  as  used  by  CalcuRep,  for  the  stress  in  the 
parent  structure  under  the  repair  in  a  generic  form.  The  formulation  is  given  by 
equation  43  and  assumes  the  parent  structure  to  be  isotropic  and  the  patch  to  be 
orthotropic,  both  materials  having  different  Poisson's  ratios.  Using  equation  43, 
C5033  assumes  that  both  the  parent  structure  and  the  patch  are  isotropic  and  consist 
of  the  same  Poisson's  ratio,  simplifying  the  generic  form  of  equation  43  down  to 
equation  8.  Therefore  C5033  uses  a  simplified  form  of  the  expression  used  by 
CalcuRep  (Equation  43),  and  explains  the  good  correlation  of  results  between  the  two 
approaches  as  shown  in  Tables  3  to  6. 

•  Stress  Under  the  Repair  in  the  Parent  Structure  due  to  Thermal  Residual 
Stresses  (ctt) 

C5033  and  CalcuRep  (corrected)  followed  the  expected  trends  for  all  temperature 
cases  as  shown  in  tables  3  to  6  above.  (Note:  negative  residual  stresses  in  the  parent 
structure  under  the  repair  indicate  compressive  stresses.  Compressive  stresses  result 
in  the  167F  (Table  6)  temperature  case  since  the  operational  thermal  cycle  is  more 
dominant  in  compression  than  the  cure  cycle  is  in  tensile.)  The  most  startling  result 
(Table  6)  is  produced  by  the  CalcuRep  (original)  code,  which  indicates  a  tensile 
residual  stress  for  the  167  degree  Fahrenheit  case.  This  problem  is  predominantly  due 
to  the  errors  identified  and  rectified  (CalcuRep  corrected)  in  sections  3.2.3.2. 

Although  C5033  and  CalcuRep  (corrected)  correlated  well  in  trend,  any  discrepancies 
between  them  can  be  explained  by  the  assumptions  used  in  their  respective  thermal 
residual  stress  models.  Consider  the  following  assumptions: 

•  C5033  Thermal  Residual  Stress  Model 

One-dimensional  thermal  residual  stress  model 

Infinitely  large  structure 

Unrestrained  at  the  structures  boundary 

Uses  an  effective  coefficient  of  thermal  expansion  which  reduces  the  actual  thermal 
expansion  coefficient  of  the  parent  structure,  due  to  cooler  surrounding  material 
resisting  the  expansion  of  the  hotter  material  near  the  heated  region.  Equation  5  is 
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used  to  calculate  the  effective  coefficient  of  thermal  expansion  and  is  based  on  a  fixed 
edge  infinite  circular  plate  [20], 

•  CalcuRep  (corrected) 

Two-dimensional  thermal  residual  stress  model 

Finite  patch  and  finite  parent  structure 

Fully  restrained  at  the  parent  structure's  boundary 

Uses  a  more  complex  expression  for  the  effective  coefficient  of  thermal  expansion  of 
the  parent  structure  accounting  for  a  finite  width  and  the  effects  of  heat  transfer  from 
the  warmer  material  to  the  cooler  material. 

Assumes  the  parent  structure  to  be  isotropic  and  the  patch  to  be  orthotropic. 

Obviously  these  assumptions  can  explain  any  discrepancies,  but  regardless  of  the 
differences,  the  results  correlated  well  in  trend. 

Consider  the  new  approach  developed  in  reference  [8]  for  single-sided  supported 
repairs,  see  Appendix.  The  thermal  residual  stress  model  used  in  the  Appendix  was 
validated  by  a  fully  three-dimensional  finite  element  analysis  [22].  Running  the  -25 
degree  Fahrenheit  case  for  the  F-111C  Case  Study  has  produced  n  good  agreement 
with  CalcuRep  for  thermal  residual  stress  under  the  repair.  When  the  edge  of  the 
structure  is  fully  restrained  at  the  boundary,  the  MathCad  module  [8]  can  produce 
results  which  can  differ  by  a  factor  of  two  (higher)  than  C5033  for  this  stress 
component.  This  suggests  that  the  thermal  residual  stress  calculated  by  the  C5033 
approach  may  not  be  conservative. 

•  Go  Total  Stress  Under  the  Repair  in  the  Parent  Structure  (applied+thermal) 

The  total  stress  under  the  repair  in  the  parent  structure  is  a  summation  of  the  applied 
and  thermal  residual  stress  components.  Therefore  errors  due  to  these  two 
components  will  be  accumulated  at  this  point.  As  expected.  Tables  3  to  6  show  that 
C5033  and  CalcuRep  (correcteded)  follow  the  expected  trends  for  all  temperature 
cases.  CalcuRep  (original)  produced  an  odd  result  for  the  167  Fahrenheit  case, 
indicating  errors  in  analysis  tools  as  previously  discussed. 

•  gia  Stress  in  the  Parent  Structure  at  the  Edge  of  the  Patch  due  to  Applied 
Stresses 

Reference  [11]  derives  expressions,  as  used  by  CalcuRep,  to  calculate  the  stress 
attracted  into  the  repaired  region,  equation  (34),  in  a  generic  form.  The  formulation 
assumes  the  parent  structure  to  be  isotropic  and  the  patch  to  be  orthotropic,  both 
having  different  Poisson's  ratios.  The  stress  attracted  into  the  repaired  region  is  then 
related  to  the  stress  in  the  parent  structure  at  the  edge  of  the  patch  via  equations  41 
and  42.  C5033  does  not  calculate  the  stress  attraction  factor  but  rather  uses  a  constant 
value  of  1.2. 
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All  three  Mathcad  runs  produced  consistent  results  for  the  stress  in  the  parent 
structure  at  the  edge  of  the  patch,  demonstrating  that  C5033's  stress  attraction  factor 
of  1.2  is  a  valid  approximation. 

•  arr  Stress  in  the  Parent  Structure  at  the  Edge  of  the  Patch  due  to  Thermal 
Residual  Stresses 

Consistent  results  were  produced  between  C5033  and  CalcuRep  (corrected))  for  all 
temperature  cases.  On  the  other  hand,  CalcuRep  (Original)  calculated  compressive 
thermal  residual  stresses  instead  of  tensile  for  the  -40F,  -25F  and  75F  temperature 
cases. 

Errors  of  this  natural  are  mostly  due  to  incorrect  algebraic  manipulations  (Shown  in 
section  3.2.3.2)  in  the  analysis  tools  and  one  major  programming  error.  The 
programming  error  occurs  in  the  operational  thermal  cycle  procedure.  The 
temperature  change  used  in  the  original  software  is  as  follows:  AT=T0p-Tamb,  where 
Top  and  Tamb  are  the  operational  and  ambient  temperature,  respectively.  The  aircraft 
is  travelling  from  ambient  temperature  to  its  operational  temperature,  therefore  the 
change  in  temperature  should  be  given  as:  AT=  Tamb-  T0p  using  CalcuRep's  sign 
convention. 

•  ai  Total  Stress  in  the  Parent  Structure  at  the  Edge  of  the  Repair 

v 

The  total  stress  in  the  parent  structure  at  the  edge  of  the  patch  is  a  summation  of 
applied  and  thermal  residual  stresses  at  this  location.  As  a  result  of  previous 
calculations,  C5033  and  CalcuRep  (corrected)  correlate  reasonably  and  CalcuRep 
(Original)  correlated  inconsistently  with  the  trends  shown  in  tables  3  to  6. 

•  YMax  Maximum  Shear  Strain  in  the  Adhesive 

C5033,  CalcuRep  (Original)  and  CalcuRep  (corrected)  use  identical  expressions  for 
the  maximum  shear  strain  in  the  adhesive,  equations  (10)  and  (12).  Equations  (10) 
and  (12)  are  a  function  of  the  stress  in  the  parent  structure  under  the  repair,  therefore 
discrepancies  in  the  maximum  shear  strain  between  the  codes  can  be  explained  by 
the  discrepancies  between  the  stress  under  the  patch.  As  a  result,  C5033  and 
CalcuRep  (corrected)  have  produced  consistent  results  while  CalcuRep  (original)  has 
produced  inconsistent  results. 

•  ICo  Upper  Bound  or  Limiting  Stress  Intensity  Factor  Post  Repair 

C5033,  CalcuRep  (Original)  and  CalcuRep  (Corrected)  use  identical  expressions  for 
the  upper  bound  stress  intensity  factor  post  repair  (equations  14  and  15).  Equations 
14  and  15  are  again  a  function  of  the  stress  in  the  parent  structure  under  the  repair, 
therefore  discrepancies  in  the  upper  bound  stress  intensity  factor  between  the  codes 
can  be  explained  by  discrepancies  in  the  stress  under  the  repair.  As  a  result,  C5033 
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and  CalcuRep  (corrected)  have  produced  consistent  results  and  CalcuRep  (original) 
has  produced  inconsistent  results. 

Comparing  the  results  with  the  Finite  Element  Analysis  in  Table  7  shows  good 
agreement  between  C5033  and  CalcuRep  (corrected).  CalcuRep  (original)  has 
produced  an  odd  result  for  the  167F  temperature  case,  which  is  due  to  the  over 
estimate  of  the  stress  under  the  patch  in  the  parent  structure. 


Although  comparable,  some  differences  have  occurred  between  the  FE,  C5033  and 
CalcuRep  (corrected).  Firstly,  when  the  Finite  Element  Analysis  was  performed,  it 
was  done  considering  the  cure  and  operation  cycle  as  one  cycle,  (ie  AT=TC-T0p).  The 
current  C5033  standard  now  includes  two  independent  thermal  cycles,  a  cure  and  an 
operational  cycle,  see  section  3.1.3.  Secondly,  the  Finite  Element  Analysis  did  not 
assume  the  adhesive  to  be  elastic/ perfectly  plastic.  In  stead  it  treated  the  adhesive  as 
follows: 


Actual  Curve 


Assumed  Curve 


C5033  and  CalcuRep  have  treated  the  adhesive  as  elastic/ perfectly  plastic  as  per  the 
C5033  standard,  see  Figure  11.  Lastly,  the  stress  attraction  factor  used  in  reference  [6] 
for  the  applied  stress  component  was  1.44,  which  was  derived  from  considerations 
out  of  reference  [14].  This  value  could  not  be  used  in  the  above  analysis  since  it 
would  have  required  modifying  the  CalcuRep  code.  As  a  result,  CalcuRep  was 
allowed  to  calculate  the  stress  attraction  factor,  which  resulted  in  a  value 
approximately  equal  to  1.2,  similar  to  that  used  by  C5033.  Regardless  of  these  issues, 
the  results  are  consistent. 


•  Op  Max  Maximum  Patch  Stress  Assuming  the  Presence  of  a  Crack 

C5033  and  CalcuRep  (corrected)  have  produced  consistent  results  for  the  maximum 
patch  stress  for  all  temperature  cases.  CalcuRep  (original),  on  the  other  hand  has 
correlated  poorly  with  the  other  two  codes.  This  is  directly  due  to  the  incorrect 
direction  calculation  of  the  stress  in  the  parent  structure  at  the  edge  of  the  patch  for 
the  thermal  stress  component,  and  the  errors  listed  in  section  3.2.I.2. 
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6  CONCLUSIONS  AND  RECOMMENDATIONS 


The  verification  of  mathematical  analysis  tools  used  in  bonded  repair  technology  is 
essential  in  providing  the  means  for  good  repair  design.  Designing  a  repair  with 
unverified  design  tools  or  software  can  lead  to  an  unreliable  repair.  This  report  has 
provided  a  means  of  verification  for  the  USAF  bonded  repair  design  software 
package  called  "CalcuRep  Version  1.0". 

During  initial  runs  with  the  CalcuRep  software,  a  number  of  limitations  were 
identified  as  shown  in  section  3.2.2.  These  limitations  led  to  the  development  of  a 
Mathcad  module  containing  the  exact  coding  used  in  the  original  software.  This 
enabled  most  limitations  to  be  lifted  and  provided  the  flexibility  of  material  data 
entry  for  the  module. 

A  Boron/Epoxy  composite  repair,  which  was  implemented  on  a  RAAF  F-111C 
aircraft  [16]  in  1994,  provided  the  means  for  a  verified  case  study  to  compare 
analytical  results.  This  repair  was  designed  and  specified  using  the  RAAF 
Engineering  Standard  C5033,  which  was  verified  by  Finite  Element  Analysis  and 
testing  [6  &  7].  Repair  efficiency  was  demonstrated  at  three  temperature  conditions  (- 
40F,  -25F  and  167F).  An  additional  temperature  condition  is  provided  in  this  report  at 
75  Fahrenheit  to  view  post  repair  stresses  at  ambient  temperature.  The  C5033 
approach  was  also  entered  into  a  Mathcad  module,  see  Appendix.  The  RAAF 
Engineering  Standard  C5033  in  conjunction  with  the  F-lllC  Boron/ Epoxy 
Composite  repair  case  study  was  used  to  produce  baseline  results  for  stresses,  stress 
intensity  factors  and  shear  strain  post  repair. 

Running  the  F-lllC  Boron/Epoxy  composite  case  study  through  the  Mathcad 
modules  for  CalcuRep  and  C5033  produced  some  strange  findings.  Inconsistent 
results  were  produced  by  CalcuRep,  especially  pertaining  to  the  stress  in  the  parent 
structure  at  the  edge  of  the  patch.  The  stresses  produced  for  this  case  where  in  the 
wrong  directions.  This  led  to  a  further  examination  of  the  original  CalcuRep  code 
and  the  analytical  theory  it  is  based  on. 

On  further  investigation  of  CalcuRep' s  analysis  tools  and  programming  structure, 
many  errors  were  identified.  Discrepancies  were  found  in  reference  [1]  in  the  out  of 
plane  bending  stress  and  thermal  residual  stress  formulations.  Revisiting  the  analysis 
demonstrated  where  and  why  errors  occurred  and  rectified  solutions  to  the  problems 
were  offered  in  sections  3.2.3.  The  software  was  also  found  to  use  the  incorrect 
temperature  change  for  an  aircraft  operation  thermal  cycle.  This  problem  was 
identified  as  a  major  programming  error,  which  led  to  an  incorrect  sign  usage  by 
CalcuRep.  As  a  result,  another  Mathcad  module  was  developed  containing 
corrections  to  the  original  CalcuRep  module:  noted  as  CalcuRep  (corrected),  see 
Appendix. 
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Re-running  the  corrected  CalcuRep  version  produce  good  results  against  C5033. 

These  results  were  consistent  in  trend  and  good  in  magnitude  and  showed  that  by 

fixing  the  errors  in  CalcuRep,  the  software  became  consistent  with  C5033. 

A  number  of  recommendations  for  CalcuRep  were  established  throughout  this  report 

and  are  summarised  as  follows: 

•  Make  corrections  to  the  original  CalcuRep  code  as  shown  in  sections  3.2.3 

•  Make  a  correction  to  the  programming  error  shown  in  section  5. 

•  Remove  the  restraints  on  user  material  input  and  replace  it  with  a  warning  flag  in 
the  software.  An  example  is  a  small  flashing  window  which  warns  the  user  that 
the  value  they  have  entered  my  not  be  applicable  to  the  software. 

•  Allow  CalcuRep  to  change  the  adhesive  properties  with  temperature. 

•  Allow  the  user  to  input  the  aircraft  operational  temperature  directly  instead  of 
using  the  aircraft  operational  altitude. 

•  Provide  CalcuRep  with  a  method  to  produce  preliminary  patch  dimensions. 

•  Expand  CalcuRep  to  deal  with  unsupported  single-sided  repairs  and  a  range  of 
other  repair  geometry  through  the  use  of  more  advanced  bonded  repair  design 
tools. 

•  Compare  the  thermal  residual  stress  formulation  derived  in  reference  [22]  with 
the  revised  CalcuRep  analysis  from  section  3.2.3.2.  Reference  [22]  correlated 
excellently  with  Finite  Element  Analysis  methods,  therefore  a  good  correlation 
between  section  3.2.3.2  and  reference  [22]  will  indicate  the  validity  of  CalcuRep' s 
revised  thermal  residual  stress  equations. 

•  Check  if  equations  44  and  45  have  been  derived  correctly  (Note:  These  equations 
are  directly  from  [1]).  Time  constraints  did  not  allow  these  equations  to  be  re¬ 
derived  and  checked,  therefore  they  may  not  be  consistent  with  the 
corresponding  effective  coefficient  of  thermal  expansion  equations  from  [22], 
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Appendix  A:  Corrections  to  Equation  Discrepancies 


Corrections  to  Out  of  Plane  Bending  Equations 

The  following  derivation  outlines  the  problems  that  exist  in  reference  [1]  for  the  out 
of  plane  bending  analysis.  These  problems  cause  excessively  high  bending  stresses  at 
the  edge  of  the  patch.  The  following  derivation  provides  the  corrections  to  the 
analysis  (pp  D1-D7  of  [1])  and  indicates  where  these  discrepancies  occur. 

(Note:  All  constants  and  notation  listed  in  the  following  are  as  per  section  3.2.1.3 
unless  otherwise  stated) 

Derivation: 


The  geometry  of  the  repaired  structure  is  shown  in  Figure  26  below.  The  V  axis  is 
placed  about  the  plate  centroid  for  convenience. 
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Figure  26 

It  is  assumed  that  the  patch  covers  the  entire  width  of  the  parent  structure. 

Bending  Stiffness  inside  and  outside  the  inclusion 

(Note:  The  inclusion  section  refers  to  the  repair  and  plate  inside  the  repaired  region 
•  as  one  entity.) 

Inside  The  Inclusion 

The  second  moment  of  area  of  the  plate  and  repair  about  the  neutral  axis  of  the 
inclusion  is  as  follows: 
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g  _  ExRtR  ^12 ) 
Ext^-vnRvllR) 


(Al-10) 


Discrepancies  also  exist  between  Al-10  above  and  the  equivalent  equation  presented 
in  reference  [1].  Differences  occur  since  plane-strain  elastic  moduli  are  used  above 
while  plane-stress  elastic  moduli  are  used  in  reference  [1].  As  shown  in  reference  [21] 
page  2764,  plane-strain  conditions  prevail  rather  than  plane-stress  conditions. 

This  completes  the  analysis  of  all  the  required  material  parameters. 

Secondary  Bending  Formulation 


Consider  the  region  outside  the  inclusion 

Outside  the  inclusion  the  following  bending  moment  exists: 


Mp  =  Ptw{x) 


(Al-11) 


Where,  Mp  moment  outside  the  inclusion  as  a  function  of  V 
P  Applied  stress 
w(x)  Deflection  outside  the  inclusion 

For  a  one-dimensional  system,  the  deflection  is  given  as  follows: 

d2w  _  Mp 
dx2  Dp 

Solving  for  w(x)  gives  the  following: 
wix)=Aex0M  +Be~x*“ 


(Al-12) 


(Al-13) 
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where  A,  B  are  unknown  constants.  The  boundary  conditions  outside  the  inclusion 
are  as  follows: 


w(oo)=0 
w(a)=Mi/ (Pt) 


(Al-14) 

(Al-15) 


Replacing  'x'  in  equation  Al-12  with  'x-a'  ,  see  Figure  27,  and  solving  for  the 
boundary  conditions  we  find  that 

W=Mi/Pt  e(a-x)PM  (Al-16) 


This  furnishes  the  solution  of  the  deflection  outside  the  inclusion. 

Consider  inside  the  inclusion 

Inside  the  inclusion  the  following  bending  moment  exists: 


M ,  =  Ptw!  (x)  +  Pt  zx  (Al-17) 

Equation  Al-16  above  considers  the  bending  moment  due  to  deflection  and  the 
neutral  axis  offset.  The  bending  moment  due  to  the  neutral  axis  offset  is  not 
considered  in  the  moment  relation  used  in  reference  [1], 


For  a  one-dimensional  system,  the  deflection  is  given  as  follows: 

d\>L_ML 

~chF~~D^ 


(Al-18) 


Substituting  equation  Al-17  into  Al-18  above  and  solving  for  wi(x)  gives  the 
following: 


w,{x)  -  A je*0'  Jt-Bze~x^‘  -zx 
where 

o;  = 

where  are  A2,  B2  constants  to  be  determined 

The  boundary  conditions  inside  the  inclusion  are  as  follows: 

~(x  =  0)  =  0 ,  wi(a)=M2/(Pt) 
ax 


(Al-19) 


(Al-20) 
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Note:  Reference  [1]  has  an  error  in  the  stated  boundary  conditions,  but  has  used  the 
correct  boundary  condition  in  the  analysis. 

Using  the  boundary  conditions  gives  A2  =  B2  and 


A2=- 

2  2 


M, 


V 


+  Z„ 


V 


Pt  jycosh^O;) 


, ,  (m1 


+  z 
Pt  x 


/^cosh(xO/)^ 

cosh(a<t>7) 


~Zx 


This  furnishes  the  solution  of  the  deflection  inside  the  inclusion. 


(Al-21) 

(Al-22) 


The  bending  moments  Mi  and  M2  can  be  determined  from  considering  the  following 
compatibility  requirements. 


Wj(a)-  w(a)  =  0 
dx  dx 


From  these  equations  we  find  that  Mi=  M2,  and 


(Al-23) 


M\  =  M2 

1^2  _  ~  zX  ®/ tanh(a<l)/) 
Pt  +  O /  tanhfaO/) 


(Al-24) 


The  maximum  bending  stress  in  the  parent  structure  at  the  edge  of  the  repair  is  given 
by  the  following: 


a 


Bend  max 


Ext  zx  tanh(<3<I)/) 
2  +0/tanh(a<I>/) 


(Al-25) 


A.l.  Corrections  to  Thermal  Residual  Stress  Equations 

Errors  in  the  thermal  residual  stress  derivation  in  reference  [1]  are  mostly  due  to 
incorrect  algebraic  manipulation  throughout  the  entire  analysis.  The  following 
provides  the  corrections  to  the  analysis  (pp  C12-C20  of  [1])  and  indicates  where 
discrepancies  occur  in  the  original  analysis. 
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The  geometry  of  the  repaired  structure  is  as  shown  in  Figure  19.  The  repaired 
structure  is  divided  into  the  following  elements: 

Repair 

Reinforced  Plate  Region 
(Region  D)  ^ 

Figure  28 . 


Inclusion 


Surrounding  Matrix  (Region  M) 


Figure  29 


The  patch  is  assumed  to  be  rigidly  bonded  to  the  parent  structure  after  adhesive  cure 
has  occurred.  The  structure  is  then  cooled  down  to  ambient  temperature  resulting  in 
compressive  and  tensile  residual  stresses  in  the  patch  and  under  the  patch 
respectively.  In  order  to  ensure  the  rigid  bond  assumption  complies,  internal  loads 
are  applied  to  the  patch  and  region  D,  noted  as  Pi  loads,  as  shown  in  Fgure  20.  On 
cooling,  region  D  and  the  patch  want  to  experience  the  following  displacements. 


Region  D 

x-y  direction 

8  =  -aP(eff)xATa 

5  =  -aP(eff)yA  Tb 

(Al-26) 

Patch 

x-y  direction 

8  =  -axRATa 

8  =  -ayRATb 

(Al-27) 

The  corresponding  strains  produced  by  the  Pi  loads  (Figure  20)  are  as  follows: 
Region  D 
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Consider  the  Inclusion  and  Surrounding  Matrix  (M) 


The  internal  Pi  loads  in  region  D  and  the  patch  are  now  applied  to  the  inclusion.  This 
leads  to  the  following  two  expressions  for  the  displacement  the  inclusion  experiences 
due  to  enforcing  the  rigid  bond  assumption  on  cooling  post  cure. 


Inclusion  x-y  directions 


8  =  -aI(eff)xATa  +  J  —  -  ^L 
W  {EJ  EyIt 

8  =  -aI(fr)ATb  +  b\^- 

I(eff)y  EyIt  EJ 


(A1  -  39a) 


(A1  -  39b) 


where  oti(eff)x  &  ai(eff)y  are  the  CTE  of  the  inclusion  in  the  'x'  and  'y'  directions 
respectively 

The  derivation  listed  in  [1]  is  the  similar  to  except  that  it  does  not  consider  the  effect 
the  inclusion  as  one  entity  has  on  the  coefficient  of  thermal  expansion.  The  inclusion 
CTE  is  taken  from  [22]  and  is  as  follows: 
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(Al-41) 


Aj=A  +  ^A]t  (Al-42) 

Cj  =  A;1  (Al-43) 
where 

r  i  —  v .  i 


(Al-44) 


=  C,A 


+—CjAr  x 
t  1  R  av 


(Al-45) 


Here,  A,  Ar,  Ai 


Plane  stress  Elastic  Properties  of  the  plate,  repair  and 
inclusion  respectively 
Elastic  Modulus  of  inclusion  (x-direction) 
Elastic  Modulus  of  inclusion  (y-direction) 
Inclusion  Poisson's  ratio  (xy-direction) 
Inclusion  Poisson's  ratio  (yx-direction) 
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•  Surrounding  Matrix  (M) 

As  the  inclusion  experiences  a  displacement  due  to  internal  loads,  the  surrounding 
matrix  experiences  a  displacement  due  to  thermal  stress  as  follows: 


x-y  directions 


8  =  -a-P(eff>ATa 

6  =  -cXp^yATb 


(Al-46) 


The  second  step  is  to  consider  external  loads  applied  to  the  edges  of  the  inclusion 
and  matrix  due  to  cooling  post  cure,  see  Figure  21.  The  following  sections  detail  the 
inclusion  and  matrix  reactions  to  these  loads. 


Inclusion  Reaction  to  P2  Loads 


P2  loads  are  applied  to  the  inclusion,  which  are  produced  by  cooling  the  structure 
after  adhesive  cure  is  complete.  These  loads  produce  corresponding  displacements 
within  the  inclusion,  which  can  be  used  to  determine  how  the  P2  loads  are 
distributed  between  region  D  and  the  patch.  This  procedure  is  as  follows: 
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(Al-47) 
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1  ~^xyRVyxR 
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—  +  VxyR  — 

\b  a  J 


u1  and  v1  Displacements  of  the  inclusion  due  to  P2  loads  in  the 
V  and  'y'  directions  respectively 


Where, 

Reference  [1]  useD  the  following  stress-strain  relations  for  equations 


PlxRI  ~ 


1  ^ xyR^  yxR 


P  = 

r2yRI 


^  ^xyRYyxR 


<  u1  „  S' 

v  a  b  J 

f  u1  v1 ' 

VxyRpyR  ~  +  PyR  T 

v  a  b  J 


As  can  be  seen,  the  stress-strain  relations  above  are  not  correct. 

Therefore  the  total  loads  applied  to  the  inclusion  are  as  follows,  see  Figure  22: 

£P  2x1 =P  2xPI"l"P  2xRI 
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£P 2yI=P  2yPI"*"P  2yRI 


Therefore: 
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(Al-48b) 


Although  the  stress-strain  relationship  used  in  [1]  is  incorrect  it  still  produced  the 
same  two  expressions  above.  Therefore  errors  listed  in  [1]  for  the  above  section  may 
only  be  due  to  typo  errors. 

Solving  for  displacements  yields: 


c3I4c/-c,I^/ 

C2C3  -  C,C4 


a  C3 


where 


c  _  Ex{  ,  ExR(R  Cl  _  vxy>Ext  |  v yxR(RExR 
1  l-v^;  l-vjy  1  “  vxyRvyxR 

c  _  VxyExt  i  VxyREyRtR  g  _  Ext  + 

3  l-v™  4  l-v™ 


(Al-49) 


(Al-50) 


(Al-51) 


•  Matrix  Reaction  to  P2  Loads 

The  matrix  is  taken  as  a  circular  plate  with  a  hole  in  the  centre  loaded  by  P2  loads  as 
shown  in  Figure  21.  The  equilibrium  equation  for  a  plate  with  a  loaded  circular  hole 
is  as  follows: 


d°r  +  9  =  0  (Al-52) 

dr  r 


which  has  the  following  solution  for  the  displacement 
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a2)  Cl 


u  =  ClH  1 — r  + —  (r  >  a) 


(Al-53) 


where,  Cl  and  C2  which  are  determined  by  the  boundary  conditions.  For  the 


case: 


u=0  and  ar=P/t  at  r=a 

Where,  P  Resultant  loading  around  the  hole 
R  Outer  edge  radius  of  plate 

we  obtain 


(Al-54) 
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Et 
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C2  =  — - - — x  - y— - 

Ext  (  {R2  , 

2+^1_%v _1 

Adopting  the  following  definitions: 


(Al-54) 


x-direction 

y-direction 

x-direction 

y-direction 


P=P  2xM 
P=P  2yM 
R=Rx 
R=Ry 


we  obtain  the  following  expressions: 

-P2xMa(l-Vly) 


(Al-55) 


y  i 

-  PlvMb{\.-v2\  1  b2  , 

— 2J^L1 - - i - -l - -  (Al-56) 

E  t  ,  {Rl  )  K 

2  +  (l-v,  hf-l 


This  furnishes  the  solution  for  a  plate  with  a  hole  loaded  by  P2  loads. 
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Compatibility 


Compatibility  requires  that  at  the  hole  boundary,  equal  displacements  in  the  x  and  y 
directions  must  exist  between  the  inclusion  and  surrounding  matrix.  This  leads  to  the 
following  equation  for  the  x-direction: 


=  8p  +5 


M 

P* 


Here, 


Spt  Displacement  of  inclusion  due  to  Pi  loads 

Sp2  Displacement  of  inclusion  due  to  P2  load 

Spt  Displacement  of  Matrix  due  to  Pi  loads 

Sp2  Displacement  of  Matrix  due  to  P2  loads 


(Al-57) 


Equation  Al-57  yields  the  following: 
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Reference  [1]  contains  an  error  in  the  equivalent  equation.  From  Figure  21 
equilibrium  requires; 


2P  2xI='P  2xM— P 2x 


as  a  result  we  obtain: 

P2x  =  A 
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Similarly,  matching  compatibility  in  the  y-direction  yields: 

8i+8k“8J+8if  (A1‘61) 
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from  which  we  obtain: 
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Reference  [1]  contains  an  error  in  the  equivalent  equation, 


(Al-62) 


From  Figure  21  equilibrium  requires 

2P 2yI=_P 2yM=P 2y 

from  which  we  obtain: 


(Al-63) 
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from  which  we  obtain 


P2x  =  A 
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This  set  of  equations  can  be  solved  to  give: 


=  4±zAA  and  p  = 
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where 
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Reference  [1]  does  not  derive  the  stresses  produced  by  equations  Al-69  and  Al-70.  It 
only  considers  the  stress  at  the  edge  of  the  repairand  the  actual  stresses  in  the  patch, 
i.e.  not  the  maximum  possible  stress  in  the  repair  due  to  the  presence  of  a  crack. 

CalcuRep  uses  equation  (Al-69a)  to  compute  the  stress  under  the  repair,  but  does 
not  use  equations  (Al-69b)  and  (Al-70a)  to  calculate  the  stress  at  the  edge  of  the 
patch  and  the  maximum  stress  in  the  repair,  respectively.  (See  equations  55  to  57). 
•Instead  it  calculates  the  stress  at  the  edge  of  the  patch,  not  allowing  for  the  extra 
stress  (ie  Pix)  due  to  enforcing  the  rigid  bond  assumption. 

Notation  for  section  3.2.1.2 

ap(eff)x  &  ap(eff)y  CTE  of  the  parent  structure  in  the  'x'  and  'y'  directions 
a  Parent  structure  CTE 

dx  Heated  region  length  in  the  x-direction 

dy  Heated  region  length  in  the  y-direction 

Exr  Elastic  Modulus  of  repair  in  the  x-direction 
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Ex  Elastic  Modulus  of  plate  in  x-direction  (same  as  'y'  since 

Isotropic) 

EyR  Elastic  Modulus  of  repair  in  the  y-direction 

Pixp  &  Piyp  Load  applied  to  plate  to  satisfy  the  rigid  bond 

Assumption  in  the  V  &  'y'  directions  respectively 
PixR  &  PiyR  Load  applied  to  repair  to  satisfy  the  rigid  bond 

Assumption  in  the  'x'  &  'y'  directions  respectively 
Pix  &  Piy  Load  tractions  equal  in  magnitude  to  the  above 

loads  (V  &  V  directions  respectively) 

P2x  &  P2y  Loads  applied  to  the  inclusion  in  the  'x'  &  'y'  directions 

SP2xi  &  SP2yi  Sum  of  loads  applied  to  the  plate  and  repair  ('x'  &  'y'  directions) 

P2xM  &  P2yM  Loads  applied  to  the  outer  structural  matrix 

P2xpi  &  P2yPi  Distributed  inclusion  loads  in  the  plate  (V  &  'y'  direction) 

?2xRi  &  P2yRi  Distributed  inclusion  loads  in  the  repair  ('x'  &  'y'  direction) 

Rx  Half  frame  spacing 

Ry  Stringer  spacing 

Tx  &  Ty  Temperature  at  the  edge  of  the  parent  structure  at  Rx  &  Ry 

Tc  Cure  temperature 

Tamb  Ambient  temperature 

T  Thickness  of  parent  structure 

AT  Tcuie-RT 


Tamb 

Vxy 

VxyR 

VyxR 

Vp 

ui 

V1 


Temperature  at  sea  level  during  operation 
Parent  structure  Poisson's  ratio 
Poisson's  ratio  of  repair  in  the  xy-direction 
Poisson's  ratio  of  repair  in  the  yx-direction 
Poisson's  ratio  of  parent  structure 
Displacement  of  inclusion  in  the  x-direction 
Displacement  of  inclusion  in  the  y-direction 
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1.  BACKGROUND 


Bonded  repair  technology  provides  enhanced  structural  efficiency  to  extend 
the  fatigue  life  of  aircraft  structures.  It  provides  industry  with  an  efficient  and 
cost  effective  solution  to  aircraft  structural  defects.  This  has  prompted  many 
organisations  to  use  bonded  repair  technology  as  a  preferred  repair  method 
over  mechanically  fastened  repairs. 

The  USAF  has  utilised  bonded  repair  technology  to  prolong  the  structural 
lives  of  aircraft  such  as  the  F-16  and  the  C130.  But  the  verification  process  of 
such  repairs  can  be  expensive  and  time  consuming,  since  detailed  Finite 
Element  Analysis  and  testing  may  be  required  to  ensure  an  efficient  repair. 
The  use  of  validated  analytical  procedures  for  deriving  repairs  is  therefore 
desirable  to  reduce  verification  and  testing  time.  The  US  Airforce  academy 
has  developed  a  software  package  called  "CalcuRep"  for  this  purpose. 

Robert  Fredell  and  Walter  VanBarneveld  first  developed  CalcuRep  at  Delft 
University  in  the  early  1990's.  Based  on  work  done  by  Robert  Fredell  [1]  at  the 
Delft  University  of  Technology,  CalcuRep  was  designed  to  repair  cracks 
specifically  on  an  Aluminium  aircraft  fuselage.  This  initial  version  of  the 
software  was  developed  in  DOS  and  provided  the  user  with  repair  design 
tools  with  a  basic  user  interface.  This  prompted  the  development  of  CalcuRep 
for  Windows  (Beta)  at  the  US  Air  Force  Academy  in  1994,  which  consisted  of 
a  user-friendly  menu  system  and  a  graphical  interface.  CalcuRep  for 
Windows  Version  1.0  [2],  released  in  1997,  incorporated  user  requirements 
and  suggestions.  • 

The  RAAF  Engineering  Standard  C5033  [3,  4,  5]  was  released  in  1994  to 
provide  bonded  repair  design  tools  to  assist  in  the  repair  of  structural  defects 
in  aircraft  used  by  the  RAAF.  Aimed  more  towards  providing  engineers  with 
a  set  of  design  standards  rather  than  a  software  package,  C5033  provides  a  set 
of  constraints  to  ensure  quality  of  repair  design,  specifications  and  application 
processes. 

In  early  1999,  the  US  Air  Force  Academy  awarded  a  contract  through  the 
DSTO  Centre  of  Expertise  at  Monash  University  to  verify  and  extend  the 
capability  of  CalcuRep  version  1.0  for  Windows.  The  work  is  noted  as  the 
CalcuRep  review  task.  The  verification  aspect  of  the  task  (know  as  part  1  of 
the  work)  was  performed  by  comparing  the  CalcuRep  approach  with  the 
RAAF  Engineering  Standard  C5033  procedures.  A  significant  example  case  of 
a  RAAF  F-111C  wing  repair  was  used  to  compare  the  CalcuRep  and  C5033 
repair  design  techniques.  The  results  are  detailed  in  reference  [8], 

Part  2  of  the  task  was  to  extend  the  capability  of  CalcuRep  to  include 
unsupported  single-sided  repairs  and  repairs  to  stiffened  panels.  The  results 
for  part  2  are  detailed  in  this  report. 


The  work  compiled  in  this  report  establishes  a  new  design  approach  from  the 
most  recent  works  developed  at  the  DSTO  in  the  field  of  bonded  repair 
design.  Applying  these  recent  findings  in  such  a  manner  to  design  repairs  will 
prove  to  benefit  both  the  DSTO  and  the  USAF.  The  content  of  this  report  may 
also  be  considered  by  the  RAAF  in  expanding  the  current  RAAF  Engineering 
Standard  C5033  to  deal  with  single-sided  repairs  to  unsupported  structure. 
But  overall,  this  task  proves  to  be  beneficial  to  both  the  USAF  and  the  DSTO 
by  attempting  to  reduce  any  differences  between  the  US  and  Australian  repair 
design  methodologies. 

2.  INTRODUCTION 

An  evaluation  and  comparison  between  the  bonded  repair  design  software 
package  CalcuRep  version  1.0  and  the  RAAF  Engineering  Standard  C5033  for 
single  sided  fully  supported  repairs  confirmed  a  good  correlation  of  results. 
The  focus  of  this  report  is  to  expand  and  develop  CalcuRep  to  deal  with  a 
wider  range  of  repair  cases.  As  a  result,  this  has  led  to  a  revised  package, 
which  considers  the  best  aspects  of  CalcuRep  and  C5033  and  extends  its 
capability.  This  package  has  been  designated  as  the  "USA/ Australia  Bonded 
Repair  Design  Module".  For  convenience,  this  module  has  been  developed  in 
a  software  package  called  Mathcad  [9]  to  provide  users  with  visual  inspection 
of  equations  and  quick  accessibility  of  user  input.  Although  the  package  is  in 
Mathcad  form,  it  can  be  programmed  by  any  software  package  the  user 
desires,  while  the  Mathcad  sheet  provided  in  this  report  can  provide  the 
means  to  verify  the  package. 

The  package  allows  the  user  to  design  bonded  repairs  to  metallic  aircraft 
structure  under  the  following  circumstances: 

•  Double  sided  repairs  and/  or  single  sided  fully  supported  repairs  to  centre 
cracked  biaxially  loaded  metallic  structure.  Thermal  and  applied  in-plane 
loading  effects  are  included.  The  scope  of  this  work  allows  for  plastic 
behaviour  in  the  adhesive  should  the  adhesive  become  plastic. 

•  Single  sided  repairs  to  unsupported  centre  cracked  biaxially  loaded 
metallic  structure.  Thermal  and  applied  in  plane  and  out  of  plane  loading 
effects  are  included.  The  scope  of  this  work  is  limited  to  a  geometric  linear 
analysis  and  elastic  behaviour  in  the  adhesive.  Justification  of  these  points 
will  be  described  in  this  report. 

•  Repair  to  structures  with  riveted  stringers.  The  scope  of  this  work  can  be 
applied  to  two  scenarios,  a  stiffened  panel  with  a  crack  in  the  stiffener  and 
a  stiffened  panel  with  a  centre  crack  in  the  skin  beneath  the  centre  stringer. 


3.  APPROPRIATE  TOOLS  FROM  THE  CALCUREP  AND  C5033 
APPROACHES 

3.1  New  Approach 

Part  1  of  the  CalcuRep  task  [8]  outlined  the  approach  adopted  by  the  RAAF 
Engineering  standard  C5033  and  CalcuRep  Version  1.0  for  adhesive  bonded 
repair  design.  Many  limitations  were  identified  in  the  single  sided  repair 
methodologies  used  by  the  two  approaches  [8]  and  are  summarised  below: 

C5033 

•  Unable  to  deal  with  biaxial  loading  systems 

•  Limited  to  single  sided  repairs  to  fully  supported  structure 

CalcuRep 

•  Restraint  against  user  input 

•  Does  not  allow  the  user  to  change  the  adhesive  properties  with 
temperature  via  the  adhesive  data  library 

•  Does  not  allow  operational  temperatures  above  see  level  to  be  used. 

•  Does  not  provide  a  method  to  estimate  preliminary  patch  dimensions 

•  Limited  to  single  sided  repairs  to  fully  supported  structure  ' 

•  Contains  errors  in  the  out  of  plane  bending  and  thermal  residual  stress 
analysis  tools. 

Many  repair  cases  consist  of  cracks  in  unsupported  structure  accessible  from 
only  one  side  of  a  panel.  In  such  cases,  out  of  plane  bending  effects  become 
significant  and  cannot  be  ignored.  Feasible  preliminary  design  is  therefore 
difficult  using  packages  such  as  C5033  and  CalcuRep  due  to  the  limitations 
listed  above. 

This  type  of  analysis  requires  the  use  of  more  advanced  analysis  tools  to 
account  for  out  of  plane  bending  effects  inherent  in  single  sided  unsupported 
repairs.  Tools  that  deal  with  these  effects  have  been  developed  [10, 11, 12, 13] 
and  will  be  discussed  in  detail  latter  in  this  report.  The  approach  developed  in 
this  report  attempts  to  utilise  these  tools  to  provide  users  with  the  ability  to 
design  repairs  for  a  wider  range  of  cases.  This  package  has  been  named  the 
"USA/ Australia  Bonded  Repair  Design  Module".  It  can  deal  with  the 
following  repair  scenarios: 

•  Single  Sided  Fully  supported  repairs 

•  Double  Sided  repairs 

•  Single  Sided  Unsupported  repairs 

•  ..  Stiffened  Panels  with  a  crack  in  either  the  stiffener  or  the  skin  beneath  the 

stiffener 


CalcuRep  and  C5033  contain  good  design  tools  for  single  sided  supported 
repairs  and  double  sided  repair.  Since  the  new  design  module  will  account  for 
such  repair  cases,  it  is  beneficial  to  use  some  of  these  existing  tools.  The 
following  sections  will  provide  a  brief  summary  of  useful  tools  from  C5033 
and  CalcuRep  for  the  new  approach  as  listed  in  reference  [8]. 

3.2  Overview  of  Appropriate  C5033  and  CalcuRep  tools  for  the  New 
Approach 

Part  1  of  the  CalcuRep  review  task  [8]  outlined  an  overview  of  the  RAAF 
Engineering  Standard  C5033  approach.  To  summarise  C5033's  capability,  the 
standard  as  able  to  deal  with  the  following  repair  cases: 

•  Double  Sided  Bonded  Repairs 

•  Single  Sided  Supported  Bonded  Repairs 

•  Single  Sided  Supported  Bonded  Repairs  to  a  Part  Through  Thickness 
Defect 

•  Single  Sided  Supported  Bonded  Doubler 

•  Double  Sided  Bonded  Doubler 

•  Single  Sided  Supported  Bonded  Repairs  for  Stress  Corrosion  Cracking 

•  Double  Sided  Bonded  Repairs  for  Stress  Corrosion  Cracking 

•  Scarf  Repairs  to  Composite  Structures  ' 

Part  1  of  the  CalcuRep  review  task  [8]  also  outlines  an  overview  of  the 
CalcuRep  version  1.0  approach.  To  summarise  briefly,  CalcuRep  can  provide 
users  with  analysis  tools  to  deal  with  the  following: 

•  Single  sided  fully  supported  repairs  to  aluminium  structure. 

•  Biaxial  stress  fields. 

•  Orthotropic  patch  bonded  to  an  isotropic  structure. 

•  Provides  the  out  of  plane  bending  stress  component  at  the  edge  of  the 
patch  in  the  parent  structure. 

CalcuRep,  similarly  to  C5033,  is  limited  to  single  sided  fully  supported 
repairs,  but  proves  to  be  advantageous  in  dealing  with  orthotropic  patches 
and  biaxial  stress  fields. 

A  summary  of  appropriate  tools  from  C5033  and  CalcuRep  is  as  follows: 

3.2.1  Patch  Dimensions 

C5033  provides  the  means  for  determining  preliminary  patch  dimensions  and 
can  be  adapted  to  the  following  cases. 


Single-sided  fully  supported  repairs  for  uniaxial  loadings  (See  sections  3.1.1.2 
and  3.1.1.5  of  [8]  for  details). 

Double-sided  repairs  for  uniaxial  loadings. 

For  a  two-sided  repair  case,  the  stiffness  of  the  patch  is  doubled  [3,  4,  5]. 
Therefore  the  required  patch  thickness  can  be  determined  by  modifying 
equation  1  of  reference  [8]  as  follows: 

EjH_  <  t  <  (1) 

2 E0-‘°~  2  E0 

where,  Ei  Elastic  Modulus  of  the  Parent  Material 
E0  Elastic  Modulus  of  the  Patch 
h  Thickness  of  the  Parent  Material 

t0  Thickness  of  the  Patch 

The  patch  length  for  two  sided  repairs  [3,  4,  5]  can  also  be  determined  by 
modifying  equations  11  and  25  to  27  of  reference  [8]  as  follows: 


where  Lt  Total  transfer  length 

L  Overlap  Length 

le  Elastic  Load  Transfer  Length 

lp  Plastic  Load  Transfer  Length 

X  Elastic  Exponent 

G  Adhesive  Shear  Modulus 

r|  Adhesive  thickness 

ctulti  Ultimate  Strength  of  the  Parent  Material 

The  patch  width  for  a  double  sided  repair  is  determined  using  the  same 
method  as  for  a  single  sided  supported  repair  and  can  be  calculated  using 
section  3.1. 1.5  of  [8],  To  summarise,  the  patch  is  required  to  extend  at  least  1 
inch  over  both  crack  tips. 

For  a  biaxial  load  case,  load  transfer  is  assumed  in  both  directions.  The  patch 
length  calculation  remains  unchanged  and  can  still  be  found  as  shown 


previously.  But  in  such  a  case,  the  patch  width  is  no  longer  determined  from 
section  3.1.1 .5  of  [8],  but  is  calculated  as  follows: 

Patch  width  for  single  sided  supported  repair  (Biaxial  Loading) 


LT  =  K  +  lp  (6) 

'•  -  i <7) 


Patch  width  for  double  sided  repair  (Biaxial  Loading) 

Lt  =  l,  +  l„  (10) 

l.-T  (11) 


where,  Eyr  is  the  elastic  modulus  of  the  patch  in  the  y-direction 

Note:  The  total  length  of  pre-existing  bond-line  defects  must  be  taken  into 
account  in  determining  patch  dimensions  [3,8]. 

3.2.2  Evaluate  the  Rapid  Repairability  Criterion  (RRC)  C5033 

The  rapid  repairability  criterion  is  a  management  tool  to  ensure  the  adhesive 
is  never  the  critical  element  in  a  repaired  joint,  see  section  3.1.1.2  of  reference 
[8]  for  details.  Section  3.1. 1.2  [8]  outlines  the  rapid  repairability  criterion  for 
fully  supported  single  sided  repairs  and  can  be  expanded  to  double-sided 
repairs  by  replacing  equations  3  and  4  in  [8]  by  the  following: 

(14) 
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where, 

P  Load  capacity  of  the  adhesive 

<Xo  Coefficient  of  Thermal  Expansion  of  the  Patch 

ai  Coefficient  of  Thermal  Expansion  of  the  Parent  material 

aieff  Effective  Coefficient  of  Thermal  Expansion  of  the  Parent 

material 

RT  Room  Temperature 

Tcure  Cure  Temperature 
Toper  Operating  Temperature 
tp  Adhesive  Plastic  Shear  Stress  Limit 

Ye  Adhesive  Elastic  Shear  Strain  Limit 

yp  Adhesive  Plastic  Shear  Strain  Limit 


3.2.3  Check  the  Adhesive  Shear  Strain 


The  adhesive  shear  strain  is  calculated  using  equations  10  to  12  from  section 
3. 1.1.3  in  [8]  for  single  sided  fully  supported  repairs.  Both  C5033  and 
CalcuRep  use  these  equations.  To  expand  these  equations  to  deal  with  two 
sided  repairs,  equations  10  to  12  in  [8]  can  be  changed  as  follows^ 


•  Elastic  behaviour  in  the  adhesive 
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Plastic  behaviour  in  the  adhesive 
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where,  go  is  the  stress  in  the  parent  structure  under  the  repair  and  where  X  is 
given  by  equation  5. 

3.2.4  Peel  Stress 

Peel  stress  is  predominantly  due  to  load  path  eccentricity  inherent  in  single 
sided  repairs.  Although  load  path  eccentricity  effects  are  eliminated  when  a 
double-sided  repair  or  a  fully  supported  single  sided  repair  is  used,  a 
localised  peel  stress  is  still  evident  when  each  bond-line  is  considered  in 
isolation.  Section  3.1.1.3  of  [8]  outlines  the  equation  used  to  calculate  the  peel 
stress  for  both  fully  supported  single-sided  and  double-sided  repairs.  The 
equation  is  as  follows: 


(15) 


(18) 


where  ac  max  Maximum  Peel  Stress 

Ec  Effective  Transverse  Tensile  Modulus  for  the  Adhesive 

3.2.5  The  Upper  Bound  Stress  Intensity  Factor 

The  upper  bound  stress  intensity  factor  for  single  sided  supported  repairs  is 
detailed  in  section  3.1.1.3  step  5  of  reference  [8].  Both  C5033  and  CalcuRep 
use  this  analysis  tool.  Equation  13  and  14  of  [8]  can  be  changed  to  deal  with 
double-sided  repairs  as  follows: 

•  Elastic  Behaviour  in  the  adhesive 


•  Plastic  behaviour  in  the  adhesive 


4.  NEW  TECHNIQUES 

4.1  Single  Sided  Repair  Methods 

Double-sided  repair  methodology  was  initially  idealised  by  a  double-lap  joint 
system.  This  system  was  found  to  be  ideal  since  no  load  path  eccentricity 
exists.  Shortly  following,  it  was  found  that  most  repair  cases  required  patch 
application  to  one  surface  rather  than  two,  due  to  a  lack  of  accessibility  of 
both  surfaces  of  many  aircraft  structure,  ie  cracks  in  wings.  This  prompted  the 
development  of  single-sided  repair  analysis  tools.  Previous  analytical  theory 
for  two-sided  repairs  were  then  manipulated  by  cutting  the  repair  geometry 
in  half.  This  effectively  produced  the  analysis  tools  for  fully  supported  single¬ 
sided  repairs.  It  was  then  found  that  these  tools  were  not  sufficient  to  design 
for  single-sided  repair  cases  where  insufficient  localised  structural  support 
was  provided  and  load  path  eccentricity  effects  became  apparent.  As  a  result, 
repairs  were  found  to  exhibit  the  out  of  plane  bending  phenomena.  Works 
from  references  [10,  11,  14,  15]  have  developed  analytical  tools  to  deal  with 
load  path  eccentricity  effects. 

Existing  tools  to  date  can  deal  with  the  following  cases: 


•  Fully  supported  single  sided  repairs  with  applied  stresses  [16] 

•  Out  of  Plane  Bending  in  unsupported  single  sided  repairs  due  to  applied 
stresses  [10, 11] 

•  Fully  supported  single  sided  repairs  with  thermal  residual  stresses  [17] 

Recent  works  at  the  DSTO  by  C.H.  Wang  and  D.  Erjavec,  have  developed 
tools  that  consider  the  following  systems: 

•  In  plane  and  out  of  plane  stresses  in  unsupported  single  sided  repairs  due 
to  applied  stresses  [13] 

•  In  plane  and  out  of  plane  stresses  in  unsupported  single  sided  repairs  due 
to  thermal  residual  stresses  [12] 

The  following  sections  will  outline  some  of  these  most  recent  analytical 
contributions  to  bonded  repair  analysis  methodology. 

4.2  In  plane  and  Out  of  plane  loading  due  to  Thermal  stresses 

A  commonly  used  thermal  residual  stress  model  (equation  6  [8])  developed  in 
reference  [18],  assumes  a  one-dimensional  fully  supported  repair  system 
consisting  of  an  isotopic  patch  and  plate.  Such  an  approach  is  very  simplistic 
and  does  not  consider  the  effects  of  an  orthotropic  patch  in  a  two-dimensional 
stress  field,  or  the  load  attracted  into  the  repair  region  due  to  the  extra 
stiffness  provided  by  the  patch.  A  fully  supported  single  sided  repair  model 
developed  by  Wang,  C.  H.,  Rose,  L.  R.  F.,  Callinan,  R.  and  Baker,  A.  A  [17] 
deals  with  these  issues  by  the  use  of  Eshelby's  inclusion  analogy  [20].  The 
report  derives  equations  for  thermal  residual  stresses  post  cure  considering 
the  following  two  cases;  1)  an  orthotropic  patch  bonded  on  to  an  infinite 
isotropic  plate;  2)  an  orthotropic  patch  bonded  on  to  a  finite  isotropic  plate. 

Since  reference  [17]  is  limited  to  fully  supported  single-sided  repairs, 
expanding  the  analysis  to  unsupported  single  sided  repairs  would  enhance 
the  models  capability.  This  has  lead  to  the  development  of  design  tools  for 
thermal  residual  stresses  post  repair  for  unsupported  single  sided  repairs  by 
C.  FI.  Wang  and  D.  Erjavec  [12].  The  following  provides  the  results  of  the 
formulation  as  shown  in  [12]: 

Repaired  Structure  Geometry 

The  curvature  and  mean  strains  are  postulated  to  be  uniform  throughout  the 
inclusion  [12].  The  inclusion  region  refers  to  the  repair  and  the  plate  under  the 
repaired  region  as  one  entity.  A  free  body  diagram  of  the  system  is  as  follows: 


Where,  Ri 

Ro 


Figure  1 

Radius  of  patch 
Radius  of  plate 


M 


Equilibrium  requires  the  following: 


Mx=mx  +  Fxz*  (21) 


My  ~  My  +  Fy  Zy  (22) 


where, 

Mx  My  are  the  total  Moments  in  the  V  and  'y' 
directions  respectively 

mxmy 

Moment  in  the  plate  region  in  the  V  and  'y' 

directions  respectively 

Fx  Fy 

Total  Force  in  the  'x'  &  'y'  directions  respectively 

Zx  Zy 

Neutral  axis  offset  in  the  'x'  &  'y'  directions 
respectively 

The  stresses  in  the  plate  and  reinforcement  inside  the  repaired  region  are 
related  through  Hooke's  law  as  follows: 


(23) 


axx^\-  A<  £x~Kx(z~zx)~a&T 
ayy(z\  J  sy-Ky{z  -  zy)~  aAT 

I  sx-Kx(z-zx)-alRAT 

AR  _ 

£y-Ky(z-zy)-a2RAT 

Where,  axx(z)  &  a>T(z)  Stress  variation  through  the  plate  (V  &  'y' 

directions  respectively) 

CxxR(z)  &  OyyR(z)  Stress  variation  through  the  repair  ('x'  &  'y' 
directions  respectively) 

£x  &  sy  Mean  strains  inside  the  inclusion  ('x'  &  ‘y’ 

directions  respectively) 

Kx  &  Ky  Curvatures  inside  the  inclusion  ('x'  &  'y' 

directions  respectively) 

a  Coefficient  of  thermal  expansion  of  the 

plate 

aiR  Coefficient  of  thermal  expansion  of  the  repair  in  the  x- 
direction 

<X2R  Coefficient  of  thermal  expansion  of  the  repair  in  the  y- 
direction 

AT  Temperature  change 

A  In-plane  Plate  Properties 

Ar  In-plane  Repair  Properties 

The  total  forces  and  moments  can  be  determined  as  follows: 


where. 


In-plane  Inclusion  Properties 

Ai=-^—A  +  -^—Ar 

t  +  tR  t  +  tR 


(27) 
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Rigidity  of  Inclusion 
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Where,  Ex 

Parent  structure  Elastic  Modulus 

Exr 

Repair  Elastic  Modulus  in  the  x-direction 

EyR 

Repair  Elastic  Modulus  in  the  y-direction 

t 

Parent  structure  thickness 

tR 

Repair  thickness 

V12 

Poisson's  ratio  of  Parent  structure 

V12R 

Poisson's  ratio  of  repair  in  the  xy-direction 

(35) 


V21R  Poisson's  ratio  of  repair  in  the  yx-direction 

By  determining  the  displacement  and  deflection  functions  for  the  inclusion 
and  outer  structure,  the  mean  strains  and  curvatures  inside  the  inclusion  can 
be  determined  by  matching  compatibility  conditions  at  the  inclusion 
boundary.  Enforcing  these  compatibility  conditions  results  in  the  following 
formulation. 

Sx 

£y 


(36) 
.  2  J 


(45) 


Where,  Ti  Cure  temperature 

To  Temperature  in  parent  structure  at  a  distance  R0 

Having  determined  the  mean  strains  and  curvatures  of  the  inclusion,  the 
corresponding  stress  in  the  plate  and  repair  is  given  by  equations  23  and  24. 
The  results  above  reduce  to  the  results  presented  in  reference  [17]  by  setting 
the  curvatures  in  equation  36  to  zero,  which  yields  the  following: 


The  above  analysis  was  compared  to  a  fully  3-D  finite  element  repair  case 
using  Abaqus  6.8  for  a  number  of  repair  material  variations.  The  results 
correlated  well  against  the  FE  results. 

4.3  In  Plane  and  Out  of  Plane  Loading  due  to  Applied  Stresses 

As  discussed  previously,  out  of  plane  bending  is  due  to  load  path  eccentricity 
inherent  in  single  sided  repairs.  Its  effects  cannot  be  neglected  when  a  single¬ 
sided  repaired  structure  is  unsupported.  Previous  works  conducted  on  this 
subject  can  be  found  in  references  [14,  15].  One  of  the  most  commonly  used 
out  of  plane  bending  formulations  was  derived  by  L.R.F  Rose  [14]  which 
considers  a  one-dimensional  strip  model  of  a  repaired  structure  under  an 
uniaxial  applied  stress  and  is  summarised  in  section  3.2.1.3  of  reference  [8].  It 
derives  the  stress  variation  through  the  thickness  of  the  parent  structure  and 
repair  due  to  a  neutral  axis  offset  using  a  non-linear  analytical  formulation. 

The  most  recent  works  conducted  in  this  area  are  listed  in  references  [10, 11] 
and  also  consider  a  one-dimensional  model  with  applied  stresses,  but 
introduces  additional  applied  bending  moments  to  account  for  acoustic 
fatigue  problems.  Both  reports  comprise  of  a  linear  stress  analysis  [10],  and  an 
approximate  non-linear  stress  analysis  [11],  excluding  load  attraction  effects. 

Expanding  these  methodologies  to  include  biaxial  loading  systems  and  load 
attraction  effects  would  be  beneficial  in  providing  a  set  of  comprehensive 
analysis  tools  for  this  particular  repair  case.  This  has  lead  to  the  development 
of  reference  [13]  by  C.  FI.  Wang  and  D.  Erjavec,  which  deals  with  these  effects. 
The  following  provides  a  brief  summary  of  the  findings: 

The  loading  configuration  in  the  x-direction  is  as  follows: 


Similarly  in  the  y-direction 


Where,  axR  Applied  stress  in  the  x-direction  at  the  edge  of  the 
structure 


structure 


ayR  Applied  stress  in  the  y-direction  at  the  edge  of  the 
MXR  Applied  moment  per  unit  length  in  the  x-direction  at  the 

edge  of 

the  structure 

MyR  Applied  moment  per  unit  length  in  the  y-direction  at  the 

edge  of 

the  structure 

Free  Body  Diagram  of  Repaired  Region 
M 


Figure  5 


Equilibrium  requires  the  following: 


Mx=mx+Fxzx  (52) 

Similarly  in  the  y-direction 

My=my+FyZy  (53) 


The  stresses  in  the  plate  and  reinforcement  inside  the  repaired  region  are 
related  through  Hooke's  law  as  follows: 
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Sy-K  (Z-Zy) 


(54) 


S.-K^Z-Z.)  (M) 
Sy-Ky(z~Zy ) 


(57) 


1 _ 

“  l4/]j  *}  +  (t  + 

4 

; 

*  2  ^  (?  +  )[^ 

f*. 

L 

LAJ  L  . 

Vs  >\ 

By  determining  the  displacement  and  deflection  functions  for  the  inclusion 
and  outer  structure,  the  mean  strains  and  curvatures  inside  the  inclusion  can 
be  determined  by  matching  the  compatibility  conditions  at  the  inclusion' 
boundary.  Enforcing  these  compatibility  conditions  results  in  the  following 
formulation. 
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The  remaining  constants  are  given  in  section  4.2  by  setting  the  temperature 
terms  to  zero. 

Verification  of  this  formulation  will  be  provided  on  release  of  reference  [13]. 
In  the  mean  time,  a  repair  case  study  has  been  provided  in  section  6  of  this 
report  and  provides  a  preliminary  means  of  verifying  section  4.3  (To  be 
discussed  later). 


4.4  Stress  intensity 


One  of  the  most  commonly  used  formulations  of  a  stress  intensity  factor  was 
derived  by  L.R.F  Rose  [14,  16]  and  is  given  by  equations  19  and  20.  These 
expressions  can  only  be  applied  to  fully  supported  single-sided  repairs  and 
double-sided  repairs,  since  they  do  not  take  load  path  eccentricity  into 
account.  Two  such  reports  that  deal  with  stress  intensity  factors  due  to  load 
patch  eccentricity  are  references  [10]  and  [11].  Reference  [10]  determines  the 
limiting  strain  energy  release  rate  for  an  unsupported  single-sided  one¬ 
dimensional  repair  due  to  applied  loads.  In  such  a  case,  load  path  eccentricity 
causes  a  non-uniform  stress  under  the  repair,  thus  resulting  in  a  varying 
stress  intensity  factor  through  the  thickness.  The  report  found  that  when  load 
path  eccentricity  effects  were  included,  the  stress  intensity  factor  still 
exhibited  asymptotic  properties  as  crack  length  increased,  as  shown  in  figure 
12  of  reference  [8].  This  report  provides  expressions  for  maximum,  minimum 
and  root  mean  square  limiting  stress  intensity  factors  in  the  bounds  of  a 
geometrical  linear  analysis. 

Reference  [11]  provides  a  crack-bridging  model  that  replaces  the  patch  by  a 
set  of  distributed  springs  bridging  the  two  cracked  faces.  Similarly, 
considering  a  one-dimensional  unsupported  single-sided  repair,  the  report 
considers  a  geometrically  linear  analysis  and  an  approximate  geometrical 
non-linear  analysis.  The  findings  produced  an  equivalent  spring  compliance 
matrix  representing  the  patch  and  determined  the  corresponding  variations  of 
stress  intensity  factor  through  the  thickness  using  Kirchhoff's  Poisson  plate 
bending  theory  and  Reissner's  shear  deformation  theory. 

Repairs  that  are  unsupported  can  yield  peel  delamination  shortly  after  the 
adhesive  become  plastic  (See  the  strap  joint  in  the  figure  below).  Therefore  the 
scope  of  references  [10]  and  [11]  retain  the  adhesive  in  the  elastic  regime. 


Formulation 

Using  reference  [10]  and  [11],  a  solution  can  be  determined  for  the  limiting 
minimum,  maximum  and  root  mean  squared  stress  intensity  factors  as 
functions  of  the  mean  and  bending  stresses.  The  following  describes  the 
analysis: 


The  total  mean  force  and  bending  moment,  assuming  no  crack  to  be  present, 
can  be  given  as  follows  [10]: 


N  =  crt 

(65) 

u  =  ct‘,! 

(66) 
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(68) 

As  the  crack  propagates,  work  is  extracted  to  relax  the  first  stage  stresses  to 
zero,  (ie,  first  stage  stresses  refer  to  the  stress  calculated  in  the  parent 
structure  under  the  repair  from  sections  4.2  &  4.3)  The  strain  energy  release 
rate  can  be  found  by  the  following  equation  [10]: 

(69) 

where.  Goo*  Asymtopic  strain  energy  release  rate 
u0  Crack  opening  displacement 

Go  Crack  face  rotation 

Using  the  results  from  [11],  the  crack  opening  displacement  and  crack  face 
rotation  can  be  expressed  as  functions  of  the  mean  force  and  bending  moment 
as  follows: 


OBP 

Figure  7 

Mean  stress  component 
Bending  stress  component 

Stress  at  the  top  of  the  plate  in  the  major  direction 
(V  direction) 

Stress  at  the  bottom  of  the  plate  in  the  major 
direction  (V  direction) 


where. 


The  corresponding  asymptotic  maximum  and  minimum  stress  intensity 
factors  (ie  due  to  a  varying  stress  through  the  thickness)  can  be  determined  as 
follows: 
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4.5  Maximum  Adhesive  Shear  Strain  Under  the  repair 


Equations  10  [8],  12  [8],  16  and  17  were  developed  assuming  a  uniform-stress 
state  through  the  thickness  of  the  parent  structure.  They  do  not  consider  out 
of  plane  bending  effects  and  therefore  cannot  be  used  to  evaluate  the 
maximum  adhesive  shear  strain  in  unsupported  single  sided  repairs. 
Reference  [11]  derives  the  maximum  adhesive  shear  strain  likely  to  occur 
under  the  repair  including  load  patch  eccentricity  effects.  The  formulation  is 
in-terms  of  the  mean  force  and  bending  moment  (equations  65  and  66),  which 
follows  on  from  the  results  obtained  from  the  previous  section,  4.4.  The 
maximum  adhesive  shear  strain  under  the  repair  is  given  by  the  following 
[11]: 
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The  shear  strain  occurring  at  the  crack  location  is  given  by  -yAMax- 

4.6  Peel  Stresses 

Hart-Smith  [21]  derived  the  maximum  peel  stress  at  the  end  of  the  repair 
(Equation  18),  as  used  by  C5033  [3].  This  expression  is  derived  considering  a 
one-dimensional  lap  repair  joint  subject  to  applied  loads.  Considering  vertical 
equilibrium  requirements  allow  the  peel  stresses  to  be  solved  as  a  function  of 
position  from  the  centre  of  the  patch. 

More  recent  works,  reference  [11],  derives  a  more  complex  expression  for  a 
strap  joint,  which  is  more  representative  of  a  bonded  cracked  plate  than  a  lap 
joint  configuration  (See  figure  6).  The  patch  is  replaced  by  distributed  springs, 
which  are  used  to  formulate  the  effective  spring  compliance  matrix  (Equation 
(70)).  The  peel  strain  expression  is  determined  from  the  following  [11]: 


sA  =  (dcos(/a)  +  5sin(/a))e  **  (84) 


where  'A'  and  'B'  can  be  determined  from  the  following: 
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(Note:  Constants  used  in  equations  above  are  previously  defined) 

Knowing  the  peel  strain  distribution  from  equation  84,  the  peel  strain  near  the 
crack  and  at  the  end  of  the  patch  can  be  found  by  substituting  x=0  and  x=a  (ie 
a=half  patch  length)  respectively.  The  corresponding  peel  stress  is  given  by 
the  following: 


=  Easa 

(86) 

Ea=2Ma(1  +  va) 

(87) 

where. 

Ea  Elastic  modulus  of  adhesive 

va  Poisson's  Ratio  of  adhesive 

5.  Recommended  Approach  for  an  Improved  Design  ' 

C5033  and  CalcuRep  provide  designers  with  analysis  tools  to  develop 
appropriate  repair  designs  for  a  range  of  cases.  These  tools  have  a  number  of 
limitations  as  previously  discussed.  The  development  of  more  advanced  tools, 
as  discussed  in  section  4,  have  removed  many  of  these  limitations.  Therefore 
it  would  be  beneficial  to  apply  these  tools  in  an  appropriate  manner  to 
provide  the  means  to  design  a  repair. 

This  section  entails  the  development  of  a  design  approach  using  new  and 
already  established  design  tools  to  provide  a  more  versatile  package.  The 
package  deals  with  the  following  cases: 

•  Single  Sided  Fully  supported  repairs 

•  Double  Sided  repairs 

•  Single  Sided  Unsupported  repairs 

•  Stiffened  Panels  with  a  crack  in  either  the  stiffener  or  the  skin  beneath  the 
stiffener 

The  following  will  outline  the  steps  and  tools  required  to  produce  such  a 
package. 


5.1  Repair  Methodology  for  Single  Sided  Unsupported  Repairs 

5.1.1  Patch  Dimensions  (Single  Sided  Unsupported  Repairs) 

Determine  the  patch  dimensions  using  section  3.2.1 

Note:  The  load  transfer  length  used  by  C5033,  i.e.  equation  (11)  [8]  and  (5),  has 
been  derived  assuming  plane  stress  conditions.  References  [10]  and  [11]  show 
that  plane  strain  conditions  prevail  rather  than  plane  stress  conditions,  see  [8] 
for  details.  Therefore  the  load  transfer  length  used  for  the  following  sections 
will  be  changed  to  equation  76. 

5.1.2  Determining  the  Stress  Under  and  at  the  edge  of  the  Repair  (Single 
Sided  Unsupported  Repairs) 

•  Applied  Stresses 

Using  section  4.3,  calculate  the  stress  at  the  top  and  bottom  of  the  parent 
structure  under  the  repair  in  the  x-direction  and  denote  their  values  as  otp 
and  obp  respectively.  Calculate  the  total  force  in  the  parent  structure  at  the 
edge  of  the  repair  in  the  x-direction,  denoted  its  value  as  Fx.  ' 

•  Thermal  Stresses  , 

The  thermal  residual  stress  variation  through  the  thickness  under  the  repair  is 
effected  by  two  thermal  cycles.  The  first  is  the  cure  cycle,  while  the  second  is 
the  operational  cycle.  Consider  the  effect  these  two  cycles  have  on  the 
structure  outside  the  repaired  region  and  inside  the  repaired  region: 

Outside  the  Repaired  Region 

Cure  Cycle 

During  the  cure  cycle  the  repaired  region  is  heated  to  the  cure  temperature 
(Ti).  At  the  edge  of  the  structure  away  from  the  patch,  the  temperature  is  a 
lower  temperature  T0.  Due  to  this  non-uniform  temperature  distribution, 
constraint  is  provided  from  the  cooler  material  that  effectively  resists  the 
expansion  of  the  warmer  material.  This  phenomenon  effectively  reduces  the 
thermal  expansion  coefficient  of  the  plate,  given  by  equation  50.  It  can  be 
assumed  that  T0  is  at  room  temperature  (RT)  provided  the  edge  of  the  parent 
structure  is  sufficiently  remote  from  the  heated  region. 

After  cure,  the  heated  region  of  the  structure  is  cooled  down  to  room 
temperature.  Therefore  a  negative  temperature  change  occurs,  which  is  the 
reversal  of  the  heating  cycle.  As  a  result,  Ti  and  T0  become  negative. 


Operational  Cycle 

During  this  cycle  it  can  be  assumed  that  a  uniform  temperature  change  exists 
throughout  the  repaired  structure.  This  can  occur  in  aircraft  fuselage,  as  the 
aircraft  gains  altitude.  Therefore  T0/Ti  can  be  set  to  unity  in  equation  50. 

The  Inclusion  Region  (Inside  the  repaired  region) 

Cure  Cycle 

During  the  cure  cycle,  the  inclusion  experiences  the  following  temperature 
change 


AT=RT-Ti 


Operational  Cycle 

During  the  operational  cycle,  the  inclusion  experiences  the  following 
temperature  change. 

AT=Ti-RT 

Where,  AT  Temperature  Change 

Toper  Operational  Temperature  } 

RT  Room  Temperature 

Using  section  4.2,  calculate  the  stress  at  the  top  and  bottom  of  the  parent 
structure  under  the  repair  in  the  x-direction  and  denote  their  values  as  otp 
and  <7BP  respectively.  Calculate  the  total  force  in  the  parent  structure  at  the 
edge  of  the  repair  in  the  x-direction,  denoted  its  value  as  Fx.  This  procedure 
should  be  completed  twice,  for  both  the  cure  and  operational  thermal  cycles. 
Then  add  the  corresponding  stresses  from  these  two  cycles  to  produce  the 
total  thermal  stress  components. 

5.1.3  Finding  the  Total  Stress  (Single  Sided  Unsupported  Repairs) 

Add  the  following  components,  cttp,  obp  and  Fx,  calculated  from  the  applied 
and  thermal  stress  procedures.  The  margin  of  safety  of  safety  as  per  [3,  4,  5] 
for  the  stress  in  the  parent  structure  at  the  edge  of  the  repair  is  as  follows: 

MOS  =  ^ 1-1 
'  FJt 


where 
MOS  >  0.2 


(88) 


5.1.4  Determining  the  Maximum  Adhesive  Shear  Strain  (Single  Sided 
Unsupported  Repairs) 

The  maximum  adhesive  shear  strain  can  be  determined  using  section  4.5.  The 
margin  of  safety  as  per  [3,  4,  5]  for  the  maximum  adhesive  shear  strain 
retained  in  the  elastic  regime  is  as  follows: 


MOS  =  ^f — 1 
Yu* 

where  (89) 

MOS  >  0.2 

5.1.5  Determining  the  upper  bound  stress  intensity  factor  (Single  Sided 
Unsupported  Repairs) 

Using  the  values  of  cttp,  obp  and  Fx  calculated  from  section  5.1.3,  determine  the 
maximum,  minimum  and  root  mean  squared  upper  bound  stress  intensity 
factor  using  section  4.4.  The  margin  of  safety  of  safety  as  per  [3,  4,  5]  for  the 
upper  bound  stress  intensity  factor  is  as  follows: 

MOS  = -4s- — 1 
K 

where 
MOS  >  0.2 

5.1.6  Determining  the  Maximum  Peel  Stress  at  the  Crack  and  at  the  Edge 
of  the  Patch  (Single  Sided  Unsupported  Repairs) 

The  peel  stress  at  the  crack  and  at  the  edge  of  the  repair  can  be  determined 
using  section  4.6.  The  margin  of  safety  as  per  [3,  4,  5]  for  the  maximum 
adhesive  shear  strain  is  as  follows: 

cta  <  6000  psi  (Brittle  Adhesive  Systems) 
oa  — 10000  psi  (Ductile  Adhesive  Systems) 

5.1.7  Calculating  the  Maximum  Stress  in  the  Patch  (Single  Sided 

Unsupported  Repairs) 

It  is  assumed  that  all  the  stress  at  the  edge  of  the  repair  is  transferred  into  the 
patch.  This  is  valid  providing  the  actual  crack  length  is  greater  than  the 
characteristic  crack  length  (See  [8]  for  details).  This  results  in  the  following 
equation: 


\ 

(90) 


(91) 


The  margin  of  safety  as  per  [3,  4,  5]  for  the  maximum  patch  stress  is  as 
follows: 


MOS  =  1 

^ Max 

where  (92) 

MOS  >  0.2 

where,  ctulto  Patch  ultimate  strength 

5.2  Repair  Methodology  for  Single  Sided  Fully  Supported  Repairs 

5.2.1  Patch  Dimensions  (Single  Sided  Fully  Supported  Repairs) 

Determine  the  patch  dimensions  using  section  3.2.1. 

5.2.2  Adhesive  Management 

Check  the  Rapid  Repairability  Criterion  using  section  3.2.2.  The  margin  of 
safety  as  per  [3, 4, 5]  for  the  adhesive  load  capacity  is  as  follows^ 

P  >  1 .2  <?uLritj  v  (93) 

5.2.3  Stress  Under  and  at  the  Edge  of  the  Repair  (Fully  supported  single 
sided  repairs) 

•  Applied  Stress 

For  a  repair  to  fully  support  structure,  a  uniform  stress  will  result.  Therefore 
by  setting  the  curvatures  to  zero,  equation  58  in  section  4.3  reduces  to  the 
following: 


(94) 


Using  section  4.3  with  zero  curvatures  and  equation  94,  calculate  the  stress  in 
parent  structure  under  the  repair  in  the  x-direction  and  denote  its  value  as  a0. 
Calculate  the  total  force  in  the  parent  structure  at  the  edge  of  the  repair  in  the 
x-direction  using  section  4.3,  denoted  its  value  as  Fx. 

•  Thermal  Stresses 


Since  restraint  is  provided  against  out  of  plane  bending  no  curvatures  exist. 
Therefore  equation  36  in  section  4.2  reduces  to  equation  51: 


Using  section  4.2  with  zero  curvatures  and  equation  51,  calculate  the  stress  in 
the  parent  structure  under  the  repair  in  the  x-direction  and  denote  its  value  as 
a0.  Calculate  the  total  force  in  the  parent  structure  at  the  edge  of  the  repair  in 
the  x-direction  using  section  4.2,  denoted  its  value  as  Fx.  This  procedure 
should  be  completed  twice,  for  both  the  cure  and  operational  thermal  cycles. 
Then  add  the  corresponding  stresses  from  these  two  cycles  to  produce  the 
total  thermal  stress  components. 

5.2.4  Finding  the  Total  Stress  (Single  Sided  Fully  Supported  Repairs) 

Add  the  values  of  a0  and  Fx,  determined  from  the  applied  and  thermal  stress 
procedures.  The  margin  of  safety  of  safety  as  per  [3,  4,  5]  for  the  stress  in  the 
parent  structure  at  the  edge  of  the  repair  is  given  by  equation  88. 

5.2.5  Determining  the  Maximum  Adhesive  Shear  Strain  (Single  Sided 
Fully  Supported  Repairs) 

The  maximum  adhesive  shear  strain  can  be  determined  using  section  3.2.3. 
The  margin  of  safety  as  per  [3,  4,  5]  for  the  maximum  adhesive  shear  strain  is 
given  by  the  following: 

MOS  =  ^-l 

Ymox 

where  '  (95) 

MOS  >  0.2 

5.2.6  Determining  the  upper  bound  stress  intensity  factor  (Single  Sided 
Fully  Supported  Repairs) 

Calculate  the  upper  bound  stress  intensity  factor  using  section  3.2.4.  The 
margin  of  safety  of  safety  as  per  [3,  4,  5]  for  the  upper  bound  stress  intensity 
factor  is  given  by  equation  90. 


5.2.7  Determining  the  Maximum  Peel  Stress  at  the  Edge  of  the  Patch 
(Single  Sided  Fully  Supported  Repairs) 

The  peel  stress  at  the  edge  of  the  repair  can  be  determined  using  section  3.2.4. 
The  margin  of  safety  as  per  [3, 4,  5]  for  the  maximum  adhesive  shear  strain  is 
as  follows: 


<j a  <  6000  psi  (Brittle  Adhesive  Systems) 
ct a  <  10000  psi  (Ductile  Adhesive  Systems) 


5.2.8  Calculating  the  Maximum  Stress  in  the  Patch  (Single  Sided  Fully 
Supported  Repairs) 

It  is  assumed  that  all  the  stress  at  the  edge  of  the  repair  is  transferred  into  the 
patch.  This  is  valid  providing  the  actual  crack  length  is  greater  than  the 
characteristic  crack  length  (See  [8]  for  details).  Calculate  the  maximum  stress 
transferred  into  the  patch,  given  by  equation  91.  The  margin  of  safety  as  per 
[3, 4, 5]  for  the  maximum  patch  stress  is  given  by  equation  92. 


5.3  Repair  Methodology  for  Double  Sided  Repairs 

5.3.1  Patch  Dimensions  (Double  Sided  Repairs) 

Determine  the  patch  dimensions  using  section  3.2.1. 

5.3.2  Adhesive  Management 


Check  the  Rapid  Repairability  Criterion  using  section  3.2.2.  The  margin  of 
safety  as  per  [3, 4, 5]  for  the  adhesive  load  capacity  is  given  by  equation  93 

5.3.3  Stress  Under  and  at  the  Edge  of  the  Repair  (Double  sided  repairs) 

\ 

•  Applied  Stress 

Including  an  additional  patch  will  correspondingly  change  some  of  the 
inclusion  properties.  In  such  a  case,  equations  27  and  30  of  section  4.2  become: 


AI=—A  +  -^~  A. 


t  +  t  n 


t  +  tD 


a. 


u 


\a. 


■2i , 


t  +  t  jn 


(96) 

(97) 


Using  section  4.3  with  zero  curvatures,  noting  equations  94,  96  and  97, 
calculate  the  stress  in  parent  structure  under  the  repair  in  the  x-direction. 
Denote  its  value  as  a0.  Calculate  the  total  force  in  the  parent  structure  at  the 
edge  of  the  repair  in  the  x-direction  using  section  4.3,  denoting  its  value  as  Fx. 

•  Thermal  Stresses 

Using  section  4.2  with  zero  curvatures,  noting  equations  51,  96  and  97, 
calculate  the  stress  in  the  parent  structure  under  the  repair  in  the  x-direction. 
Denote  its  value  as  a0.  Calculate  the  total  force  in  the  parent  structure  at  the 
edge  of  the  repair  in  the  x-direction  using  section  4.2,  denoting  its  value  as  Fx. 
This  procedure  should  be  completed  twice/for  both  the  cure  and  operational 


thermal  cycles.  Then  add  the  corresponding  stresses  from  these  two  cycles  to 
produce  the  total  thermal  stress  components. 

5.3.4  Finding  the  Total  Stress  (Double  Sided  Repairs) 

Add  the  values  of  g0  and  Fx,  determined  from  the  applied  and  thermal  stress 
procedures.  The  margin  of  safety  of  safety  as  per  [3,  4,  5]  for  the  stress  in  the 
parent  structure  at  the  edge  of  the  repair  is  given  by  equation  88. 

5.3.5  Determining  the  Maximum  Adhesive  Shear  Strain  (Double  Sided 
Repairs) 

The  maximum  adhesive  shear  strain  can  be  determined  using  section  3.2.3. 
The  margin  of  safety  as  per  [3, 4,  5]  for  the  maximum  adhesive  shear  strain  is 
given  by  equation  95: 

5.2.6  Determining  the  upper  bound  stress  intensity  factor  (Double  Sided 
Repairs) 

Calculate  the  upper  bound  stress  intensity  factor  using  section  3.2.4.  The 
margin  of  safety  of  safety  as  per  [3,  4,  5]  for  the  upper  bound  stress  intensity 
factor  is  given  by  equation  90. 

5.3.7  Determining  the  Maximum  Peel  Stress  at  the  Edge  of  the  Patch 
(Double  Sided  Repairs) 

The  peel  stress  at  the  edge  of  the  repair  can  be  determined  using  section  3.2.4. 
The  margin  of  safety  as  per  [3,  4,  5]  for  the  maximum  adhesive  shear  strain  is 
as  follows: 


ga  ^  6000  psi  (Brittle  Adhesive  Systems) 
ga  <  10000  psi  (Ductile  Adhesive  Systems) 

5.3.8  Calculating  the  Maximum  Stress  in  the  Patch  (Double  Sided  Repairs) 

It  is  assumed  that  all  the  stress  at  the  edge  of  the  repair  is  transferred  into  the 
patch.  This  is  valid  providing  the  actual  crack  length  is  greater  than  the 
characteristic  crack  length  (See  [8]  for  details).  Calculate  the  maximum  stress 
transferred  through  the  patch,  given  by  equation  91.  The  margin  of  safety  as 
per  [3, 4, 5]  for  the  maximum  patch  stress  is  given  by  equation  92. 

5.4  Unsatisfactory  Margin  of  Safety 

The  margin  of  safeties  listed  in  the  previous  section  are  provided  to  ensure 
the  proposed  repair  design  is  appropriate  in  repairing  the  cracked  region 
efficiently.  The  inability  to  satisfy  any  of  these  criteria  may  lead  to  an 
ineffective  repair.  If  one  or  more  of  the  criteria  cannot  be  met  by  the  designer. 


there  are  a  number  of  options  available  to  improve  the  design.  C5033  and 
CalcuRep  contain  a  number  of  options  to  reduce  stresses  and  other  specified 
variables  in  order  to  meet  design  criteria.  The  following  is  a  list  of  advice 
given  by  C5033  and  CalcuRep  to  the  user. 

C5033 

Insufficient  overlap  length 

•  Use  a  higher  modulus  patch  material 
Insufficient  load  capacity 

•  Select  an  adhesive  with  a  higher  plastic  strain  capacity 

•  Match  repair  stiffness  with  parent  structure  stiffness 

•  Match  thermal  expansion  coefficients  of  repair  and  parent  structure 

•  Increase  the  adhesive  thickness 

If  the  upper  bound  stress  intensity  factor  is  too  high 

•  Use  an  adhesive  with  a  high  shear  modulus 

•  Increase  the  stiffness  of  the  patch  (Note:  do  note  exceed  1.2  Eiti) 

Peel  stress  too  high 

•  Taper  the  ends  of  the  repair 

•  Use  a  higher  modulus  patch  material 

•  Use  serrated  ply  ends 

CalcuRep 

If  the  upper  bound  stress  intensity  factor  is  too  high,  try  the  following: 

•  Use  a  thicker  or  shorter  patch 

•  Use  a  patch  material  with  a  higher  coefficient  of  thermal  expansion 
If  the  patch  stress  is  too  high 

•  Use  a  thicker  or  wider  patch 

•  Use  a  shorter  overlap  length 

•  Stronger  patch  material 

If  the  stress  in  the  parent  structure  at  the  edge  of  the  patch  is  too  high 

•  Use  a  thinner  or  wider  or  tapered  patch ' 

•  Use  a  shorter  overlap  length 


•  Use  a  patch  material  with  a  lower  coefficient  of  thermal  expansion 
If  the  adhesive  Shear  Strain  is  too  high 

•  Use  a  thicker  patch 

•  Use  a  patch  material  with  a  higher  coefficient  of  thermal  expansion 
If  the  load  transfer  is  too  high 

•  Use  a  longer  or  thinner  patch 

•  Use  a  less  stiff  patch 


6.  Case  Study  Comparison 

A  finite  element  model  case  study  was  developed  in  reference  [22]  and 
demonstrated  the  significance  of  out  of  plane  bending  effects  in  unsymmetric 
bonded  repairs.  The  results  produced  from  this  study  were  alarming  and 
suggested  that  the  stress  intensity  factor  for  a  single  sided  unsupported  repair 
increased  indefinitely  with  increasing  crack  length.  An  ^earlier  model 
developed  in  reference  [15]  supported  these  findings.  As  a  result,  serious 
doubt  was  cast  over  the  validity  of  existing  bonded  repair  technology. 
Therefore  an  independent  study  was  conducted  in  references  [10]  and  [11] 
using  this  case  study  to  confirm  or  otherwise  the  findings  produced  by 
references  [15]  and  [22].  The  results  of  the  independent  study  showed  that  the 
stress  intensity  factor  in  single  sided  unsupported  repairs  approaches  an 
asymptote  value  for  increasing  crack  length. 

The  case  study  used  above  is  a  standard  repair  case  consisting  of  a  boron 
patch  bonded  to  a  centre  cracked  aluminium  plate  with  a  uniaxial  stress 
applied.  The  patch  covers  the  entire  width  of  the  plate,  therefore  load 
attraction  effects  are  neglected.  The  repair  case  configuration  and  material 
properties  are  listed  in  section  6.1. 

Verification  of  previous  single  sided  supported  and  double-sided  repair 
analysis  tools  used  by  C5033  and  CalcuRep  has  been  provided  in  reference 
[8].  Verification  of  analyse  tools  listed  in  sections  4.2  [12],  4.4  [10, 11],  4.5  [11], 
4.6  [11]  and  4.7  [21]  have  previously  been  confirmed  by  the  references  shown 
above.  Complete  verification  of  the  analysis  tools  shown  in  section  4.3  have 
yet  to  be  provided  and  will  be  demonstrated  on  release  of  reference  [13].  In 
the  mean  time,  a  preliminary  verification  process  will  be  provided  by  the  use 
of  the  finite  element  cases  study  discussed  previously. 

The  following  section  is  denoted  as  a  preliminary  verification  process  for  the 
following  two  reasons: 


(1)  The  case  study  used  in  reference  [8]  and  [9]  excludes  load  attraction 
effects.  The  procedure  listed  in  section  5.1  (which  includes  the  use  of 
section  4.3)  includes  load  attraction  effects  and  cannot  be  modified  to 
exclude  these  effects. 

(2)  Comparing  the  results  produced  from  section  5.1  with  those  obtained 
from  reference  [8]  and  [9]  will  not  provide  a  completely  fair  comparison 
due  to  the  point  above.  It  will  rather  demonstrate  whether  the  results 
approximately  coincide  and  how  the  differences  listed  in  point  1  above 
effect  the  results. 


6.1  Single  Sided  Unsupported  Repair  Configuration 


Part 

Elastic 

Modulus 

(MPa) 

Poisson's 

Ratio 

Thicknes 

s 

(m) 

Length 

(m) 

Width 

(m) 

Plate 

71000 

0.3 

0.003 

0.5 

0.5 

Repair 

207000 

0.3 

0.00102 

0.1 

0.5 

Adhesive 

1890 

0.33 

0.000203 

0.1 

0.5 

(Table  1) 


The  above  repair  configuration  is  a  single  sided  unsupported  repair  with  an 
applied  stress  at  the  plate  boundary.  No  thermal  effects  have  been  taken  into 
account.  The  load  attraction  can  be  assumed  to  be  zero,  since  the  patch  covers 
the  entire  width  of  the  plate.  Reference  [8]  and  [9]  provide  results  for  the 
upper  bound  minimum,  root  mean  squared  and  maximum  stress  intensity 
factors  due  to  an  applied  stress. 


6.2  Results 


The  methodologies  discussed  in  section  5  of  this  report  have  been 
programmed  into  Mathcad  [9].  The  Mathcad  output  listed  in  the  appendices 
are  as  follows:  Single  Sided  Unsupported  repairs  section  5.1  (Al),  Single 
Sided  Fully  supported  repairs  section  5.2  (A2),  Double  Sided  repairs  section 
5.3  (A3),  Single  Sided  Unsupported  repairs  section  5.1  (A4),  Stiffened  Panels 
with  a  crack  in  either  the  stiffener  or  the  skin  beneath  the  stiffener  section  5.4 
(A5). 

The  appendices  listed  above  detail  the  following: 

Al  Section  5.1  applying  the  case  study  above. 

A2  to  A5  The  remaining  appendices  relay  examples  of  the  methods  listed 
in  sections  5.2  to  5.4  for  the  F-lll  Boron  repair  discussed  in 
reference  [8].  This  is  provided  for  two  reasons;  1)  to  allow  the 
user  to  view  how  the  mathematical  expressions  in  section  5  are 
arranged,  2)  to  provide  a  means  to  verify  any  future  software 
developments  against  the  methods  discussed  in  section  5. 


Using  the  single  sided  unsupported  repair  design  mathcad  module  (ie  Section 
5.1  Appendix  Al)  and  running  the  case  study  produced  the  following  results. 


Reference 

KWa“  (Vm) 

K'rms/ct”  (Vm) 

KWa“  (Vm) 

[10&11] 

0.102 

0.282 

Section  5.1 

0.125 

0.261 

| 

Table  2 


Relative  Asymptote  Stress  Intensity  Factor  Comparison  Post  Repair 
Including  and  Excluding  Load  Attraction 


Figure  9 


7.  Main  Discussion  ' 

Applying  the  case  study  using  section  5.1,  the  minimum,  root  mean  squared 
and  maximum  stress  intensity  factors  were  calculated.  The  results  are  shown 
in  table  2  and  figure  9  above.  Figure  9  demonstrates  that  load  attraction 
reduces  the  root  mean  squared  and  maximum  stress  intensity  factors.  This 
phenomenon  can  be  explained  by  considering  the  following  equation  from 
[13]: 


Rearranging  equation  98  to  make  the  mean  strains  in  the  inclusion  the  subject 
gives  the  following: 

r  x 

<£x  >  =  L?  R2-L4  j*'j  (99) 

Applying  a  stress  to  the  boundary  of  an  unsupported  single  sided  repair  will 
cause  a  positive  curvature  inside  the  inclusion  section  (See  figure  6).  This  is 
predominantly  due  to  the  neutral  axis  of  the  plate  shifting  towards  the  neutral 
axis  of  the  inclusion  section.  With  this  concept  in  mind,  consider  equation  99 


for  positive  curvatures  and  positive  U.  A  positive  curvature  in  the  inclusion 
section  will  result  in  a  reduction  of  the  mean  strains  inside  the  inclusion.  If  the 
curvatures  are  large  enough,  as  in  this  case  study,  load  is  extracted  from  the 
repaired  region  rather  than  attracted  into  the  repaired  region. 

In  establishing  the  above  point,  now  consider  the  strain  energy  release  rate 
provided  by  equation  43  of  reference  [10],  listed  as  follows: 


(q<t  ^(l  ^12)^2 
(1  +  S)SfiE, 


(100) 


Where, 


G'»  Upper  bound  strain  energy  release  rate 

a”  Applied  stress 

S  Stiffness  Ratio 

P  Inverse  of  the  load  transfer  length 

Ep  Plate  Modulus 

co  Bending  correction  factor 

Q  Load  attraction  factor  (Not  included  in  reference  [10]) 


Note:  The  upper  bound  strain  energy  release  rate  provided  in  section  4.4  is 
derived  considering  crack  bridging  by  springs  and  compares  well  with 
equation  100.  The  above  equation  is  used  in  this  case  to  provide  a  simpler 
mean  to  justify  the  following  argument,  rather  than  using  the  more  complex 
form  of  equation  68  in  section  4.4. 

The  upper  bound  root  mean  squared  stress  intensity  factor  is  therefore  given 
by  equation  79  as  following: 


K 


Lc c.RMS 

or 


K. 


=  £lca 

& " 


(1  +  S)SJ3 


(101) 


Equation  101  demonstrates  that  a  decrease  in  load  attraction  will  result  in  a 
decrease  in  the  upper  bound  root  mean  square  stress  intensity  factor,  thus 
confirming  the  load  extraction  effect  previously  discussed  and  shown  figure 
9. 

Consider  the  gradients  of  the  two  curves  plotted  in  figure  9.  The  difference  in 
the  gradient  between  the  two  results  is  predominantly  due  to  the  'R'  ratio 
effect  shown  in  equation  81  and  82.  For  larger  'R'  values  (See  equation  80)  the 
maximum  upper  bound  stress  intensity  factor  will  decrease  and  the  minimum 
upper  bound  stress  intensity  factor  will  increase.  The  actual  'R'  value 


calculated  by  reference  [10]  and  section  5.1  are  0.236  and  0.333  respectively 
and  confirms  the  gradient  trend  shown  in  figure  9. 

Consider  the  dominating  factors  effecting  the  'R'  ratio,  equation  80: 


(102) 


For  larger  mean  to  bending  stress  ratios  (— ),  a  larger  'R'  value  is  produced. 

The  mean  to  bending  stress  ratios  produced  by  reference  [10]  and  section  5.1 
are  1.62  and  2  respectively.  From  these  two  figures  it  can  be  seen  that  the  load 
extraction  from  section  5.1  has  significantly  reduced  maximum  bending  stress 
(ob).  This  consequently  results  in  a  lar'ger  'R'  value  and  is  the  predominant 
factor  producing  the  gradient  differences  in  figure  9. 

Consider  the  effects  of  thermal  stresses  in  the  analysis.  Conclusions  from 
reference  [12]  demonstrate  the  following: 

•  Load  is  attracted  into  the  repaired  region 

•  The  stress  at  the  top  of  the  plate  under  the  repair  has  increased  from  the 
mean  value. 

•  The  stress  at  the  bottom  of  the  plate  under  the  repair  has  decreased  from 
the  mean  value 

The  above  results  have  demonstrated  the  following  for  an  unsupported  repair 
with  a  uniaxial  load  applied: 

•  Load  is  extracted  from  the  repaired  region 

•  The  stress  at  the  top  of  the  plate  under  the  repair  has  decreased  from  the 
mean  value. 

•  The  stress  at  the  bottom  of  the  plate  under  the  repair  has  increased  from 
the  mean  value 

These  points  show  that  out  of  plane  bending  effects  produced  by  applied 
stresses  counteract  out  of  plane  bending  effects  caused  by  thermal  residual 
stresses  in  single  sided  unsupported  repairs.  Further  work  would  be  required 
to  establish  which  of  these  two  applications  may  be  dominant  and  under 
which  conditions. 


8.  CONCLUSION  AND  RECOMMENDATIONS 


Bonded  repair  technology  is  widely  becoming  the  preferred  method  of  repair 
over  mechanically  fastened  repairs.  An  efficient  and  effective  repair  requires 
the  use  of  validated  analytical  formulations  that  describe  events  occurring  in 
reality.  Such  analytical  tools  exist  and  are  becoming  more  developed  as  the 
technology  becomes  the  preferred  repair  method.  A  summary  of  these 
existing  tools  as  used  by  C5033  and  CalcuRep  has  been  provided  in  CalcuRep 
review  task  1  [8]. 

CalcuRep  review  task  2  has  been  centred  toward  developing  a  new  approach 
that  includes  the  best  tools  from  both  CalcuRep  and  C5033,  as  well  as  utilising 
the  most  recent  developments  in  this  area.  The  key  to  extending  CalcuRep' s 
capability  has  been  to  logically  arrange  existing  and  new  design  tools  to 
develop  an  appropriate  repair  design  approach  for  a  wider  range  of  repair 
cases.  The  methods  used  for  this  procedure  are  listed  in  section  5  of  this 
report.  This  report  has  been  successful  in  providing  uses  with  a  new  approach 
that  deals  with  the  following  repair  scenarios: 

•  Single  Sided  Fully  supported  repairs  , 

•  Double  Sided  repairs 

•  Single  Sided  Unsupported  repairs 

•  Stiffened  Panels  with  a  crack  in  either  the  stiffener  or  the  skin  beneath  the 
stiffener 

This  report  has  also  provided  mathcad  modules  in  appendices  A1  to  A5 
which  detail  the  new  approach  (section  5)  for  each  case  shown  above. 
Designers  using  these  approaches  can  either  use  the  mathcad  modules 
provided  in  the  appendices  or  program  a  software  package  ensuring 
verification  against  these  modules. 

Section  6  of  this  report  details  a  case  study,  which  was  used  to  provide  a 
preliminary  validation  of  the  approach  developed  in  section  4.3.  This  case 
study  consisted  of  an  unsupported  single  sided  repair  with  a  uniaxial  stress 
applied  to  the  plate  boundary.  The  results  of  this  study  produced  the  upper 
bound  stress  intensity  factor  post  repair  excluding  load  attraction  effects  by 
the  use  of  Finite  Element  Analysis  [10, 11].  Using  the  new  approach  described 
in  section  5.1  for  this  particular  repair  case,  corresponding  upper  bound  stress 
intensity  factors,  taking  into  account  load  attraction  affects,  were  calculated. 
Discrepancies  in  the  results  were  directly  due  to  the  inclusion  of  load 
attraction  effects,  which  reduced  the  root  mean  and  maximum  upper  bound 
stress  intensity  factors.  Reductions  of  this  nature  were  found  to  be 
predominantly  due  to  load  being  extracted  from  the  repaired  region,  rather 
than  attracted  into  the  repair  region.  The  findings  also  demonstrated  that  out 
of  plane  bending  stresses  produced  by  -thermal  residual  stresses  can  be 
beneficial  in  counteracting  those  produced  by  applied  stress. 


Further  work  can  be  developed  on  the  following  points: 

•  Complete  verification  of  the  approach  listed  in  section  5.1  has  yet  to  be 
performed  by  a  finite  element  comparison  allowing  for  load  attraction 
effects.  This  can  be  done  as  follows:  Develop  a  2-stage  finite  element 
model  of  an  unsupported  finite  circular  plate  bonded  by  a  finite  circular 
orthotropic  patch.  The  first  stage  of  the  analysis  is  to  determine  the 
stresses  in  the  repair  region  including  thermal  and  applied  stress  effects 
(Biaxial)  assuming  no  crack  present.  Compare  these  results  against  the 
method  described  in  section  5.1.  Introduce  a  centre  crack  in  the  finite 
element  model  and  determine  the  upper  bound  minimum,  root  mean 
squared  and  maximum  stress  intensity  factor.  Compare  against  results 
from  section  5.1. 

•  Confirm  the  approach  listed  in  section  5.1  and  the  finite  element  analysis 
by  the  use  of  experimental  test  data. 

•  Expand  the  analysis  to  consider  an  elliptical  patch  and  plate.  This  allows 
the  user  to  adjust  the  patch  and  plate  length  and  width  rather  than 
corresponding  radii,  giving  the  user  more  scope  in  the  design  method. 

•  Develop  a  software  package  that  includes  the  recommended  approaches 
listed  in  section  5  of  this  report. 

•  Determine  whether  out  of  plane  bending  effects  due  to  applied  stresses  are 
more  dominant  than  those  produced  by  thermal  stresses,  or  visa  versa. 
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SUMMARY 

There  are  several  methods  available  for  obtaining  a  first  estimate  for  the  design  of 
composite  repairs  to  cracks  in  thin  metallic  skins.  The  program  Calcurep  is  based  one 
such  technique,  which  builds  on  an  analogy  with  the  problem  of  a  bonded  symmetric 
lap  joint. 

The  present  report  begins  by  comparing  the  results  of  a  3D  finite  element  analysis  of 
various  composite  repairs  with  these  analytical  estimates.  Finite  element  studies  of 
composite  repairs  to  (cracked)  rib  stiffened  panels  are  then  presented.  This  study 
reveals  the  stress  intensity  factors  to  be  a  weak  function  of  the  stiffener  spacing.  The 
current  analytical  formulae  are  then  extended  to  allow  for  repairs  to  rib  stiffened 
panels  and  verified  by  comparison  with  the  results  of  a  detailed  experimental  test 
program. 

To  complete  this  section  a  simple  test  methodology,  based  on  the  ASTM  thick 
adherend  test,  ASTM  D-1002,  for  determining  the  true  visco-plastic  behaviour  of 
structural  film  adhesives,  is  also  presented.  A  methodology  for  accounting  for  the 
observed  visco-plasticity  is  then  presented. 

The  analytical  formulae  make  extensive  use  of  the  analytical  solution  for  a  simple  lap 
joint.  Unfortunately,  this  solution  is  restricted  to  uniformly  thick  adhesives, 
monotonic  loading  and  an  idealised  inelastic  response.  To  extend  the  range  of 
applicability  of  these  formulae  a  range  of  new  solutions  are  developed  that  include: 

i)  variable  adhesive  thicknesses; 

ii)  random  'variations  in  the  adhesive  thickness; 

iii)  more  representative  adhesive  behaviour,  e.g.  visco-plastic  response. 

As  a  result  of  this  work  we  find  that,  even  for  extensive  plasticity,  even  if  there  is 
variability  in  the  thickness  of  the  adhesive  bond  the  energy  field  Wmax  can  be 
computed  using  only  the  elastic  solution  for  the  repair.  Consequently,  we  show  that 
for  repairs  with  significant  variability  in  the  thickness  of  the  adhesive  bond  the 
strength  of  the  repair,  provided  failure  is  due  to  failure  of  the  adhesive,  can  be  then 
estimated  from  a  purely  linear  elastic  analysis.  This  will  occur  at  the  load  level  when 
Wmax  =  WCrit .  This  finding  significantly  simplifies  the  design  process  for  repairs  with 
a  variable  thickness  adhesive. 

It  is  also  mentioned  that  the  majority  of  commercially  available  finite  element 
programs,  e.g.  ABAQUS,  MSC-Nastran,  NEnastran,  MARC,  ANSYS  etc,  use  iso¬ 
parametric  elements.  When  analysing  three  dimensional  problems  involving 
composites  the  use  of  elements  with  mid-side  nodes  is  essential,  i.e.  20  noded  hex 
(brick)  elements,  15  noded  wedge  elements,  etc.  Unfortunately,  ABAQUS,  MSC- 
Nastran,  MARC,  ANSYS  use  formulations  based  on  that  given  in  Zienkiewicz  [23]. 
This  means  that  for  accurate  results  the  aspect  ratio’s  are  limited  to  less  than  —5:1, 
even  if  reduced  integration,  i.e.  the  use  of  2  x  2  x  2  Gauss  points  in  formulating  the 
element  stiffness  matrix,  is  used.  In  contrast  NEnastran  use  an  enriched  20  noded 
element,  where  the  basic  formulation  is  extended  to  include  additional  shape 
functions.  This  means  that  the  NEnastran  3D  (mid-side  nodded)  elements  can  be  used 
with  dramatically  greater  aspect  ratio’s  and  as  such  are  better  suited  to  the  analysis  of 
composite  repairs. 
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1.  NUMERICAL  ANALYSIS  &  EXTENSION  OF  2D  FORMULAE  TO 
ALLOW  FOR  3D  EFFECTS 

Externally  bonded  composite  patches  have  proved  to  be  an  effective  method  of 
repairing  cracked,  or  damaged,  structural  components,  [1-3,  6-12,  16,  17,  19].  A 
variety  of  approaches  are  now  available  for  the  design  of  composite  repairs  to  cracks 
in  thin  metallic  skins,  i.e.  typical  thickness  less  than  ~  3  mm.  One  such  technique  is 
based  on  the  fact  that  the  stress  intensity  factor  for  a  patched  crack  approaches  a 
constant  (limiting)  value,  defined  as  Koo,  as  the  crack  length  increases.  This  approach 
was  based  on  the  premise  that,  for  a  sufficiently  long  crack  in  a  structure,  which  is 
subjected  to  a  remote  uniform  stress  field,  the  central  region  of  the  patch,  over  the 
crack,  behaves  like  an  overlap  joint  [3,  2],  From  this  premise  it  follows  that  the  stress 
distribution  in  this  central  region  and  the  stress  intensity  factor  should  become 
independent  of  crack  length  . 

As  a  result  of  this  analogy  it  has  been  found  that  the  problem  of  a  bonded  symmetric 
lap  joint  can  be  used  in  the  initial  design  process.  Indeed,  the  resultant  analytical 
formulae  are  particularly  easy  to  use  and  provide  a  first  estimate  for  the  patch  design. 

It  is  possible  to  increase  the  accuracy  of  the  current  approximate  2D  formulae, 
initially  developed  by  Rose  in  [2,  3],  for  the  limiting  stress  intensity  factor  K«>  for  a 
crack  repaired  with  an  externally  bonded  composite  repair,  by  (partially)  accounting 
for  through-the-thickness  effects.  To  this  end  the  value  of  Kq  can  be  approximated  by 
the  formulae; 

\ 

Koo  =  Y  ClLa0^fnl  (1) 

\ 

where  g0  =  aEptp/(Eptp+Ertr)  (2) 

Y  is  a  geometry  factor,  which  accounts  for  repairs  to  centre  or  edge  cracks; 

Y  is  a  geometry  factor  =1  for  a  repair  to  a  centre  crack 

=  0.9  for  a  repair  to  an  edge  crack 

and  Ql  is  the  load  attraction  factor.  For  long  uni-directional  fibre  patches  it  has  been 


found  that  the  term  QL  can  be  approximated  as  follows 

r=3/8  1  ./(1+X  E,-  tr/(B  Ep  tp))-  Ep  tp  /(  Er  tr  +  EP  tP))  (3) 

QL  =  (Eptp/(Ertr  +  Eptp)  +  r)/CT0  (4) 

TtX.  =  ^Eptp/(3(1  +tpEp/Ertr)  (5) 

and 

P  =  (ta/Ga+tr/3Gr+tp/3Gp)/(ta/Ga+3tr/8Gr  +  3tp/8Gp)2  (6) 


Here  X  is  the  width  of  the  patch,  B  is  the  width  of  the  plate/skin,  ta,  tp  and  tr  are  the 
thickness’  of  the  adhesive,  plate,  and  patch  respectively,  G  and  E  denote  the  shear  and 
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Young’s  modulus  and  the  subscripts  a,  p  and  r  denote  their  values  for  the  adhesive, 
plate,  and  patch  respectively. 

This  expression1  has  been  found  to  produce  results  in  good  agreement  with  three 
dimensional  finite  element  analysis,  see  Section  1.1.  The  program  utilising  this 
approximation  will  be  termed  CalcuRep*. 

Until  now  we  have  ignored  bending  effects.  However,  even  if  the  wing  skin  is  in  a 
state  of  plane  stress,  the  location  of  the  neutral  axis  of  the  patch-adhesive-skin  section 
will  differ  from  the  neutral  axis  of  the  wing  skin  itself.  Hence,  forces  applied  to  the 
skin  will  result  in  an  out  of  plane  bending,  which  will  reduce  the  efficiency  of  the 
repair. 

There  are  several  methods  that  can  be  used.  One  approach  is  presented  in  Chapter  1 . 
The  method  developed  at  Northrop,  see  [14]  for  more  details,  can  also  be  used  to 
account  for  this  out  of  plane  bending.  In  this  method,  the  apparent  stress  intensity 

factor  K*  at  the  mid  surface  of  the  sheet  is  given  by: 

Kp  =  (1  +BC)  Kp  (7) 

where  BC  is  a  bending  correction  factor.  Here, 

BC  =  a  ymax  (1  -  Kp/ Ks)  (tp  +  tr)II  (8) 

where  Ks  is  the  value  of  the  stress  intensity  factor  before  patching,  tp and  tR  are  the 
thicknesses  of  the  sheet  and  patch  respectively,  ymax  is  the  distance  of  the  lower 
unpatched  surface 'of  the  plate  from  the  neutral  axis  of  the  section  (i.e.,  sheet  plus 
patch).  I  is  the  moment  of  inertia  of  the  section  and  a  is  the  crack  half  length. 

The  value  of  K*p  can  be  related  to  J,  the  energy  release  rate  for  self  similar  crack 
growth,  in  the  usual  way,  viz: 

J  =  Kf/E  (9) 

When  growth  is  non  self-similar  with  the  maximum  growth  occurring  at  the  lower 

H/t  *  , 

free  surface,  it  is  best  not  to  use  Kp  or  J  but  to  design  on  the  basis  of  the  maximum 
stress  intensity  factor  Kp™  which  is  given  by: 

Kmax  =  (1  +2BC)Kp  (10) 

1.1  Repair  of  thin  skins 

To  illustrate  the  accuracy  of  these  formulae  let  us  first  consider  an  externally  bonded 
composite  repair  (patch)  to  a  centre  crack  in  a  3mm  thick  aluminium  panel  with  E= 
72,000  MPa  and  v  =  0.33.  Two  patch  thicknesses  were  analysed  viz:  1.0  and  1.85 


1  The  derivation  is  both  cumbersome  and  tedious  and  hence  for  simplicity  will  be 
omitted. 
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mm,  and  the  patch  was  assumed  have  the  following  mechanical  properties:  En 
=208,000  MPa,  E22  =  25432  MPa,  v12  =  0.183,  and  Gi2  =  Gi3  =  G23  =  7241  MPa.  The 
plate  was  assumed  to  have  dimension  of  200  (length)  x  290  (width)  x  3  mm  and  the 
patch  was  assumed  to  have  a  planform  of  100  (length)  mm  x  82  (width)  mm.  The 
adhesive  was  taken  to  be  0.25  mm  thick  with  a  shear  modulus  of  375  MPa  and  v  = 
0.33.  The  analysis  used  20  noded  isoparametric  3D  brick  elements  with  2  layers  of 
elements  through  each  of  the  plate,  patch  and  adhesive.  The  plate  was  subjected  to  a 
remote  uniform  stress  of  229.8  MPa  and  global  bending,  after  patching,  was 
prohibited. 

In  this  analysis  symmetry  was  used  and  only  a  quarter  of  the  structure  was  modeled. 
The  resultant  mesh  consisted  of  4,889  nodes  and  1132  elements  and,  to  simulate  the 
crack  tip  singularity,  the  near  tip  elements  had  the  mid  side  nodes  moved  to  the  14 
points.  The  resultant  solution  was  well  conditioned  and  the  accuracy  of  the  solution 
was  evaluated  by  performing  two  separate  analyses.  The  first  used  optimum 
integration,  i.e.  2x2x2  Gaussian  quadradture  points,  to  form  the  stiffness  matrices 
of  the  elements  whilst  the  second  used  full  integration,  i.e.  3x3x3  Gaussian 
quadradture  points,  to  form  the  stiffness  matrix.  Both  analyses  gave  results,  which 
agreed  to  within  1%. 

Following  this  analysis  the  composite  repair  to  edge  notch  cracks  in  a  200  (length)  x 
145  (width)  x  3  mm  plate  was  then  considered.  In  this  study  the  patch  was  assumed  to 
be  1.85  mm  thick  and  have  a  plan  form  of  100  (length)  mm  x  41  (width)  mm.  The 
results  both  these  repair  configurations  are  presented  in  Tables  1.1 -1.4. 

\ 

The  third  study  involved  similar  repair  configurations  where  skin  thicknesses  of  2mm, 

3.15  and  adhesive  shear  modulii  of  700  MPa  were  also  considered. 

\ 

For  the  first  two  repair  configurations  the  differences  between  the  present 
approximation,  which  we  refer  to  as  Calcurep*,  and  those  used  in  C5033  and  Calcurep 
are  summarised  in  Tables  1.4.  The  differences  for  the  third  study  are  summarised  in 
Table  1.5. 

The  C5033  solution  assumed  equal  stresses  in  the  patch  and  the  plate.  This  was  not 
the  case.  The  (corrected)  CalcuRep+  solution  underestimated  the  solution  by  ~  18%  . 
Accounting  for  tbrough-the-thickness  effects  CalcuRep*  reduced  this  error  to  within 
the  error  bound  of  the  finite  element  solution,  see  Table  1.4.  Although  the  error  in  the 
C5033  solution  was  similar  to  that  of  CalcuRep+  its  other  assumptions  were  further 
from  the  true  values.  For  comparison  the  Calcurep’  results  for  repairs  to  centre  notch 
and  edge  notch  panels  are  shown  in  Table  1.5  along  with  the  corresponding  finite 
element  solutions.  In  each  case  the  agreement  is  shown  to  be  quite  good. 


Table  1 . 1 :  K  for  a  repaired  edge  notch  panel,  patch  thickness  1 ,85mm. 


Crack 

Length 

5  mm 

6  mm 

7  mm 

10  mm 

20  mm 

patched 

surface 

10.6 

10.8 

10.9 

10.9 

10.4 

middle 

11.6 

11.7 

11.9 

12.1 

12.9 

bottom 

11.7 

11.9 

12.0 

12.2 

12.7 
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Table  1.2:  K  values  for  repaired  centre  notch  panel,  patch  thickness  1 .85mm 


Crack  half 
Length 

5  mm 

6  mm 

7  mm 

10  mm 

20  mm 

patched 

surface 

11.2 

11.5 

11.8 

12.0 

11.6 

middle 

11.9 

12.4 

12.7 

13.3 

14.3 

bottom 

12.0 

13.16 

12.9 

13.4 

14.1 

Table  1.3:  K  for  a  repaired  centre  notch  panels  patch  thickness  1  mm. 


Crack  half 
Length 

5  mm 

6  mm 

7  mm 

10  mm 

20  mm 

patched 

surface 

14.6 

15.2 

15.6 

16.1 

15.7 

middle 

15.5 

16.3 

16.8 

17.8 

19.4 

bottom 

15.6 

16.4 

17.0 

17.9 

19.1 

Table  1.4:  Comparison  of  solutions  for  for  a  patch  thickness  1 .85  mm 


Crack 

Configuration 

Crack 

length 

C5033 

CalcurepT 

Calcurep 

K 

Finite  element 

Centre  crack 

20  mm 

17.4 

11.8 

14.1 

-14.1  to  14.3 

Edge  crack 

20  mm 

13.8 

10.3 

12.7 

-12.7  to  12.9 

T  Initial  corrected 

and  modi: 

Bed  version  of  CalcuRep 

CalcuRep*  incorporates  through  thickness  effects. 


Table  1.5:  Comparison  of  prediction  for  K«  for  various  patch  configurations 


Crack 

Configuration  ' 

Crack 

length 

Configuration 

Calcurep 

K 

Finite  element 

Edge  crack 
cj  =  344.8  MPa 

20  mm 

SkinA  2  mm, 
patch  -=  0.75mm 

Ga  =  375  MPa 
ta  =0.25  mm 

21.6 

17.8  -22.4 

Centre  crack 
a  =  344.8  MPa 

20  mm 

SkinA  2  mm, 
patch  -=  0.75mm 

Ga  =  375  MPa 
ta  =0.25  mm 

23.6 

19.8  -24.0 

Centre  crack 
a  =  229.8  MPa 

20  mm 

SkinA  3.  mm, 
patch  =1.00  mm, 

Ga  =  375  MPa 
tj  =0.25  mm 

19.2 

15.7-19.1 

Centre  crack 
a  =  229.8  MPa 

+  ^  -  - : — r~" 

20  mm 

i  i  •  r*  i 

SkinB  3.  mm, 
patch  =  0.889mm, 

Ga  =  750  MPa 
ta  =0.165  mm 

•  r  r' t  _  i 

16.2 

15.7-15.8 

*  CalcuRep*  incorporates  through  thickness  effects. 

A  The  planform  of  the  skin  and  patch  were  as  previously  described. 

B  Planform  of  plate  was  320  mm  x  150  mm  (wide)  with  a  semi-circular  patch,  radius 
=  80mm 
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The  analytical  formulae  are  also  quite  accurate  for  the  case  when  bending  is  allowed, 
see  Tables  1.6,  1.7  and  1.8.  However,  in  this  case  the  accuracy  of  the  simple  analytical 
formulae  decreases  when  the  crack  length  approaches  a  half  of  the  patch  width  (X). 


Table  1.6:  Comparison  of  prediction  for  bending  problems,  Case  1. 


Table  1.7:  Comparison  of  prediction  for  bending  problems,  Case  2. 


Edge  20  |  2  72000  |  0.75  j  208030  '  .25  '  375  345  |  200  [  200 

predictd  85.2  _ 


f.e.  |  83.4  _ _ 

_ T able  1.8:  Comparison  of  prediction  for  bending  problems,  Case  3 .  _ _ 

Kmax  Edge  a  tp  Ep  tr  Gr  h  Ga  a  X  B 

MPa  Vm  Centre  mm  mm  MPa  mm  MPa  mm  MPa  MPa  mm  mm 


"Centre  |  20  |3  71000  0.89  [  208000  .17  750  ^  230  41  145 

predictd  62.7 


2.  DEVELOPMENTS  IN  THE  STRESS  ANALYSIS  OF  BONDED 

REPAIRS 

2.1  Background . 

When  designing  bonded  repairs,  or  adhesively  bonded  joints,  the  stress/strain 
response  of  the  adhesive  plays  a  central  role  in  determining  both  the  load  carrying 
capacity  and  the  fatigue  performance  of  the  repair  or  joint,  see  [4,  9].  To  this  end  we 
need  a  test  methodology  capable  of  accurately  and  consistently  reproducing  this 
behaviour. 

In  this  context  Chiu  et  al.  [13]  have  developed  a  variant  of  the  ASTM  thick-adherend 
short  over-lap  adherend  test  specimen  for  characterising  the  stress/shear  behaviour  of 
thin  film  adhesives.  The  results  revealed  that  the  properties  of  the  film  adhesive  FM73 
exhibited  significant  visco-plasticity,  even  at  room  temperature.  Chiu  et  al.  [13]  also 
used  a  unified  constitutive  model  to  describe  the  visco-plastic  behaviour  of  the 
adhesive. 

This  section  also  assesses  the  structural  significance  of  this  visco-plastic  behaviour,  in 
a  realistic  symmetric  double  lap  joint,  via  a  detailed  finite  element  analysis.  Here  we 
show  that  the  stress/strain  behaviour  of  the  adhesive  in  a  bonded  joint  is  dependent  on 
both  the  loading  rate  and  the  load  history.  When  the  joint  is  subjected  to  monotonic 
loading  with  a  constant  loading  rate  the  energy  density,  which  is  the  primary  adhesive 
design  variable  in  a  bonded  repair,  at  the  worst  point  in  the  joint  is  essentially  time 
independent.  This  is  very  important  as  it  allows  the  existing  design  tools  to  be  used 
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even  though  the  adhesive  exhibits  significant  visco-plastic  behaviour,  and  as  a  result 
is  strongly  rate  dependent. 


Unfortunately,  the  peak  adhesive  shear  stresses  and  strains  are  strongly  effected  by  the 
loading  history.  Since  the  shear  stresses  are  continuous  across  the  interface  an 
accurate  knowledge  of  the  interlaminar  stresses  in  the  composite  repair  requires  these 
visco-plastic  effects  to  be  incorporated  in  the  analysis.  To  this  end  a  simple  fortran 
program  has  been  developed,  which  uses  the  Glinka  hypothesis  [14],  and  negates  the 
need  for  a  complex  non-linear  finite  element  analysis. 

2.2  Structural  Response  Of  Thin  Film  Adhesives 

As  has  been  previously  mentioned  the  stress/strain  behaviour  of  the  adhesive  plays  a 
central  role  in  determining  both  the  load  carrying  capacity  and  the  fatigue 
performance  of  a  bonded  repair,  see  [4,  9,  10,  11].  This  section  presents  a 
methodology  for  determining  the  stress/strain  response  using  a  thick  adherend  test 
specimen  and  illustrates  it  by  considering  the  commonly  used  adhesive  FM73. 

2.2.1  Specimen  preparation 

The  results  presented  in  this  section  were  obtained  using  the  symmetric,  thick 
adherend,  short  overlap,  shear  test  specimen.  The  dimensions  of  the  lap  shear 
specimen  are  shown  in  Figure  2.2.  These  specimens  were  made  from  hardened  tool 
steel. 

The  bond  surfaces  of  the  specimen,  after  being  grit  blasted,  were  immersed  in  a  1.0% 

silane  solution  for  about  10  minutes.  They  were  then  dried  in  an  oven  at  80°C.  The 
dried  specimens  were  subsequently  assembled  with  a  FM73  film  adhesive  at  room 
temperature.  These  assembled  specimens  were  then  vacuum  bagged  and  placed  in  an 
autoclave  where  the  adhesive  was  cured  under  a  pressure  of  105  kPa  and  at  a 

temperature  of  120°C.  The  bond  line  thickness  of  this  joint  was  approximately  0.2 
mm  (i.e.  this  is  a  typical  adhesive  thickness  for  2  layers  of  FM73). 


▲ 

19  mm 

T 

p 1  1 '  1  rJ 

- 

^Adhesive  ~  0.2  mm  thick 

190  mm 

◄ — -  - ► 

Figure  2.2  Schematic  diagram  of  the  ASTM  D  1002  thick  adherend  test  specimen. 


2.2.2  Experimental  procedure 

The  specimens  were  loaded  under  strain  control,  in  a  lOkN  Instron  test  machine  and 
the  shear  strain  was  measured  using  a  pair  of  linear  variable  displacement  transducers 
(LVDT).  A  clip  gauge  was  attached  to  the  LVDT  to  provide  a  feedback  signal  for  the 
control  unit  in  the  Instron.  This  clip-gauge  was  necessary  because  the  control  signal 
from  the  LVDT  was  too  noisy  for  the  Instron  control  circuit. 
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The  load  cell,  LVDT's  and  the  clip  gauge  were  calibrated  prior  to  each  experiment. 
Readings  from  the  LVDT  and  the  load  cell  were  acquired  using  an  IBM  Personal 
Computer.  The  sampling  frequency  of  the  acquisition  unit  was  varied  according  to 
the  strain  rates  used  in  each  experiment. 

The  results  presented  in  this  paper  were  obtained  from  a  set  of  experiments  performed 
at  room  temperature.  The  nominal  temperature  in  the  test  laboratory  was 

approximately  19°C. 

It  is  known  that  the  stress  distribution  in  the  adhesive  layer  of  the  lap  shear  joint  is 
complex,  especially  around  the  comer  of  the  joint.  The  material  discontinuity  (i.e.  at 
the  adhesive  and  steel  adherend  interface)  results  in  a  region  of  high  peel  stress.  This 
region  of  high  stress,  although  highly  localised,  is  a  potential  failure  location.  To 
avoid  failure,  the  experiments  were  terminated  at  a  shear  strain  of-  0.4. 

2,2.3  Experimental  results 

To  determine  the  visco-plastic  response  of  the  adhesive  the  experiments  were 
conducted  at  a  range  of  strain  rates,  ~1.0e-2,  1.0e-3,  and  1.0e-4  and  temperatures,  viz: 
23C  (room  temperature)  and  60C.  A  set  of  initial  tests  was  also  performed  to  test 
repeatability  and  the  deviation  in  the  results  was  found  to  be  in  the  order  of  5%,  which 
can  be  attributed  to  various  sources  of  experimental  error. 

The  results  of  these  tests,  together  with  an  analytical  prediction  based  on  the 
analytical  formulation  outlined  in  section  3  are  presented  in  Figures' 2.3  and  2.4. 
These  diagrams  reveal  that  the  FM73  adhesive  exhibits  significant  visco-plasticity, 
which  in  these  tests  was  reflected  in  strain-rate  dependency,  even  at  room 
temperature,  (i.e.  well  below  its  glass  transition  temperature).  It  should  be  noted  that 
it  was  difficult  to  control  the  test  at  the  slowest  strain  rate.  This  was  the  result  of 
endeavouring  to  control  the  experiment  at  a  strain  rate  close  to  the  lower  limit  of  the 
machine,  in  terms  of  cross-head  speed.  The  yield  stresses  tabulated  in  Table  2.1  are 
calculated  from  the  experimental  data  using  Consideres'  constmction.  This  Table  also 
shows  a  consistent  value  of  the  shear  modulus  obtained  from  this  set  of  experiments. 

Since  the  yield  stress  of  the  adhesive  is  dependent  on  the  strain  rate,  it  follows  that  the 
adhesive  properties  must  have  a  strong  dependence  on  the  loading  history.  One 
method  of  revealing  this  time  dependency,  is  via  a  stress-relaxation  test  (i.e.  strain 
hold  test).  Figure  2.5  shows  the  result  of  loading  the  specimen  past  yield  at  a 
particular  strain-rate  and  then  holding  the  strain  at  a  constant  value.  This  figure 
reveals  a  very  large  stress  drop,  to  its  "threshold"  value,  as  a  result  of  the  strain  hold. 
The  level  of  the  stress  relaxation  seen  in  these  figures  reveals  that  the  "threshold"  to 
which  the  stress  relaxes  is  independent  of  the  value  of  the  plastic  strain,  see  [9,  13]. 
The  threshold  value  is  also  independent  of  the  strain  rate  and  determines  the  inelastic 
(irreversible),  and  thus  the  fatigue,  threshold  of  the  adhesive,  i.e.  the  level  beneath 
which  irreversible  effects  do  not  occur. 


Shear  Stress 

>  * 

Strain  rate 

— I - j - x - 1 - 

Experiment  Predicted 

50  - 

10-2 

10-3 

10-4 

* 

a  - 

*  ❖  o 

*  O  o  Q 

+  + 

-  +  4.  + 
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30  *  t 

a  □ 

a  o  ui 
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Figure  2.3  Stress  strain  curves  for  FM73  at  room  temperature. 
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Figure  2.4  Stress  strain  curves  for  FM73  at  60C. 
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The  stress  drop  during  stress-relaxation  to  its  "threshold"  value  results  in  a  conversion 
of  elastic  strain  into  inelastic  (irreversible)  strain,  with  the  total  strain  remaining 
constant.  This  results  in  recoverable  (i.e.  elastic)  strain  energy  being  converted  into 
irreversible  energy,  and  as  such  is  a  direct  reflection  of  the  degradation  of  the 
adhesive. 

•  This  degradation  can  be  avoided  by  designing  such  that  the  adhesive 
stresses/strains  are  beneath  the  "threshold"  value.  In  this  case  inelastic 
(irreversible)  strains,  and  hence  damage,  will  not  occur,  see  [10,1 1]. 

For  the  adhesive  under  consideration  this  "threshold"  stress  is  approximately  25  MPa, 
which  corresponds  to  a  "threshold"  energy  Wf  of  0.416  MPa.  As  can  be  seen  from 

Table  2. 1  the  "threshold"  stress  corresponds  to  the  lower  limit  of  the  yield  (shear) 
stress,  i.e.  the  value  of  the  yield  (shear)  stress  at  vanishingly  small  strain  rates.  This 
finding  is  intuitively  valid. 


Table  2.1 :  Yield  shear  stress  and  shear  modulus  as  a  function  of  strain  rate. 


Strain  Rate  (/sec) 

Yield  Stress  (MPa) 

Shear  Modulus  (MPa) 

8.68e-5 

26.78 

753 

8.78e.4 

32.81 

746 

8.67e.3 

35.94 

766 

8.80e-2 

39.06 

758 

2.2.4  Discussion 

The  section  has  shown  how  to  use  the  ASTM  thick  adherend  tests  to  measure  the 
structural  response  of  a  typical  thin  film  adhesive,  viz:  FM73,  at  both  room 
temperature  and  at  elevated  temperature. 

The  stress-relaxation  resulting  from  a  strain  hold  was  significant  and  reveals  the 
"threshold"  state,  below  which  irreversible  effects,  and  hence  fatigue  damage,  do  not 
occur.  This  "threshold"  value  is  independent  of  strain  rate  and  of  the  degree  of  prior 
plastic  deformation.  The  extensive  stress-relaxation,  which  occurs  in  a  matter  of 
seconds  or  minutes,  can  severely  degrade  the  performance  of  the  adhesive.  This 
degradation  is  due  to  the  conversion  of  recoverable  elastic  energy  to  irreversible 
inelastic  energy,  and  can  be  avoided  by  designing  such  that  the  design  state  is  below 
the  "threshold"  value. 

An  analytical  methodology  enabling  this  time-dependent  inelastic  deformation  of  the 
adhesive  to  be  included  in  the  design  process  is  given  in  Section  2.4. 
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2.4  Governing  Differential  Equations  For  Bonded  Joints/Repairs 

In  the  previous  section  we  have  outlined  a  simple  test  methodolgy,  based  on  the 
ASTM  thick  adherend  test  specimen,  for  measuring  the  structural  response  of  thin 
film  adhesives.  This  has  then  been  used  to  determine  the  structural  response  for  the 
adhesive  FM73,  which  is  commonly  used  in  composite  repair  technology.  From  the 
test  results  presented  we  have  seen  that  this  response  exhibits  significant  visco-plastic 
effects.  The  question  thus  arises: 

How  do  we  account  for  this  in  both  the  analysis  and  design  of  a  composite 
repair/joint?. 

To  this  end  let  us  consider  a  symmetric  double  lap  joint,  see  Figure  2.6,  subjected  to  a 
remote  load  P  where  the  lower  and  the  upper  adherend  thicknesses  and  Moduli  are  Tj, 

E|  and  T2,  E2  respectively  and  where  the  adhesive  has  a  thickness  t  and  an  elastic 
shear  modulus  of  Ga.  The  governing  differential  equations  relating  the  stresses  in  the 
adherends  <7  \  ,  <7  2  >  where  the  subscript  differentiates  between  the  lower  and  the 
upper  adherend,  to  the  adhesive  shear  stress  x,  in  the  joint  are: 


dax 

3c 


(11) 


&.+-L-0'  (12) 

&  Ti  ) 

If  the  adherends  are  elastic  then  Hooke's  law  applies  and  the  displacements  in  the 
lower  and  the  upper  adherents,  defined  as  uj  and  U2  respectively,  are  related  to  the 

stresses  through  the  Young's  moduli ;  viz: 


Ei 


3c 


a 


(13) 

(14) 


whilst  the  shear  strain  y  in  the  adhesive  is  related  to  the  relative  displacements  in  the 
adherends;  viz: 


y  =  (u!  -u2)/t 


(15) 


This  system  of  equations  results  in  the  following  relationship  between  the  adhesive 
shear  stresses  and  shear  strains;  viz: 


dx1 


=  rO-ltEiTi  +  lltEiTi) 


(16) 


where  from  equilibrium  considerations  the  load  P  applied  to  the  upper  adherend  must 
equal  the  integral  of  the  adhesive  shear  stresses;  viz: 
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p=  1  T  dx 

0 


(17). 


where  1  is  the  overlap  length  of  the  joint.  To  complete  this  system  of  equations  it  is 
necessary  to  prescribe  the  relationship  between  the  adhesive  shear  strain  y  and  the 
adhesive  shear  stress  t .  Since  this  relationship  is  usually  in-elastic  it  best  to  consider 
a  general  functional  form;  viz: 


(18) 


where  the  inelastic  shear  strain  rate  y  1  is  considered  to  be  an  arbitrary  function,  fix  , 
y  *)  of  the  shear  stress  t  and  the  shear  strain  rate  y  \ 

This  formulation  extends  the  approach,  used  by  Hart-Smith  [4],  to  allow  for  the 
visco-plastic  response  of  the  adhesive. 

The  solution  to  this  set  of  non-linear  time  dependent  equations  can  be 'quite  messy  and 
time  consuming.  To  over  come  this  we  use  Glinka’s  [14]  technique  for  estimating  the 
peak  stresses  and  strains  in  the  adhesive. 
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Figure  2.6:  A  schematic  diagram  of  a  symmetric  double  lap  joint 

2.4.1  Glinka’s  Approach  For  Estimating  The  Peak  Stresses  and  Strains 

Structural  components,  albeit  composite  repairs,  bonded  joints  or  stiffener  runnouts, 
are  frequently  subjected  to  complex  loading  spectra.  These  alternating  loads  tend  to 
initiate  fatigue  cracks  at  notches  and  at  other  regions  of  high  stresses.  Historically  the 
field  of  fatigue  has  been  classified  into  a  number  of  specific  areas;  viz:  high-cycle  and 
low-cycle  fatigue;  fatigue  of  notched  members;  the  initiation  and  propagation  of 
cracks  and  fatigue  life  extension  techniques.  Fatigue  initiation  and  crack  growth 
programs  require  an  accurate  knowledge  of  the  local  notch  tip  stresses  and  strains. 
These  quantities  can  be  determined  in  several  ways,  viz:  via  direct  strain  gauge 
measurements,  using  finite  element  analysis  or  by  using  approximate  methods,  such 
as  the  Glinka’s  approach  [14],  that  relate  local  stresses  and  strains  to  their  remote 
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values.  To  this  end  the  present  report  briefly  outlines  the  Glinka  approach  for 
calculating  the  peak  shear  stresses  and  shear  strains.  Attention  is  then  focused  on 
developing  a  simple  method,  which  combines  modem  constitutive  theory  with  the 
Glinka  approach  to  calculate  the  peak  (visco-plastic)  adhesive  stresses  and  strains. 

The  design  of  composite  repairs  requires  a  knowledge  of  the  interlaminar  shear 
stresses  and  strains.  These  quantities  are  necessary  to  design  against  interlaminar 
failure  in  the  composite  repair.  At  the  moment  these  stresses  can  only  be  determined 
using  finite  element  analysis. 

Let  us  first  define  the  peak  adhesive  shear  stresses  and  strains  obtained  via  an  elastic 
solution  as  x  and  y  (=x/G)  respectively.  Let  us  next  define  the  peak  adhesive  shear 
stresses  and  strains  obtained  via  an  in-elastic  solution  as  x  and  y  respectively.  In  this 
case,  Hooke's  law  cannot  be  used  to  relate  the  peak  shear  stress,  x,  to  the  peak  strain, 
y.  Instead  Glinka’s  rule  [2]  can  be  used  to  compute  the  peak  stresses  and  strains. 
According  to  this  hypothesis  the  peak  in  the  strain  energy  density  field  obtained  via  in 
an  elastic  plastic  analysis  is  the  same  as  the  peak  strain  energy  obtained  via  a  purely 
elastic  analysis;  viz: 


1/2 1 1  =  1 2  /2G  =  Jajj  dsy  =  lx  dy  (19) 

It  should  be  noted  that  this  relationship  also  follows  from  the  Hart-Smith  solution  [4], 
As  such  the  Hart-Smith  solution  is  contained  as  a  special  case  of  Glinka’s  hypothesis. 
However,  in  [4]  visco-plastic  effects  were  not  considered. 

The  solution  process  thus  involves  first  solving  for  the  peak  elastic  stresses  and  strains 
in  the  repair.  Once  a  valid  stress  strain  relationship  for  the  epoxy  is  known  we  then 
use  equation  (19)  together  with  the  elastic  solution  to  determine  the  peak  (in-elastic) 
adhesive  stresses  and  strains.  A  simple  fortran  based  computer  code  has  been 
developed  to  do  this.  This  program  was  given  to  the  USAF  at  the  June  99  meeting. 

In  the  past  few  years  there  has  been  an  increasing  interest  in  the  use  of  unified 
constitutive  models  for  predicting  the  inelastic  response  of  structural  materials,  i.e. 
visco-plasticity,  creep,  stress  relaxation  etc.  These  models  overcome  many  of  the 
deficiencies  of  the  classical  approaches  to  the  inelastic  behaviour  of  materials.  The 
present  work  uses  the  formulation  presented  in  [9,  13]  to  represent  the  visco-plastic 
response  of  thin  film  adhesives. 


This  formulation  uses  two  internal  state  variables,  back  stress  Qtj  (deviatoric)  tensor 
and  drag  stress  Z  to  define  the  state  of  the  material.  Here  the  inelastic  deformation  is 
driven  by  the  over  stress  which  is  defined  as  (SyQy).  Similar  types  of  theories  have 
been  widely  used  to  describe  the  behaviour  of  thermoplastics. 

In  this  approach,  the  flow  equation  is  defined  as: 


s\j  =  D  exp 


-A 


(  72 


V3  K2j 


('Sy  Qg) 

' 


(20) 


where  is  the  inelastic  strain  rate  tensor,  Stj  is  the  deviatoric  stress  and  K2  = 
}r(Sy -Clij)(Sy -Qy)  .  Here  D  A  and  n  are  material  constants.  The  term 
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(Sy-Cly)/  in  equation  (20)  is  the  normalised  overstress  and  defines  the 

straining  direction.  The  drag  stress  is  history  dependent  and  evolves  with  the  effective 
inelastic  strain  s'  ■  The  initial  drag  stress  is  Zq,  final  Z\,  and  mg  controls  the  rate  of 
evolution.  The  “growth”  of  the  drag  stress  is  defined  by  equation  (21). 

Z  =  Zl+(Z0+Z1)e~mo£'  (21) 

The  back  stress  is  initially  zero  and  its  evolution  is  controlled  by  an  evolution 
equation,  viz: 


£W2V/i4-  <22) 

^  ^max 

where  e‘  =  is  the  effective  inelastic  strain  rate  .  Here,  fj,  f2  and  Qmax  are 

material  constants.  The  effective  inelastic  strain  rate  is  defined  in  the  normal  fashion. 


2.4.2  Validation 

Consider  the  symmetric  double  lap  joint  shown  in  Figure  2.6.  This  joint  has  an 
overlap  length  of  90  mm,  and  consists  of  two  identical  aluminium  adherends  with  an 
elastic  modulus  E  =  70000  MPa,  Mid  t=1.5  mm.  The  aluminium  adherends  and 
doublers  were  assumed  to  be  bonded  together  using  FM73.  The  adhesive,  FM73  film, 
was  0.2  mm  thick  and  had  a  shear  modulus  G  of  750  MPa.  Its  visco-plastic  behaviour 
is  described  using  the  state  variable  formulation  described  above.  The  material 
parameters  for  the  constitutive  law  for  FM73  at  room  temperature  are  given  in  Table 
2.2. 

This  joint  was  analysed  using  both  the  Glinka  approach,  described  above,  and  the 
ABAQUS  analysis  program  together  with  the  associated  UMAT  material  subroutine. 
In  each  case  a  remote  stress  of  401  MPa  was  assumed  to  be  monotonically  applied  to 
the  ends  of  the  inner  adherend  with  the  load  increasing  from  0  to  its  maximum  load  in 
either  0.1  sec,  1  sec,  10  sec  or  100  sec.  A  comparison  of  the  calculated  peak  shear 
strains  y  in  the  joint  is  shown  in  Table  2.3  . 

2.4.3  Application  to  the  analysis  of  double  lap  joints 

Having  thus  validated  the  Glinka  approximation  a  further  series  of  finite  element 
analysis  of  various  double  lap  joints  (Figure  2.7)  were  then  performed.  The  joints 
analysed  in  this  paper  had  (L/2t)  ratios  of  7.5  and  15  and  a  range  of  load  histories 
were  evaluated,  see  (Figure  2.7).  These  joints  had  overlap  lengths  of  90  mm  (for 
L/2t=15),  and  45  mm  (for  L/2t=7.5)  respectively.  The  elastic  modulus  of  the 
adherends  was  E=70000  MPa.  The  thicknesses  were  3  mm  (for  the  inner  adherend) 
and  1.5  mm  (for  the  doubler).  The  aluminium  adherends  and  doublers  were  assumed 
to  be  bonded  together  using  FM73.  The  FM73  adhesive  was  0.2  mm  thick  and  had  a 
shear  modulus  of  750  MPa.  These  joints  were  loaded  to  a  maximum  remote  stress  of 
401  MPa.  with  the  load  going  from  0  load  to  its  maximum  value  in  either  0.1  sec,  1 
sec,  10  sec  or  100  sec.  The  stress/strain  relationship  at  the  worst  (critical)  point  in  the 
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adhesive  layer  during  these  load  up  sequences  is  shown  in  Figures  2.8  (L/2t=15)  and  4 
(L/2t=7.5). 


Table  2.2:  Typical  material  parameters  for  FM73  film  at  room  temperature,  see  [3] 


A 

D 

fl 

(MPa) 

h 

n 

m0 

^max 

zo 

(MPa) 

Zl 

(MPa) 

Room 

temp. 

800 

0.2 

0.85 

10 

12.5 

120 

350  j 

These  figures  illustrate  that  the  stress/strain  behaviour  of  the  adhesive  exhibits 
significant  visco-plasticity  and,  as  a  results,  is  dependent  on  the  loading  rate  (the 
higher  the  loading  rate,  the  higher  is  the  apparent  yield  stress).  The  maximum  stress 
level  (apparent  yield  stress)  in  the  adhesive  is  higher  when  the  joint  is  loaded  more 
rapidly.  These  figures  also  show  that  the  strain  level  in  the  adhesive  is  dependent  on 
the  loading  history.  Here  the  higher  loading  rate  resulted  in  lower  adhesive  strain 
levels. 


Table  2:3  Comparison  of  peak  shear  strain  y  with  ID  Glinka  formulation,  L/2t  =  15 


Dimension 

Loading  Time  (sec) 

0.1 

1 

10 

100 

ID  algorithm 

0.143 

0.152 

0.163 

0.177 

Finite  element  2D/3D 

0.152 

0.165 

0.179 

'  0.189 

%  Diff 

5.9 

7.9 

6.3 

In  contrast  to  the  stress  or  strain  based  approaches  Hart-Smith’s  [4]  analysis  of  a 
double  lap  joint  revealed  that  the  load  carrying  capacity  of  the  bond  can  be  estimated 
from  the  energy  density  at  the  worst  point  in  the  adhesive.  This  design  methodology 
was  based  on  the  assumption  that  the  adhesive  behaviour  can  be  approximated  using  a 
single  bi-linear  stress/strain  curve.  This  simplification  is  not  necessary.  With  the 
extended  analysis  given  in  this  report  the  adhesive  can  be  allowed  to  undergo  both 
creep  and  visco-plastic  behaviour.  As  shown  in  Figures  2.8  and  2.9  the  stress/strain 
relationship  at  any  location  along  the  adhesive  layer  will  be  governed  by  the  local 
adhesive  strain  rate  (i.e.  the  higher  the  local  adhesive  strain  rate,  the  larger  will  be  the 
apparent  yield  stress  of  the  adhesive). 


To  determine  how  the  visco-plastic  behaviour  of  the  adhesive  affects  the  load  carrying 
capacity  of  the  joint  the  energy  density  at  the  worst  point  in  the  adhesive  layer  was 
also  calculated.  The  results  are  shown  in  Table  2.4.  The  first  observation,  apparent 
from  this  Table,  is  that  the  energy  density  at  the  worst  point  is  slightly  higher  when 
the  (L/2t  ratio)  is  7.5.  Hence  one  would  expect  the  joint  with  a  lower  (L/2t)  ratio  to 
have  a  slightly  lower  load  carrying  capacity.  More  importantly  we  find  that,  whilst  the 
local  stress/strain  behaviour  of  the  adhesive  are  largely  dependent  on  the  loading  rate, 
the  energy  density  at  the  worst  point  in  the  adhesive  is  relatively  independent  of  the 
loading  rate.  There  is  only  approximately  8%  variation  (refer  to  Table  2.4)  in  the 
energy  density  from  the  slowest  to  the  fastest  loading  rate. 
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The  worst  point  in  the  joint 


Figure  2.7:  A  quarter  of  a  symmetric  double  lap  joint 

This  work  thus  confirms  that  the  design  hypothesis  first  presented  by  the  Hart-Smith 
[4]  for  estimating  the  load  carrying  capacity  of  a  double  lap  joint  can  be  used  for  thin 
film  adhesives,  such  as  FM73  and  FM300,  which  exhibit  significant  visco-plastic 
effects  provided  that  failure  occurs  in  the  adhesive. 

The  Hart-Smith  approach  does  not  cover  interlaminar  failure  in  the  composite  patch, 
see  Section  2.5.  In  this  case,  since  the  shear  stresses  are  continuous  across  the  patch 
adhesive  interface  it  is  necessary  to  determine  the  stress  state  in  the  adhesive  which, 
as  we  have  seen,  requires  an  allowance  for  the  visco-plastic  response  of  the  adhesive. 


Table  2.4:  Effect  of  loading  rate  on  energy  density  both  L/2t=7.5  and  15  at  room  temp 


.Energy  Density 
(MPa) 

L/2t 

Loading  Time 
(sec) 

0.1 

1 

10 

% 

Diff 

7.5 

Adhesive  layer 

42.6 

44.4 

15 

Adhesive  layer 

41.6 

43.7 

2.4.4.  Remarks 

This  analysis  reveals  that,  at  room  temperature,  when  the  loading  history  is  simply 
monotonic  the  visco-plastic  nature  of  the  adhesive  in  the  double  lap  joint  has  limited 
effects  on  the  strain  energy  density,  and  thus  on  the  load  carrying  capacity  of  the 
adhesive  bond  (repair),  see  Figures  2.8  and  2.9. 

In  contrast  the  peak  adhesive  shear  stresses  and  strains  are  strongly  effected. 
Consequently,  since  the  shear  stresses  are  continuous  across  the  composite  patch- 
adhesive  interface,  for  an  accurate  knowledge  of  the  interlaminar  stresses  in  the 
composite  repair  it  is  important  that  these  visco-plastic  effects  to  be  incorporated  in 
the  design.  This  can  be  achieved  by  using  Glinka’s  approach.  To  this  end  a  simple 
fortran  program  has  been  developed  thereby  negating  the  need  for  a  complex  non¬ 
linear  finite  element  analysis.  This  formulation  has  been  validated  and  was  given  to 
the  USAF  at  the  June  meeting.  In  this  fashion  we  are  now  able  to  assess  for  failure  in 
the  adhesive,  the  metal  skin,  fibre  failure,  or  interlaminar  failure  of  the  composite 
patch. 
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2.5  Interlaminar  failure 

In  the  previous  sections  we  have  concentrated  on  failure  due  to  failure  of  the  adhesive 
bond.  However,  when  selecting  the  damage/degradation  to  be  considered  in  the  repair 
design/substantiation,  the  following  failure  modes,  unique  to  composite  materials, 
should  also  be  included;  viz: 

(a)  Failure  of  doubler/repair  by  fibre  failure, 

(b)  Failure  in  the  adhesive  (cohesion  failure), 

(c)  Failure  at  the  adhesive/substrate  interface(s)  (adhesion  failure),  and 

(d)  Failure  of  the  doubler/repair  by  interlaminar  failure. 

a  more  comprehensive  review  of  the  potential  failure  mechanisms  is  given  in  [15]. 

For  repairs  to  primary  structures  interlaminar  failure  is  particularly  common  as  can  be 
seen  in  Figure  2.10  which  shows  interlaminar  failure  of  a  composite  doubler  on  an 
Fill  aircraft,  see  [16]  for  details  of  this  repair. 
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Figure  2.9:  Adhesive  shear  stress-strain  curves  at  a  critical  point  in  the  joint,  L/2t=7.5. 

This  interlaminar  failure  mechanism  is  also  frequently  seen  in  laboratory  tests,  see 
Figure  2.11  [15,  17,  18],  on  boron  epoxy  repairs  to  surface  flaws.  This  test  program, 
which  is  documented  in  [  17],  dealt  with  the  repair  of  surfaces  cracks  40  mm  long  and 
5.7  mm  deep  in  an  11.2  mm  x  106.8  x  304.8  aluminium  alloy  specimen  which  was 
subjected  to  constant  amplitude  loading  with  a  peak  stress  of  69.4  MPa.  with  a  stress 
ratio  R  of  0.01.  In  these  tests  the  failure  mechanism  was  generally  similar  to  that  seen 
in  the  F-111C  repair  program  in  that  failure  was  in  the  boron  epoxy  repair.  This  is 
clearly  shown  in  Figure  2.10  where  boron  fibres  can  (again)  be  seen  attached  to  the 
surface  of  the  aluminium.  This  failure  mechanism  is  also  evident  in  the  test  program 
reported  in  section  4.4. 

It  is  thus  clear  that  the  major  design  considerations  that  need  to  be  addressed  for  a 
composite  repair  are,  viz: 

1.  The  maximum  stress  intensity  factor  should  be  as  low  as  possible  and 
preferably  below  the  critical  value  for  fatigue  crack  growth  Kth  in  the 
material. 

,2.  The  maximum  adhesive  design  state  should  be  below  the  value  at  which 
fatigue  damage  accumulates  in  the  adhesive.  To  minimise  errors  in 
measuring  and  computing  the  adhesive  stresses  and  allowables  it  is  best  to 
follow  the  Hart-Smith  methodology  and  use  the  energy  in  the  adhesive  W, 
=1/2  a,j  s  jj  =  Jay  ds  jj,  as  the  design  variable.  These  critical  design  state 
measurements  can  be  obtained  using  the  ASTM  thick  adherend  test,  ASTM 
D  1002,  see  above.  This  approach,  which  was  originally  proposed  in  [9],  has 
been  validated  in  [  10, 11]. 
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3.  The  rate  of  delamination  growth  in  the  adhesive  can  be  determined  from  a 
knowledge  of  the  peak  ahesive  shear  stresses  and  strains.  For  FM73  one  such 
formulae  relating  the  delamination  length  b  to  the  number  of  cycles 
experienced  N  is  given  in  [1 8]  as 

db/dN  =  5.43  105  (Ay ) 1096 

4.  The  composite  patch  must  not  experience  failure  by  interply  delamination. 
This  can  be  checked  by  ensuring  that  the  polynomial  failure  criteria  is  not 
greater  than  one.  The  commonly  used  failure  criteria  are:  Tsai-Hill,  Hoffman 
and  Tsai-Wu.  These  failure  criteria  are  generally  written  in  the  form: 

Tsai-Hill  criterion:  In  this  case  failure  is  assumed  to  occur  when 


(23) 


Here  the  material  is  assumed  to  have  equal  strengths  in  tension  and 
compression,  i.e.  Xt  =  Xc=  X.andY,  =  Yc  =  Y. 


Hoffman  criterion:  Failure  is  assumed  to  occur  when 


1  1  ")  1  1  ^  of  g\  g\2  g}g2 

xt  X  J  1  U  Yj  XtXc  Y,YC  S2  XtXc 


(24) 


Tsai-Wu  criterion:  Failure  is  assumed  to  occur  when 


(25) 


The  coefficient  Fn  is  experimentally  determined  from  test  specimens 
under  bi-axial  loading  and  Fn  must  satisfy  a  stability  criterion  of  the 
form 


creates  some  complication  in  the  use  of  this  theory.  It  has  been 
suggested  that  Fn  be  set  to  zero. 


(26) 


Disbonded  area  has  been 
cut  away. 

Note  fibres  adhering  to  the 
surface  of  the  wing  pivot 
fitting. 


Figure  2.10.  Close  of  delaminated  doubler  on  Fill  aircraft,  looking  from  above  the 
repair,  showing  boron  fibres  still  attached  to  the  wing  pivot  fitting. 


-350  mm 


510mm 


3  mm  aluminium  honeycombe 


Bonded  boron  patch 


Aluminium  plate  11.2  mm  thick 


Aluminium  plate  11.2  mm  thick 


Bonded  boron  patch 

Figure  2.1 1  Geometry  of  the  surface  flaw  specimen 


Figure  2.12  Close  up  of  failure  surface  showing  the  1st  layer  of  the  repair  still  attached 
to  the  specimen 

The  symbols  used  in  Equations  (23)  to  (26)  are  defined  as: 

Allowable  tensile  stress  in  the  principal  x  (or  Indirection  of  the 
material 

Xc  Allowable  compressive  stress  in  the  principal  x  (or  Indirection 
of  the  material 

Yt  Allowable  tensile  stress  in  the  principal  y  (or  2)-direction  of  the 

material 

Yc  Allowable  compressive  stress  in  the  principal  y  (or  2)-direction 
of  the  material 

S  Allowable  shear  stress  in  the  principal  material  system 

The  interlaminar  stresses  can  be  estimated  using  the  Glinka  formulation  outlined 
earlier  in  this  report. 


For  relatively  thin  repairs  the  average  stress  through  the  thickness  of  the  boron  patch 
should  not  exceed  1000  MPa.  This  failure  rule  is  a  modification  of  the  average  stress 
failure  criterion,  which  is  commonly  used  for  composite  materials. 

It  must  be  stressed  that  for  repairs  to  primary  structures  a  3-D  finite  element  analysis 
should  be  performed.  To  ensure  continued  airworthiness  this  analysis  should  include  a 
damage  tolerant  assessment  of  both  the  structure  and  the  composite  repair. 
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At  this  stage  it  should  be  noted  that  the  majority  of  commercially  available  finite 
element  programs,  e.g.  ABAQUS,  MSC-Nastran,  NEnastran,  MARC,  ANSYS  etc, 
use  iso-parametric  elements.  Furthermore,  when  analysing  three  dimensional 
problems  involving  composites  the  use  of  elements  with  mid-side  nodes  is  essential, 
i.e.  20  noded  hex  (brick)  elements,  15  noded  wedge  elements,  etc.  Unfortunately, 
ABAQUS,  MSC-Nastran,  MARC,  ANSYS  still  use  formulations  based  on  that  given 
in  Zienkiewicz  [23],  which  means  that  for  accurate  results  the  aspect  ratio’s  are 
limited  to  less  than  —5:1,  even  if  reduced  integration,  i.e.  the  use  of  2  x  2  x  2  Gauss 
points  in  formulating  the  element  stiffness  matrix,  is  used.  In  contrast  NEnastran  use 
an  enriched  20  noded  element,  where  the  basic  formulation  is  extended  to  include 
additional  shape  functions.  This  means  that  the  NEnastran  3D  (mid-side  nodded) 
elements  can  have  dramatically  greater  aspect  ratio’s.  Consequently,  when  analysing 
composite  repairs,  and  consequently  performs  better  when  analysing  composite 
repairs.  Equally  important  is  the  fact  that,  unlike  MSC  NASTRAN,  the  NEnastran  20 
noded  bricks  can  be  used  when  analysing  repairs  in  which  the  adhesive  is  behaving 
plastically  and  when  there  may  be  non-linear  geometric  effects.  This  feature  is 
particularly  important  when  determining  the  residual  strength  of  a  repair  with 
delamination  damage  in  the  adhesive. 


3.  ANALYSIS  OF  REPAIRS  TO  RIB  STIFFENED  PANELS 

3.1  Large  Panel  Fatigue  Test  Specimens 

After  consultation  with  the  USAF  it  was  decided  to  use  narrow  single  stiffener  panels. 
These  panels  are  currently  being  fabricated  for  shipment  to  Australia.  Testing  is 
awaiting  the  fabrication  and  delivery  of  these  specimens. 

3.2  Finite  Element  Analysis 

A  finite  element  study  for  the  composite  repair  of  cracked  rib  stiffened  panels  has 
now  been  completed.  In  this  study  the  stiffeners  were  assumed  to  be  riveted  to  the 
skin.  Symmetry  considerations  enabled  only  a  %  of  the  structure  to  be  model.  Several 
typical  examples  are  shown  in  Figures  3. 1-3.5.  In  this  study  the  crack  length,  patch 
thickness,  skin  thickness,  stiffener  width,  stiffener  spacing,  and  stiffener  depth  were 
varied  to  determine  their  effects  on  the  stress  intensity  factors. 


Figure  3.1  Schematic  of  three  dimensional  %  finite  element  model  of  stiffened  panel. 
Note  the  definitions  with  respect  to  stiffener  dimensions. 


Figure  3.2.  Plan  view  of  a  finite  element  model  shown  in  Figure  3.1  showing  mesh 
spacing  in  and  around  the  crack  tip.  8mm  crack  length,  4mm  skin,  1mm  patch 
thickness,  5.5mm  stiffener  width,  75mm  stiffener  spacing  and  12  mm  stiffener  depth 


Figure  3.3  Schematic  of  a  finite  element  model  with  20mm  crack  length,  3.16mm 
skin,  1.4mm  patch  thickness,  3.5mm  stiffener  width,  127.5mm  stiffener  spacing  and 
37.5  mm  stiffener  depth 


Figure  3.5  Schematic  of  another  finite  element  model  with  40mm  crack  length, 
2.5mm  skin,  1mm  patch  thickness,  3.5mm  stiffener  width,  150mm  stiffener  spacing 
and  12mm  stiffener  depth. 

To  create  the  required  solution  space  this  process  was  automated  and  a  large  number 
of  geometrical  configurations,  ~  1500,  were  evaluated.  Each  run  took  ~  6  minutes  of 
cpu  time  on  an  SGI  R10000  workstation..  ' 

The  resultant  solutions  are  summarised  below.  It  is  hoped  that  the  differences  (if  any) 
in  the  stress  intensity  factors  caused  by  the  change  in  geometry  may  be  explained  by 
simple  analytical  relationships.  The  results  from  the  finite  element  analysis  will 
subsequently  be  used  to  predict  the  behaviour  of  composite  repairs  to  cracks  in  rib 
stiffened  panels. 

3.3  Results  for  stress-intensity  calculations  from  the  finite  element  analysis  of 
a  boron  epoxy  composite  repair  of  a  crack  over  a  stiffening  element 

In  this  study  the  adhesive  thickness  (ta)  was  kept  constant  at  0.2mm  and  the  following 
labelling  convention  was  adopted  : 

“File  name”  =  Skin-thickness.. Stiffener-width..Patch-thickness.. Stiffener¬ 
spacing.  .  Stiffener-depth 

Thus  a  designation  with  3. 16. . 3. 5. .1.. 127.5. .12  describes  a  stiffened  skin  with  : 

3.16  mm  skin  thickness, 

3.5  mm  stiffener  width, 

1  mm  patch  thickness 

127.5  mm  stiffener  spacing 
12mm  stiffener  depth 

For  the  analyses,  the  following  conditions  were  considered  : 
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A  centre  crack  in  a  2.5,  3.16  or  4mm  thick  aluminium  stiffened  panel  with  E=  72,000 
MPa,  and  v  =  0.33.  Three  patch  thicknesses  were  considered,  viz:  1, 1.4  and  1.85  mm 
thick  and  the  modulii  of  the  patch  were  taken  to  be  En  =208,000  MPa,  E22  =  25432 
MPa,  V12  =  0.183,  G12  =  G13  =  G23  =  7241  MPa. 


The  panel  had  dimension  of  375  mm  (length)  x  400  mm  (width)  and  the  plan  form  of 
the  composite  repair  was  275  mm  (length)  x  400  mm  (width).  The  adhesive  was  taken 
to  be  0.20  mm  thick  with  a  shear  modulus  of  375  MPa  and  v  =  0.33.  The  analysis 
used  20  noded  isoparametric  3D  bricks  with  single  layer  of  elements  through  the  skin, 
adhesive  and  patch.  15  noded  isoparametric  3D  wedge  shaped  bricks  were  used  at  the 
crack  tip  location  through  the  entire  thickness  of  the  specimen.  The  skin  in  the  panel 
was  subjected  to  a  remote  uniform  stress  of  229.8MPa. 

The  solution  was  well  conditioned  and  the  accuracy  of  the  solution  was  evaluated  by 
performing  two  separate  analyses.  One  used  optimum  integration  to  form  the  stiffness 
matrix,  i.e.  2  x  2x  3  Gauss  quadrature  points,  the  second  used  full  integration.  Both 
analyses  gave  results  within  less  than  1%  of  each  other. 

To  illustrate  the  asymptotic  effects  of  the  solution  the  results  for  repairs  to  centre 
notch  cracks  will  now  be  presented.  To  this  end  various  plots  of  stress  intensity 
factors  vs.  crack  length  are  presented.  They  are  used  to  illustrate  the  effects  of  varying 
the  various  individual  geometrical  parameters  on  Ki. 

The  values  presented  in  these  analyses  are  taken  from  the  bottom  surface  of  the  skin 
(just  above  the  stringer).  As  expected,  these  gave  consistently  higher  values  of  Ki 
when  compared  to  the  stress  intensity  values  at  the  skin-adhesive  interface. 

3.3.1  Comparison  of  results  for  several  different  stiffener  spacing  and  two 
different  skin  thicknesses. 

This  study  focused  on  determining  the  effect  of  stiffener  spacing  on  the  stress 
intensity  factor.  The  geometrical  parameter  considered  are  shown  in  Table  3.1.  The 
results  are  summarised  in  Figures  3.6  and  3.7  where  in  each  case  the  stiffener  width 
and  depth  was  3.5mm  and  12mm  respectively.  From  these  figures  we  see  that  the 
stress  intensity  factor  Ki  is  a  weak  function  of  the  stiffener  spacing. 


Table  3.1  Geometrical  parameters  considered. 


Skin  thickness 
(mm) 

Stiffener  spacing 
(mm) 

Patch  thickness 

2.5 

75 

1.0 

3.6 

127.5 

1.8 

150 
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Figure  3.6  :  Plot  of  Ki  vs.  a  for  2.5mm  stiffened  skin  with  three  different  stiffener 
spacings.  Data  files  analysed  were  :  2. 5. . 3. 5. .1.. 75. .12,  2.5. . 3. 5. .1.. 127. 5. .12,  and 
2.5. . 3.5. .1.. 150.. 12. 


3.16mm  skin 


Figure  3.7  :  Plot  of  Ki  vs.  a  for  3.16mm  stiffened  skin  with  three  different  stiffener 
spacings.  Data  files  analysed  were:  3. 16. . 3. 5. .1..' 75. .12,  3. 16. . 3. 5. .1.. 127. 5. .12,  and 
3.16..3.5..1..150..12 

3.3.2  Comparison  of  results  for  several  stiffener  spacings  and  stiffener  depths 
and  two  different  skin  thicknesses. 

This  study  focused  on  confirming  the  effect  of  stiffener  spacing  on  the  stress  intensity 
factor.  In  this  study  the  stiffener  depths  were  varied,  see  Table  3.2,  and  the  results  are 
summarised  in  Figures  3.8  and  3.9  and  we  again  see  that  the  stress  intensity  factor  Ki 
is  a  weak  function  of  the  stiffener  spacing. 
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Table  3.2  Geometrical  parameters  considered. 


Skin  thickness  (mm) 

Stiffener  spacing 
(mm) 

Patch  thickness 

Stiffener  depth 
(mm) 

2.5 

75 

1.0 

12 

3.16 

127.5 

50 

150 

Figure  3.8  :  Plot  of  Ki  vs.  a  for  2.5mm  stiffened  skin  with  three  different  stiffener 
spacings  and  two  different  stiffener  depths.  Data  files  were  :  2. 5. . 3. 5. .1.. 75. .12, 

2.5.. 3  5..1..127.5..12,  2.5..3.5..1..150..12,  2.5..3.5..1..75..50,  2.5..3.5..1..127.5..50, 

2.5.. 3.5..1..150..50. 

3.3.3  Comparison  of  results  for  various  stiffener  spacing  and  patch  thickness 
for  two  different  skin  thicknesses. 

This  study  focused  on  confirming  the  effect  of  stiffener  spacing  on  the  stress  intensity 
factor.  In  this  study  the  patch  thickness  was  varied,  see  Table  3.3,  and  the  results  are 
summarised  in  Figures  3.10  and  3.11  and  we  again  see  that  the  stress  intensity  factor 
Ki  is  a  weak  function  of  the  stiffener  spacing. 

Figures  3.6  to  3.11  show  that  the  stress  intensity  factors  Ki  exhibit  an  asymptotic 
behaviour  as  the  crack  length  increases,  see  Figure  3.18,  although  as  the  skin 
thickness  increases  the  crack  length  required  to  achieve  this  asymptote,  K«  also 
increases.  Indeed,  for  the  2.5  mm  sheet  thickness  with  a  1.  mm  patch  the  limiting 
value  of  Km  has  essentially  been  reached  at  a  crack  length  of  between  20  -  30mm,  see 
Figure  3.10. 
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Figure  3.9:  Plot  of  K]  vs.  a  for  3.16mm  stiffened  skin  with  three  different  stiffener 
spacings  and  two  different  stiffener  depths.  Data  files  were  :  2. 5. . 3. 5. .1.. 75. .12, 

2.5. .  3.5. .1.. 127.5. .12,  2.5..3.5..1..150..12,  2.5..3.5..1..75..50,  2.5..3.5..1..127.5..50, 

2.5.. 3.5..1..150..50. 


Figure  3.10:  Plot  of  Ki  vs.  a  for  2.5mm  stiffened  skin  with  three  different  stiffener 
spacings  and  two  different  patch  thicknesses.  Data  files  were:  2. 5. . 3. 5. .1.. 75. .12, 

2.5.. 3.5..1.. 127.5. .12,  2.5..3.5..1..150..12,  2.5..3.5..1.8..75..12,  2.5..3.5..1.8..127.5..12, 

2.5. . 3.5. .1.8. .150.. 12 

In  contrast  for  the  3.16  mm  sheet  thickness  with  a  1.  mm  patch  the  limiting  value  has 
not  been  reached  by  40  mm,  see  Figure  3.11.  Similarly  the  crack  length  required  to 
reach  K«,  is  also  a  function  of  the  patch  thickness.  For  example  in  the  case  considered 
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above  when  the  3.16  mm  sheet  thickness  has  a  1.  8  mm  thick  patch  this  limiting  value 
is  essentially  achieved  with  a  crack  length  of  (say)  20  mm,  see  Figure  3.11. 

This  asymptotic  process  is  also  a  function  of  the  depth  of  the  stiffener.  For  example 
for  the  2.5  mm  skin  by  the  time  the  crack  length  is  40mm  the  limiting  value  of  K«,  has 
been  reached  for  stiffener  depths  of  both  12  and  50mm,  see  Figure  3.8.  In  contrast  for 
the  3.16  mm  thick  skin  with  a  crack  length  of  40  mm  the  solution  associated  with  a 
50  mm  deep  stiffener  has  reached  its  asymptote  whilst  that  for  a  12  mm  deep  stiffener 
has  not,  see  Figure  3.9. 


Table  3.3  Geometrical  parameters  considered. 


Skin  thickness  (mm) 

Stiffener  spacing 
(mm) 

Patch  thickness 

Stiffener  depth 
(mm) 

2.5 

75 

1.0 

12 

3.16 

127.5 

1.8 

150 

3.16mm  skin 


50.00 

40.00 

%  30.00 
» 

™  20.00 

2 

*  10.00 
0.00 


75mm  spacing,  1mm 
patch 


127.5mm  spacing, 
1mm  patch 


...  150mm  spacing, 
1mm  patch 

— 75mm  spacing, 
1.8mm  patch 

■3K — 127.5mm  spacing, 
1.8mm  patch 


— 150mm  spacing, 
1.8mm  patch 


Figure  3.11:  Plot  of  Ki  vs.  a  for  3 . 1 6mm  stiffened  skin  with  three  different  stiffener 
spacings  and  two  different  patch  thicknesses.  Data  files  were:  3. 16. . 3. 5. .1.. 75. .12, 

3. 1 6. .3. 5..  1..  127.5..  12  3. 1 6.. 3. 5. .1..  150..  12,  3.16..3.5..1.8..75..12, 

3. 16. .3. 5..  1.8..  127. 5. .12,  3. 16. .3. 5. .1.8. .150..  12 


3.3.4  Comparison  of  results  for  various  stiffener  widths  and  depths  for  two 
different  skin  thicknesses. 

This  study  focused  on  determining  whether  the  stress  intensity  factor  Ki  was  more 
strongly  effected  by  the  stiffener  depth  than  the  stiffener  width.  The  parameters 
considered  are  shown  in  Table  3.4  and  the  results  of  this  study  are  shown  in  Figures 
3.12-3.14.  This  study  found  that  K]  was  more  strongly  influenced  by  the  depth  than 
by  the  width  of  the  stiffener.  Indeed,  as  can  be  seen  from  Figures  3.13  and  3.14  the 
stiffener  width  had  little  effect  on  stress  intensity  factor(s).  It  thus  appears  that  the 
stress  intensity  factor  Ki  is  more  strongly  influenced  by  the  moment  of  inertia  I  than 
by  the  cross  sectional  area  of  the  stiffener. 
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Table  3.4  Geometrical  parameters  considered. 


2.5mm  skin 


Figure  3.12:  Plot  of  K  vs.  a  for  2.5mm  stiffened  skin  with  two  different  stiffener 
widths  (w)  and  four  different  stiffener  depths  (d).  Data  files  were:  2.5. . 3. 5. .1.. 75. .12, 

2. 5. .  3.5. .1.. 75. .23.24,  2.5..3.5..1..75..37.5,  2.5..3.5..1..75..50,  2.5..5.5..1..75..12, 

2. 5. .  5.5. .1.. 75. .23 .24,  2.5..5.5..1..75..37.5,  2.5..5.5..1..75..50 


Figure  3.13:  Plot  of  K  vs.  a  for  3.16mm  stiffened  skin  with  two  different  stiffener 
widths  (w)  and  four  different  stiffener  depths  (d).  Data  files  were:  3. 16. . 3. 5. .1.. 75. .12, 

3.16.. 3.5..1..75..23.24,  3.16..3.5..1..75..37.5,  3.16..3.5..1..75..50  3.16..5.5..1..75..12, 

3. 16. .  5.5. .1.. 75. .23. 24,  3.16..5.5..1..75..37.5,  3.16..5.5..1..75..50 
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Figure  3.14:  Plot  of  K  v’s.  a  for  a  2.5mm  stiffened  skin  with  three  different 
stiffener  widths.  Data  files  were  2. 5. . 3. 5. .1.. 75. .12,  2. 5. . 4.6. .1.. 75. .12, 
2.5. . 5.5. .1.. 75. .12 


Figure  3.15:  Plot  of  K  vs.  a  for  3.16mm  stiffened  skin  with  three  different  stiffener 
widths  Data  files  were  3.16..3.5..1..75..12,  3.16..4.6..1..75..12,  3.16..5.5..1..75..12 

3.3.5  Relationship  between  the  ratio  of  the  patch  thickness  to  skin  thickness 
and  the  stress  intensity  factor  Kj  for  5  different  crack  lengths. 

This  study  has  focused  on  determining  the  effect  of  the  ratio  of  patch  thickness  to  skin 
thickness  (tp/ts).on  the  stress  intensity  factor  Ki.  To  this  end  nine  different  geometrical 
configurations  with  three  different  skin  thicknesses  were  considered,  see  Table  3.5, 
and  results  of  this  study  are  shown  in  Figure  3.16. 
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This  study  found  that  to  a  first  approximationKi  was  a  linear  function  of  the  ratio  of 
patch  thickness  to  skin  thickness  (tp/ts). 


Table  3.5  Geometry  configurations  studied 


Panel/model  geometry 
descriptors 

tp/ts 

2.5. . 3. 5. .1..  127.5. .23 .24 

0.40 

2.5. .3. 5. .1.4..  127.5. .23 .24 

0.56 

2.5. .3. 5. .1.8..  127.5. .23. 24 

0.72 

3. 1 6. . 3. 5. .1..  127.5. .23. 24 

0.32 

3. 16. .3. 5. .1.4. .127. 5. .23.24 

0.44 

3. 16. .3.5. .1.8. .127.5. .23.24 

0.57 

4..3.5..1..127.5..23.24 

0.25 

4..3.5..1.4..127.5..23.24 

0.35 

4.. 3. 5. .1.8..  127.5. .23.24 

0.45 
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Figure  3.16:  Plot  of  Ki  vs.  tp/ts  ratio  keeping  stiffener  width,  stiffener  spacing  and 
stiffener  depth  the  same. 

3.3.6  Effect  of  the  interlaminar  shear  modulus 

The  initial  studies  used  an  interlaminar  shear  modulus  G23  =  Gi3=  G12.  To  investigate 
the  effect  of  this  assumption  analyses  were  performed  using  G23  =  V2  (G13  +  G12).  The 
results  of  these  analyses  are  shown  in  Figures  3.17  and  3.18  and  reveal  that  the  stress 
intensity  factor  is  a  weak  function  of  the  interlaminar  shear  modulii. 


36 


These  Figures  also  highlight  how  the  stress  intensity  factor  asymptotes  to  a  limiting 
value  as  the  crack  length  increases. 


Figure  3.17:  Effect  of  interlaminar  shear  modulus  on  Ki. 


Data  files  3.16..4.6..1..127.5..23.24 


0  10  20  30  40  50 

Crack  length  a  (mm) 

Figure  3.18:  Effect  of  interlaminar  shear  modulus  on  Ki. 


3.4  Summary  of  the  finite  element  analysis 

From  these  results  it  appears  that 

i)  The  solution  for  the  stress  intensity  factors  K]  exhibit  an  asymptotic  behaviour 
as  the  crack  length  increases,  see  Figures  3.8,  3.16  and  3.18. 

ii)  "  The  stress  intensity  factor  Ki  is  a  weak"  function  of  the  stiffener  spacing,  see 

Figures  3.6  and  3.7. 
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iii)  The  stress  intensity  factor  Ki  is  more  strongly  effected  by  stiffener  depth  than 
stiffener  width,  see  Figures  3.12  to  3.15.  As  such  Ki  is  essentially  a  function 
of  the  moment  of  inertia  I  for  the  stiffener  rather  than  the  area  of  the  stiffener. 

iv)  As  a  first  approximation  the  asymptotic  value  of  Ki  appears  to  decay  linearly 
with  increasing  values  of  the  ratio  of  the  patch  thickness  to  the  skin  thickness, 
see  Figure  3.16.  In  comparison  to  the  other  geometrical  panel  cross-section 
parameters,  the  ratio  of  patch  to  skin  thickness  has  the  largest  affect  on  the 
asymptotic  value  of  the  stress  intensity  factor. 

v)  The  stress  intensity  factor  is  a  weak  function  of  the  interlaminar  shear  modulii 
of  the  Boron  epoxy,  see  Figures  3.17  and  3.18. 

3.5  Load  flow  in  a  rib  stiffened  panel 

In  the  previous  sections  we  have  seen  that  the  stress  intensity  factor  Ki  is  a  weak 
function  of  the  stiffener  spacing.  This  can  be  explained  by  understanding  the  load 
flow  in  a  typical  rib  stiffened  panel.  To  this  end  let  us  consider  the  case  of  a  (rib 
stiffened)  panel  containing  a  centrally  located  10  mm  diameter  hole.  The  skin  is  1mm 
thick,  the  stiffeners  are  5  mm  thick  and  are  100mm  apart.  For  simplicity  the  stiffeners 
are  not  tapered  at  the  ends.  The  panel  is  loaded  via  the  end  stingers,  at  locations  A  and 
B,  which  are  each  subjected  to  a  uniform  displacement.  At  this  location  the  skin  is 
assumed  to  be  completely  cracked,  i.e.  the  load  is  only  being  carried  b^  the  stiffeners. 
The  stiffeners  are  uncracked.  The  stiffeners  stop  before  the  opposite  end  of  the  panel 
and  the  load  is  taken  out  of  the  skin  at  points  C  and  D,  see  Figure  3.19.  This 
represents  an  extreme  loading  case  and  is  chosen  to  illustrate  the  load  flow  from  the 
stiffeners  into  the  skin. 

To  visualise  the  load  flow  from  points  C  and  D  to  A  and  B,  after  which  the 
load  is  completely  carried  in  the  stiffeners,  we  plot  the  energy  vector  E  where  the 
components  Ej ,  i=  1,3,  are  defined  as: 


Ei  =  (  J  ay  deij)  ,  E2  =  (  J  a2j  de2j),  E3  =  (  \  cr3j  de3j),  summed  over  j ,  (27) 

Here  cry  and  ey  are  the  stress  and  (mathematical)  strain  tensors  respectively  and  the 
summation  is  over  j,  j=l,3.  For  simplicity  the  components  of  E  are  evaluated  in 
principle  stress  space  and  the  sign  of  the  components  is  taken  as  that  of  the  associated 
principal  stresses. 

Figures  3.19  and  3.20  presents  a  view  of  the  flow  field,  the  E  vector,  in  the  panel.  In 
Figure  3.20  the  view  looks  directly  down  so  that  the  stiffeners  appear  as  a  (thick  black 
)  line.  Both  Figures  clearly  show  the  stress/load  concentration  at  points  A  and  B 
where  the  load  is  entering  the  panel  and  the  stiffener.  A  more  detailed  view  of  this 
flow  is  given  in  Figure  3.21.  Here  we  see  how  the  load/energy  is  flowing  from  the 
stiffener  to  the  upper  skin. 

At  the  edges  of  the  panel,  points  E,  F,  G  and  H  in  Figure  3.20,  we  see  well  defined 
recirculation  zones,  which  indicate  that  the  material  there  is  redundant.  The 
disturbance  of  the  flow  field  caused  by  the  central  hole  is  also  very  clear,  see  Figures 


Note  distort  ion  of  the 
flow  due  to  the  presence 
of  t  he  hole 

•/ 


Figure  3.20.  Plan  view  of  the  rib  stiffened  panel. 


Figure  3.2 1  Detailed  view  of  the  flow  from  locations  A  and  B  into  the  skin. 

Past  the  hole  the  flow  field  again  becomes  dominantly  uni-axial  as  the  disturbance 
caused  by  the  hole  disappears.  The  flow  from  the  stiffener  occurs  close  to  the  end  of 
the  skin,  i.e.  at  points  A  and  B.  Away  from  this  zone,  which  in  Figure  3.20  is 
approximated  by  the  line  L  M,  the  stiffener  has  little  effect  on  the  flow  in  the  panel, 
which  is  dominantly  axial.  Thus,  if  the  crack  is  contained  between  the  stiffeners  and 
has  a  reasonable  edge  distance  from  the  stiffeners,  then  as  a  first  approximation  the 
stiffeners  can  be  ignored,  except  for  estimating  the  local  bending  effects,  when 
estimating  the  residual  strength  of  repaired  rib  stiffened  panels.  Of  course  if  the 
stiffener  was  broken  one  would  have  to  account  for  the  increased  stress  levels  in  the 
panel. 

The  precise  fashion  in  which  the  load/energy  flows  from  the  stiffeners,  at  ends  I  and  J, 
to  points  C  and  D  can  be  seen  in  Figures  3.19,  and  3.20.  More  detailed  views  of  this 
region  are  given  in  Figure  3.22  and  3.23.  Indeed,  Figure  3.23  clearly  highlights  the 
flow  in  the  stiffener,  which  is  dominantly  axial  until  near  the  ends,  i.e.  points  I  and  J, 
where  there  is  a  small  recirculation  zone. 


Figure  3.22.  The  normalised  flow  field  in  a  rib  stiffened  panel. 


Figure  3.23  The  normalised  flow  field  in  a  rib  stiffened  panel. 


This  study  has  implications  for  the  experimental  test  program.  In  this  test  program  it  is 
proposed  to  load  via  the  skin.  We  are  thus  faced  with  the  realisation  that,  other  than 
bending  effects,  the  stiffener  will  only  marginally  effect  the  load  distribution  in  the 
skin.  Thus  the  prime  role  of  stiffeners  will  be  to  induce  bending,  due  to  the  shift  in  the 
neutral  axis. 

This  study  supports/explains  the  previous  numerical  findings  on  the  stress  intensity 
factors  associated  with  repaired  cracks  in  rib  stiffened  (wide)  panels  where  the 
stiffeners  have  been  found  to  have  little  effect. 

3.6  Simple  analytical  formulae  for  predicting  the  limiting  behaviour 

The  finite  element  result  would  indicate  that  as  a  first  approximation  the  simple 
analytical  formulae  for  the  case  when  bending  was  allowed  but  the  stiffeners  were 
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ignored  can  be  used  to  obtain  a  first  estimate  of  the  limiting  stress  intensity  factors 
(Kco)  for  composite  repairs  to  cracked  rib  stiffened  panels. 


c.f.  33.9  for  a  50  mm  deep  stiffener, 
c.f.  38.3  for  a  12  mm  deep  stiffener. 

Table  3.7:  Comparison  of  prediction  for  bending  problems,  Case  3. 
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c.f.  41.9  for  a  50  mm  deep  stiffener, 
c.f.  46.2  for  a  12  mm  deep  stiffener. 

A  f.e.  analysis  found  that  this  value  asymptotes  to  47.4  as  the  stiffener  depth 
approached  zero. 


The  accuracy  of  this  approach  can  be  seen  in  Figure  3.24  which  compares  the  values 
of  K«,  with  the  value  predicted  using  the  formulae  given  in  Section  1.  This 
dramatically  simplifies  the  analysis  and  design  methodologies.  In  this  figure  the  cases 
presented  are  as  given  in  Table  3.5.  In  general  the  maximum  magnitude  of  the 
differences  between  this  simple  approach  and  the  finite  element  results  was 
approximately  10%.  This  can  be  seen  from  Tables  3.6  and  3.7. 

The  accuracy  of  this  fit  can  be  improved  by  allowing  for  the  effect  of  the  stiffener  in 
equation  (8).  In  this  fashion  it  is  proposed  to  replace  equation  (8)  by  the  following 
expression 


i 
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BC  =  a  ymax  (1  -  Kp/ Ks)  (tp  +  tr)/F  (28) 

Where  F  =  F(I,  Is)  where  Is  is  the  I  of  the  stiffener.  One  simple  form  under  evaluation 
takes  the  form 

F  =  I  +  a  Is  (29) 

Where  is  a  constant  to  be  determined  from  the  numerical  results.  This  improved 
formulae  will  be  presented  in  the  next  section 

3.7  Effect  of  riveted  stiffeners 

As  part  of  further  improvements  to  the  CalcuRep  project,  the  effect  of  riveted 
stiffeners  on  the  stress  intensity  factors  for  single  sided  repairs  to  an  aluminium 
substrate  was  investigated.  The  data  presented  in  Section  3  of  Chapter  3  was  analysed 
to  explain  what  the  addition  of  the  riveted  stiffeners  had  on  the  limiting  stress 
intensity  values  of  the  repaired  skin. 

For  the  geometry  shown  in  Figure  3.25,  the  addition  of  the  riveted  stiffeners 
significantly  increases  the  second  moment  of  area  (I)  of  the  panel  cross  section.  For 
the  range  of  geometries  used  in  this  study,  the  increase  in  I  of  the  total  section  of 
between  80  to  1900  times  the  skin  sections  without  stiffeners. 


x 

Figure  3.25  Schematic  of  rib  stiffened  panel  geometry  showing  how  stiffener  depth 
would  significantly  increase  the  second  moment  of  area  (I)  of  the  cross  section. 


In  almost  all  cases  and  geometries  considered  thus  far,  results  of  numerical  analyses 
of  these  panels  have  shown  that  stress  intensity  factor  in  a  repaired  skin  has  almost 
reached  its  limiting  value  for  crack  length  above  40mm  in  length.  Therefore,  the  stress 
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intensity  values  for  the  various  geometries  at  a  crack  length  of  40mm  were  used  to 
compare  with  the  estimates  obtained  using  the  analytical  closed  form  solution  shown 
in  equation  (30). 

Results 

Figures  3.26  and  3.27  show  the  relationship  that  exists  between  the  second  moment  of 
area  and  the  percentage  change  in  the  stress  intensity  for  the  repaired  skins.  These 
results  clearly  show  that  for  a  given  stiffness  ratio  (patch  to  skin),  the  change  in  the 
stress  intensity  factors  with  increasing  moment  of  area  (I)  is  between  5  and  25%. 
Stiffness  ratios  less  than  1  are  rarely  used  and  can  therefore  be  ignored  for  the 
purposes  of  modelling  the  behaviour  using  closed  form  solutions.  If  we  consider  more 
realistic  stif&ess  ratios  of  between  1  and  1.2,  the  limiting  analytical  solution  obtained 
using  K*  for  a  repaired  plate  without  a  stiffener  is  approximately  2  to  10%  larger  than 
that  obtained  for  a  rib  stiffened  plate.  For  a  particular  stiffness  ratio,  the  difference 
increases  with  increasing  second  moment  of  area  although  the  rate  of  change  of  K 
with  I  is  less  marked  at  the  higher  values  (i.e.  very  deep  stiffeners  (37.5  and  50mm)). 


%  Difference  from  limiting  solution  without 
stiffener 


Stiffness  Ratio 
Patch/Skin 


■  0.73  ♦  0.91 


a  1  x  1.15 


+  1.28  •  1.3 


□  1.62  a  1.65 


o  2.1 


Figure  3.26  Plot  showing  differences  in  the  value  of  K  expressed  as  percentage 
differences  in  numerical  analysis  from  closed  form  solution  for  a  range  of  stiffness 
ratios  and  second  moment  of  areas  (I). 

The  relationship  between  Is  and  difference  in  K  (which  we  will  call  the  stiffener 
correction  factor,  SCF)  for  the  single  sided  repairs  was  derived  by  fitting  logarithmic 
curves  to  the  data,  see  Figure  27.  From  these  fits,  a  table  of  coefficients  was  formed, 
see  Table  3.8. 
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SCF  = 


(31) 


to  determine  the  size  of  changes  required  to  the  “Bending  correction  factor”  described 
in  equation  (31).  The  new  bending  correction  factor  accounts  for  the  contribution  of 
stiffener  second  moment  of  area  (Is)  and  enables  a  more  accurate  representation  of  the 
stress  intensity  factor  for  the  repaired  skin  In  this  table  the  data  is  defined  according  to 
the  patch  to  skin  stiffness  ratio  Ej-tr/Eptp. 


Figure  3.27  Plot  of  percentage  difference  in  the  value  of  Kic  ,  defined  as  y,  of 
numerical  analysis  from  closed  form  solution  v’s  x,  the  second  moment  of  areas  (Is) 
for  a  range  of  stiffness  ratios. 


The  new  “bending  correction  factor”  BC,  which  included  stiffener  effects,  for  rib 
stiffened  panels  was  thus  determined  to  be  : 


I  +  SCF 


(32) 


where  I  is  the  I  of  the  unstiffened  panel  and  repair. 
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The  application  of  this  bending  correction  factor  for  the  geometries  analysed  reduces 
the  errors  in  estimation  using  the  closed  form  solution  (CalcuRep*)  to  less  than  3%. 
This  fit  has  been  validated  by  comparison  with  -2400  different  geometrical  cases. 

Table  3.8  Table  of  coefficients  “P”  and  “Q”  for  given  stiffness  ratios  to  calculate 


stiffener  correction  factor  (SCF)  with  varying  Is. 


Stiffness  ratio 
Ertr/Eptp 

“Q” 

0.73 

2.8 

10 

0.91 

2.4 

12 

1 

2.7 

18 

1.15 

1.9 

13 

1.28 

2.2 

18 

1.3 

2.5 

23 

1.62 

1.6 

17 

1.65 

1.9 

22 

2.1 

1.4 

21 

4.  THE  EXPERIMENTAL  TEST  PROGRAM 

For  the  experimental  program  it  was  initially  proposed  that  we  standardised  by  using 
essentially  the  same  configuration  as  used  in  [19],  see  Figures  4.1a)-4.1c).  However, 
these  test  specimens  would  have  riveted  stiffeners  and  not  integral  stiffeners  as  in 
[14].  The  cracks  were  to  be  cut  into  the  skin  of  each  specimen  prior  'to  fabrication. 
The  patch  would  then  be  bonded  on  the  external  surface  of  each  specimen  following 
fabrication  as  this  will  speed  up  the  test  program. 


Y 


(a)  Plan  view  of  wide  panel  test  specimen 
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Table  4.1  Specimen  Identification 


Specimen 

Number 

Skin 

Thickness 

(mm) 

Initial  crack  size 
(mm) 

Stiffener 

broken/unbroken 

1 

1.25 

40.0 

unbroken 

2 

1.25 

40.1 

unbroken 

3 

1.25 

40.14 

broken 

4 

1.25 

40.06 

unbroken 

5 

1.25 

40.24 

broken 

6 

1.25 

40.15 

broken 

7 

1.25 

20.06 

unbroken 

8 

4.15 

41.42 

broken 

9 

4.15 

40.0 

unbroken 

10 

4.15 

19.93 

unbroken 

11 

4.15 

40.19 

unbroken 

12 

4.15 

40.01 

broken 

13 

4.15 

40.03 

unbroken 

14 

4.15 

40.06 

broken 

15 

4.15 

20.01 

unbroken 

16 

4.15 

19.99 

unbroken 

The  thick,  4.15mm,  skin  panels  had  12  ply  thick  patches,  i.e.  ~  1.52mm. 
The  thin,  1.25mm,  skin  panels  had  4  ply  thick  patches,  i.e.  ~0.50  mm. 

For  40mm  cracks  the  patches  were  90mm  wide  and  280mm  long.  , 

For  20mm  cracks  the  patches  were  70mm  wide  and  250mm  long. 


Table  4.2  Specimen  Configuration 


Skin  thickness 
(mm) 

Stiffener  dimensions  (mm) 
width  depth  thickness 

1.25 

23 

50 

2.2 

4.16 

25 

50 

5. 

Following  the  June  99  meeting  with  the  USAF  it  was  decided  to  adopt  a  narrower  test 
specimen,  see  Figure  4.2.  These  specimens  were  manufactured  in  Holland.  An  outline 
of  the  proposed  test  spectra  and  the  specimen  configurations  are  given  in  Tables  4.1 
and  4.2.  In  each  case  the  loads  applied  to  the  specimens  were  chosen  to  give  a  peak 
stress  in  the  skin,  when  the  stiffener  was  intact,  of  120  MPa.  The  specimens  were 
fatigue  tested,  under  constant  amplitude  loading,  in  a  500  kN  servo-hydraulic  MTS 
fatigue  testing  machine  with  a  load  spectra  of  amax  =  120  MPa  and  amjn  =  6  MPa,  i.e. 
an  R  ratio  R  =  amin  /  cfmax  =  0.05. 


Plan  View-  Skin  only 
Dimensions  of  Relevant 
Features  Shown 


•  Dio  =  9,53 


Dia;  6,0 


All  dimensions  in  mm 


a)  View  of  the  narrow  panel  test  specimen  geometry 
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All  dimensions  in  mm  Bolt-On  Grip  -  Plan  View 


y 

b)  Details  of  the  end  fittings 


c)  Details  of  the  rivet  configuration 

Figure  4.2  Schematic  representation  of  the  proposed  narrow  panel  test  geometry 
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4.1  Alteration  to  the  narrow  panel  specimen  design 

During  pre-cracking  of  the  thinner  (1.25mm)  rib  stiffened  panels,  large  cracks  were 
observed  to  initiate  from  the  rivet  holes  closest  to  the  end  grips,  see  Figure  4.3.  This 
occurrence  arose  as  a  result  of  localised  bending  effects  induced  by  the  offset  of  the 
neutral  axis.  Remedial  measures  were  taken  on  the  remaining  thin  and  thick 
specimens.  To  reduce  the  stress  concentration  at  this  location  tapered  aluminium 
doublers  were  introduced,  in  the  vicinity  of  the  grips,  and  the  ends  of  the  vertical 
stiffeners  were  tapered,  see  Figure  4.4.  The  aluminium  doublers  were  bonded  onto  the 
aluminium  skin  using  a  room  temperature  cure  epoxy  adhesive.  These  actions 
inhibited  the  development  of  cracks  at  the  first  few  rivet  holes. 


Figure  4.3  Photograph  of  a  40mm  long  crack  forming  from  the  first  rivet  hole  in  a  thin 
skinned  specimen  with  an  unbroken  stiffener  and  a  20mm  saw  cut  in  the  skin 


4.2  Summary  of  Test  Program 

Two  different  specimen  configurations  were  tested,  see  Tables  4.1  and  4.2,  viz: 

Case  1 :  The  first  test  case  evaluated  had  a  skin  thickness  of -1.25  mm  and  a  width  of 
133mm.  The  stiffener  was  50  mm  deep,  2.2mm  thick  and  the  upper  flange  was  23mm 
wide.  The  area  of  the  skin  Ap  was  thus: 

Ap  =  [133  10'3m]x  [1.25  10'3m] 

=  1.665  lO^m2 

For  this  test  configuration  the  cross  sectional  area  of  the  stiffener  As  was 

As  =  [50  10'3  m]  x  [2.2  10'3  m]  +  [23  10'3  m]  x  [2.2  10'3  m] 

=  1.606  lO^m2 
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Total  Area  (At)  =  [1.862  lO^m2]  +  [1.606  KrW] 

=  3.468  lO^m2 

This  gave  a  total  cross  sectional  Area  (At)  of  3.271  *10"4  m2.  In  this  case  the 
maximum  Fmax  force  needed  to  be  applied  was 

FmaX  =  axAt=  39.5kN 

The  fatigue  test  program  was  performed  with  an  R  ratio  of  0.5.  Consequently,  the 
loads  applied  in  the  fatigue  test  were  Fmin  =  1.96  kN  and  Fmax  =  39.5  kN.  The 
exception  to  this  was  the  first  test  specimen  where  the  minimum  and  maximum  loads 
applied  were  Fmm  =  2.0803  kN  and  Fmx  =  41.6  kN  respectively.  In  this  case  the  stress 
amplitude,  rather  than  the  maximum  stress,  was  —120  MPa. 


Figure  4.4  View  of  patched  surface  and  underside  of  the  test  configuration  with 
cracked  stiffener,  note  adhesive  seepage  causing  the  stiffener  to  be  bonded  to  the  skin. 


Case  2:  The  second  test  case  had  a  skin  thickness  of  4.2  mm  and  a  width  of  133mm. 
The  stiffener  was  50  mm  deep,  5mm  thick  and  the  upper  flange  was  25mm  wide.  The 
area  of  the  skin  Ap  was  thus 

Ap  =  [133  10'3  m]  x  [4.15  10‘3  m]  =  5.519  10A  m2 

For  this  test  configuration  the  cross  sectional  area  of  the  stiffener  As  was 

As  =  [50  10’3  m]  x  [5  10’3  m]  +  [25  10"3  m]  x  [5  10'3  m] 

=  3.75  lO^m2 

this  gave  a  total  cross  sectional  Area  (At)  of  9.27*10'4  m2.  In  this  case  the  maximum 
force  Fmax  needed  to  be  applied  was: 
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Fmax  =  <T  X  At  =  111.2  kN 

The  fatigue  test  program  was  performed  with  an  R  ratio  of  0.5.  Consequently,  the 
loads  applied  in  the  fatigue  test  were  Fmjn  =  5.56  kN  and  Fmax  =  1 1 1.2  kN 

4.3  Test  results 

Case  1:  From  the  results  for  specimen  number  2,  see  Figure  4.5,  which  had  an 
intact  stiffener  and  was  tested  under  a  stress  amplitude  of  120  MPa,  we  found  that  the 
crack  growth  rate  was  constant  and  that 

da/dN  =  5.578  104  mm/cycle  =  21.96  10'6  in/cycle 

In  this  case  from  the  results2  given  in  [22],  page  8.9-84,  for  an  R  ratio  of  0.02  and  a 
plate  thickness  of  between  0.02”  to  0.2”,  we  found  that  the  experimental  crack 
growth  rate  corresponds  to  an  (experimental)  value  of  AK  ~  11  MPa  Vm.  This 
compares  quite  favourably  with  the  values  of  AK  of  12.2  MPa  Vm,  obtained  using  the 
semi-analytical  formulae,  and  values  of  11.4  MPa  Vm,  for  a  14mm  long  crack,  and 
12.6  MPa  Vm,  for  a  40mm  long  crack,  obtained  using  3D  finite  element  analysis. 

For  the  subsequent  tests  the  stress  amplitude  was  ~  1 14  MPa.  The  growth  rate  again 
was  constant  and  da/dN  ~  3.72  104  mm/cycle.  The  experimental  crack  growth  rate 
corresponded  to  an  (experimental)  value  of  AK  ~  10.0  MPa  Vm  compared  with  a 
predicted  value  of  11.6  and  10.9  MPa  Vm  and  12.04  MPa  Vm,  obtained  using  3D 
finite  element  analysis,  for  a  14mm  long  crack  and  a  40  mm  long  crack  Respectively. 

When  the  stiffener  was  broken  the  stress  in  the  skin  under  the  stiffener  increased  from 
120  MPa  to  ~238  MPa.  In  this  case  we  obtained  a  predicted  value  of  AK  of  23.9  MPa 
Vm,  obtained  using  the  semi-analytical  formulae,  and  values  of  21.9  MPa  Vm,  for  a  14 
mm  long  crack,  and  22.6  MPa  Vm,  for  a  40  mm  long  crack,  were  obtained  via  a  3D 
finite  element  analysis. 

The  experimental  test  results  were  quite  repeatable,  see  the  data  for  specimens  3,  5 
and  6  as  given  in  Figure  4.7.  For  specimens  3  and  6  we  again  obtained  a  constant 
growth  rate  with  a  da/dN  of  ~  2.4  *10"3  mm/cycle,  whilst  for  specimen  5  we  obtained 
a  growth  rate  of  ~  2.0  *10‘3  mm/cycle.  In  this  case  using  the  results  given  in  [22]  this 
growth  rate  gave  a  AK  of  ~  22.3  MPa  Vm,  for  specimens  3  and  6,  and  ~  20.5  MPa  Vm 
for  specimen  5.  These  values  are  in  good  agreement  with  the  numerical  predictions. 

During  the  fatigue  tests  it  was  found  that  failure  generally  initiated  at  rivet  holes 
outside  of  the  patched  area.  The  location  of  the  initiation  site  was  essentially  random 
and  appeared  to  associated  with  the  initial  fabrication  of  the  specimen.  Indeed, 
specimen  7  failed  outside  of  the  test  area  without  any  crack  growth  under  the  patch 
and  specimen  1  failed  outside  of  the  test  area  with  minimal  crack  growth  under  the 
patch,  see  Figure  4.6.  In  the  latter  case  the  growth  rates  were  distorted  by  the 
extensive  cracking  occurring  outside  of  the  patch. 


2  In  this  case  the  growth  law  can  be  approximated  as  da/dn  =  1.64  10  9  (AK)2'35 
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Figure  4.5.  Crack  growth  for  a  1.25mm  thick  skin  specimen,  number  2,  with  an  intact 
stiffener,  stress  amplitude  of  120  MPa 


Figure  4.6.  Crack  growth  for  specimens  1  and  4,  1.25mm  thick  skin  with  an  intact 

stiffener,  max  stress  of  120  MPa 
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Figure  4.7  Crack  growth  for  a  1.25mm  thick  skin  specimen  with  a  broken  stiffener, 

specimens  3,  5  and  6. 

Case  2:  For  the  second  test  case  the  test  specimens  had  a  skin  thickness  of  ~ 
4.15mm  and  a  stiffener  thickness  of  5mm.  Specimens  with  two  initial  crack  lengths, 
viz:  20  and  40  mm,  were  tested.  In  both  cases  we  (again)  saw  that  when  the  crack  was 
under  the  patch  the  crack  growth  rate  was  essentially  constant,  see  the  data  for 
specimens  9,  10  and  15  in  Figure  4.8.  Indeed,  as  was  expected  the  crack  growth  rates 
for  both  initial  crack  configurations  were  (essentially)  identical  and  gave  a  growth  rate 
of 

da/dN  ~  8.2  10"4  mm/cycle  ~  32.4  10"6  in/cycle 

For  this  plate  thickness  we  obtain  from  the  results  given  in  [22],  page  8.9-86,  for  an  R 
ratio  of  0.05  and  a  plate  thickness  of  0.2”,  i.e.  ~5mm,  which  is  close  to  the  test 
configuration  an  experimental  crack  growth  rate  of  -32.4  10'6  in/cycle  and  a 
corresponding  (experimental)  value  of  AK  -  15.95  MPa  Vm. 

The  present  analysis  predicts  a  value  of  AK  of  17.1  MPa  Vm,  obtained  using  the  semi- 
analytical  formulae.  Finite  element  analysis  gave  values  of  -12.7  MPa  Vm  and  15.6 
for  a  14mm  and  a  40mm  crack  respectively.  We  thus  see  that  the  analytical  formulae 
compare  quite  favourably  with  the  experimental  test  results. 

For  specimens  11  and  13  there  was  extensive  seepage  of  the  adhesive  through  the 
fasteners  with  the  result  that  the  stiffener  became  completely  bonded  to  the  skin,  see 
Figure  4.10.  This  meant  that  the  crack  growth  rates  for  these  specimens  were 
dramatically  reduced,  see  Figure  4.8,  and  could  not  be  related  to  the  simple  formulae 
as  the  specimen  configuration  now  differed  substantially  from  the  test  requirement. 
However,  the  crack  growth  rates  for  both,  modified,  specimens  were  consistent  and 
very  similar,  see  Figure  4.8.  The  phenomena  also  occurred  in  specimens  with  a 
cracked  stiffener,  see  Figure  4.4.  However,  in  these  cases  the  stiffener  was  not 
bridging  the  crack  and  as  a  result  did  not  transmit  load  across  the  crack. 
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In  contrast  to  the  1.2mm  thick  specimens  the  thicker  (4.2mm)  specimens  generally 
failed  due  to  crack  growth  under  the  repair,  i.e.  inside  of  the  patched  area. 

When  the  stiffener  was  broken  the  stress  in  the  skin  under  the  stiffener  increased  from 
120  MPa  to  ~204  MPa.  The  tests  were  extremely  repeatable,  see  Figure  4.9,  and  it 
was  found  that  crack  growth  rate  was  (again)  constant  with  a  da/dN  ~  4.1  10'3 
mm/cycle  ~  1.6  10‘4  in/cycle,  from  which  we  obtain  a  AK  of  ~26  MPa  Vm.  It  should 
be  stressed  that  this  value  of  AK  may  have  large  error  bounds  since  the  data  for  this 
combination  of  crack  growth  rate,  plate  thickness  and  test  frequency  is  very  scarce.  In 
this  case  we  predicted  a  value  of  AK  of  -28.1  MPa  Vm,  obtained  using  the  semi- 
analytical  formulae  for  a  long  crack.  This  compared  well  with  the  results  of  a  three 
dimensional  finite  element  analysis  which  gave  values  ranging  from  19.9  MPa  Vm, 
for  a  14  mm  long  crack,  to  22.4  MPa  Vm,  for  a  40  mm  long  crack.  This  configuration 
had  the  largest  discrepancy  between  the  analytical  and  the  finite  element  results. 


“  Si/SHSi  SUSS! 

Figure  4.10  Close  up  of  adhesive  seepage,  showing  stiffener  bonded  to  the  skin. 

4.4  Failure  surface 


For  the  thin  skinned  specimens  neither  the  patch  or  adhesive  failed  at  any  stage  dining 
the  tests.  This  was  also  true  for  those  thicker  (4.15  mm)  skinned  specimens  which  had 
an  unbroken  stiffener,  even  when  the  crack  grew  across  the  full  width  of  the 
specimen.  For  the  thicker  skinned  specimens  which  had  a  broken  stiffener  it  was 
found  that  when  the  crack  grew  across  the  full  width  of  the  specimen  the  patch  system 
failed  in  a  complex  fashion,  reminiscent  of  that  described  in  Section  2.2,  in  that  the 
first  layer  of  the  patch  remained  attached  to  the  metal  indicating  that  the  failure 
involved  interlaminar  failure  in  the  composite  patch. 
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This  can  be  seen  by  looking  at  the  underside  of  the  failed  specimen  shown  in  Figure 
4.11.  The  adhesive  used  to  bond  the  patch  to  the  specimens  was  green  in  colour.  Yet, 
as  can  be  seen  in  Figure  4.11,  the  first  ply  of  the  patch  is  still  largely  attacked  to  the 
skin  and  there  are  only  small  regions  of  green  to  be  seen  on  the  underside  of  the 
specimen.  This  (again)  illustrates  the  need  to  consider  interlaminar  (intra-ply)  failure 
when  considering  the  load  bearing  capacity  of  the  repair.  In  this  picture  the 
delaminated  patch  has  been  placed  on  the  other  side  of  the  stiffener  so  as  to  show  that 
after  (interlaminar)  failure  the  patch  was  still  in  one  piece. 


Delaminated  patch  has 
been  placed  on  the  other 
side  of  the  stiffener  so 
as  to  show  that  after 
(interlaminar)  failure 
the  patch  was  still  in 
one  piece. 


‘"r  ~ "  V?  'i  1 


Note:  Fibres  still  bonded 
and  on  most  of  the  area 
of  underside  of  the  patch 
there  is  no  evidence  of 
the  (green)  adhesive  used 
to  bond  the  patch  to  the 
specimen. 


Figure  4.1 1  Close  up  of  failure  surface  showing  the  first  ply  attached  to  the  skin  as 
well  as  the  doubler  being  in  a  single  piece  after  failure. 


4.5  Comments 


This  test  program,  when  taken  in  conjunction  with  the  results  of  a  three  dimensional 
finite  element  study  of  the  test  configurations),  has  confirmed  that  the  semi-analytical 
formulae  are  capable  of  representing  the  behaviour  of  composite  repairs  to  rib 
stiffened  panels. 

The  test  program  has  also  shown  that  after  patching  the  stress  intensity  factor  rapidly 
asymptotes  to  a  constant  value,  which  is  independent  of  crack  length,  so  that  the  crack 
growth  rate  becomes  constant  and  is  independent  of  crack  length.  This  phenomena  is 
true  regardless  of  whether  the  stiffener  is  broken  or  not.  However,  the  absolute  value 
of  the  stress  intensity  factor,  and  hence  the  growth  rate,  is  dependent  on  the  integrity 
of  the  stiffener,  i.e.  whether  or  not  it  is  broken. 


5.  EFFECT  OF  VARIABLE  ADHESIVE  THICKNESS 

Until  now  we  have  assumed  that  the  adhesive  has  a  constant  thickness.  However,  for 
field  or  depot  level  repairs  the  adhesive  thickness  may  vary  due  to  a  variety  of 
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application  and/or  fabrication  related  problems.  The  problem  of  variable  adhesive 
thickness  can  be  associated  with: 

i)  The  inability  to  maintain  an  appropriate  vacuum/pressure  throughout  the 
cure  cyle; 

ii)  A  non-uniform  temperature  distribution  over  the  region  being  heated.  This 
can  occur  due  to  the  presence  of  spars  or  ribs  which  may  act  as  a  heat  sink. 

iii)  Adhesive  seepage  through  the  fastemers; 

iv)  Geometric  discontinuities  associated  with  the  structure  being  repaired.  The 
repair  to  the  Fill  wing  pivot  fitting  [16]  is  a  good  example  of  such  a 
problem; 

v)  The  use  of  low  temperature  cure  adhesives. 

Let  us  now  address  the  problem  of  variable  adhesive  thickness.  In  this  section  we  will 
show  that  the  Glinka  extension  to  the  Hart  Smith  design  formulae  can  be  used  for 
those  cases  when  the  adhesive  thickness  is  not  constant. 

Before  examining  the  inelastic  behaviour  of  joints  with  variable  adhesive  thickness  let 
us  first  consider  the  elastic  response.  Consider  the  joint  as  shown  in  Figure  5.1,  the 
governing  equations  for  this  joint  are: 


w,  -u2  =ty . 

From  Hooke’s  law  for  the  adherends  we  find  that: 


and  for  adhesive 


7]  —  E]h]ulx , 
T2  =  E2h2u2x 


x  =  Gy  . 


,(33) 

(34) 

)  V 

(35) 


(36) 

(37) 

(38) 


where  Ev2  and  hV2  are  Young’s  modulus  and  thicknesses  of  adherends,  G  and  t  are 
the  shear  modulus  and  thickness  of  the  adhesive.  Other  parameters  are  as  shown  in 
Figure  5.1. 

Let  us  consider  the  case  when  the  adhesive  thickness  is  variable  i.e.  t  =  t(x) .  In  this 
case  equations  (33)  to  (38)  produce  the  resultant  differential  equation,  viz: 

(y'L  =P2(*h*  (39) 

where  (3 2  (x)  =  -j- ((£,/z,  )"'  +  {E2h2 )"' ) 

t\X) 

In  the  case  of  a  constant  adhesive  thickness  equation  (39)  reduces  to  the  well-known 
equation  for  the  shear  stress  in  adhesive,  viz 


61 


(40) 


where  p2  =  y((£j/Zj)  1  +  {E2h2 ) 

For  a  constant  thickness  adhesive  joint  the  concept  of  a  characteristic  length  X  is 
commonly  used.  This  distance  X  is  given  by  the  formulae: 

X  =  3/p  (41) 

and  represents  the  length  of  the  adhesive  playing  the  major  role  in  the  load  transfer 
process  transferring  the  loading. 


Ax  x 

Displacement  Equilibrium  Equations  ) 


Figure  5.1  Schematic  of  the  adhesive  joint 

5.1.  Analytical  solutions  for  variable  adhesive  thickness 

Equation  (39)  has  no  a  general  analytical  solution.  However  by  using  the 
transformation  z  =  (yt)^  j  Jt  it  can  be  reduced  to  the  special  Riccati  equation  [27]: 

z,x+az2  =  p2(x)  (42) 

with  a  =  1 .  This  equation  (14)  can  be  solved  for  a  wide  class  of  function  p(x), 
particularly  when  this  function  has  the  form  P2  (x)  =  bxm . 
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Consider  the  following  example.  Let  the  adhesive  thickness  in  the  vicinity  of  the  edge 
is  distributed  as  t(x)  =  t0xA .  In  this  case  (32  =  P2(x)  has  the  form 

(3  2(x)  =  b2x~A,  (43) 


where  b2  =  —  ^Exh\)  '  +  (£2^2) ') 

* 0 

We  thus  require  a  solution  to  the  following  equation: 


(y t)M  =  bx~4yt  at  x0  <  x  <  00  (44) 


The  solution  of  this  equation  results  in  the  following  expression  for  the  adhesive  shear 
stresses 


t(x)  =  — j  (c  ■  cosh (b/x)  +  a  sinh {bjx)) 


t0x 


(45) 


The  constants  a  and  c  in  this  solution  can  be  found  from  the  boundary  conditions, 
which  can  be  taken  in  the  following  form 


PEX\ 


Exh{  +E2h. 


•dx 


(46) 


which  represents  the  balance  of  forces  acting  on  the  adherends  and 


P 


E2h2 


= M,, 


(47) 


From  the  boundary  conditions  we  find  that 


and 


c  = 


PEA 


j EAi+E2h2  G  E2h2 


(48) 


P  +  —  sinh(i>/  x0)-c-  cosh(&/x0 ) 


a  =  ■ 


E2h2  x0 


sinh(Z)/x0 )  -  b/x0  cosh(Z>/x0 ) 


(49) 


In  the  special  case  when  £,/z,  =  E2h2  -  Eh  we  have  c  =  0  and 
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p 

a  =  - — (sinh(Z>/*0 )  -  b/x0  cosh (b/x0 ))”’ 
and  the  expression  for  the  shear  stresses  becomes 


t(x)  = 


P  G 


smh(b/x) 


Eh  t0x3  b/x o  cosh(b/^0)-sinh(h/x0) 


(50) 


(51) 


The  criterion  proposed  by  Hart-Smith  for  failure  of  adhesive  joints  can  be 
written  in  the  form 

W^<Wcr  (52) 

where  Wm &  is  the  maximum  value  of  the  strain  energy  density  in  the  adhesive,  and 
Wcr  is  the  critical  value  for  the  energy  determined  from  standard  materials  tests. 

In  this  problem  the  maximum  strain  energy  density  occurs  at  the  edge  of  the 
joint,  i.e.  x  =  x0,  where  the  shear  stress 


P__G _ 1 

Eh  t0x2  b/x0  coth(b/x0 )  - 1 


where  coth(b/x0)  is  the  hyperbolic  cotangent. 


(53) 


Consider  the  case  when  the  value  of  b/x0  »1,  this  corresponds  to  small 
changes  of  the  thickness  of  the  adhesive  on  the  characteristic  length  of  the  stress 
distribution,  then 


In  this  limiting  case  the  solution  for  reduces  to 


_  1  P2  _  1  P2 
4  Eht0xl  4  Eht 


(55) 


This  formula,  which  is  independent  of  the  stress  strain  relationship  for  the  adhesive, 
coincides  with  the  maximum  strain  energy  density  for  an  adhesive  with  a  constant 
equal  to  the  thickness  at  the  edge  x  =  x0.  This  means  that  the  present  solution 
asymptotes  to  the  solution  for  the  constant  thickness  of  the  adhesive  as  x0  reduces  to 
zero. 


Let  us  now  we  consider  the  case  when  the  adherends  are  an  equal  thickness 
aluminium  alloy  with  hx  =  h2  =  3 mm ,  and  E}  =  E2  =  73  GPa.  The  thickness  of  the 
adhesive  is  distributed  in  accordance  with  a  power  law  in  the  following  form 
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reveals  that,  as  a  first  estimate,  Glinka’s  hypothesis  is  valid,  i.e.  Wm ax,  /  W* max,  =1  see 
Figure  5.3.  Consequently,  we  can  use  the  elastic  solution  to  estimate  the  strain  energy 
density  when  the  joint  experiences  significant  plasticity.  The  validity  of  this 
approximation  for  other  joint  configurations  will  be  investigated  in  Section  5.2. 


W„ 


WT 


e 

max 


Y,/Yc 


Fig  5.3.  Ratio  of  the  maximum  strain  energy  for  the  elastic-perfectly-plastic  solution 
to  value  obtained  via  an  elastic  analysis  for  a  non-uniform  adhesive' thickness. 

This  figure  also  reveals  that  the  energy  density  is  a  function  of  the  local  geometry  of 
the  joint.  As  such  the  Hart-Smith  formulation  cannot  be  used  to  predict  the  strength  of 
the  joint.  In  this  case  we  must  use  the  Glinka  extension  to  the  Hart-Smith  formulation. 

5.2  Effective  mechanical  properties  of  micro-inhomogeneous  joints 


If  mechanical  geometric  characteristics  of  the  adhesive  joint  vary  from  point  to  point 
and  the  characteristic  scale  of  these  changes  is  much  smaller  in  comparison  with  the 
characteristic  length  of  the  joint,  then  the  joint  can  be  considered  as  micro- 
inhomogeneous  structure  (continuum).  The  problem  of  averaging  of  the  mechanical 
properties  of  micro-inhomogeneous  continua  has  long  history.  However,  its 
asymptotic  nature  has  only  recently  been  understood  and  the  corresponding 
mathematical  formulation  is  presented  in  [20-21], 

For  the  adhesive  joint  the  Lagrangian  A  can  be  written  as 


2A  = 


_ 1 _ 

(£,A.  )-'+(£  A  )■’ 


(56) 


This  expression  represents  the  sum  of  the  strain  energy  densities  for  adherend  and 
adhesive  layers.  Using  Euler’s  equation  it  can  be  shown  that  this  Lagrangian  produces 
the  same  differential  equation  for  yt ,  i.e.  equation  (39),  as  presented  previously. 
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The  effective,  or  averaged,  mechanical  characteristics  of  a  micro-inhomogeneous 
continuum  can  be  determined  by  calculating  the  Legendre  transformation,  A  (jt.p),  of 
A,  with  respect  to  (yt)  x .  In  this  case  we  find  that 

2A‘(y t,p)=  {(El  (x)hx  (x))-1  +  (E2  (x)h2  (*))"  )p2  +  f(,)>  (57) 

The  averaged  Lagrangian  (A’)  is  then  the  integral  of  A*  ,  viz: 

2(  A' )  =  ((£,*, )-  +  (E2k2 )"  y  +  (^j(ylf  (58) 


Here  the  expression  (•)  represents  the  average  value  of  a  function,  or  material 
properties,  over  a  representative  volume: 


(59) 


where  VR  is  much  bigger  than  the  characteristic  scale  of  changes  of  the  geometry  or 
material  properties  of  the  adhesive  joint. 


Calculating  the  inverse  Legendre  transformation  we  obtain  the  averaged  Lagrangian 

(A) 


2(A)  =  ((£,*,  r  +  (£A  r  )"  [M„  f  +  (f  }w2  (60) 

This  equation  reveals  how  to  calculate  effective  parameters  of  the  joint.  To  illustrate 
this  methodology  let  us  consider  a  simple  example.  Let  the  adhesive  thickness  of  the 
joint  has  two  equal  fractions  of  thicknesses  t)  and  t2  such  tx  « t2 .  Then  the  effective 
thickness  t  of  the  joint  is  thus 


t  —  l{tx  +  t2 1 )  »  2tj 


(61) 


which  reveals  that  the  joint  will  behave  as  a  joint  whose  thickness  is  double  the 
minimum  thickness  of  the  adhesive. 


Let  us  now  consider  when  the  thickness  of  the  adhesive  is  randomly  distributed  in 


accordance  with  a  uniform  probabilistic  law  between  the  values  ?min  and  tm!a.  In 
this'case  we  obtain  an  effective  thickness  t  where 


tofonaxAmin) 


(62) 


The  same  rule  is  also  applied  for  the  stiffness  of  the  adherends  Eh  .  However,  note 
that  for  the  shear  modulus  of  the  adhesive  we  have  a  different  law  for  calculating  the 
effective  value.  The  effective  shear  modulus  of  the  adhesive  can  be  found  as 
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G=(G)  (63) 

After  averaging  the  mechanical  properties  and  the  adherends  and  adhesive  thicknesses 
we  subsequently  solve  the  governing  differential  equation: 

(yO,xx  =  $ 2Jt  (64) 

where  p2  is  now  a  constant  and  represents  the  effective  (or  averaged)  value  for  the 
joint.  As  a  result  of  this  formulation  it  follows  that  for  joints  where  the  adhesive 
thickness  varies  from  point  to  point,  and  the  characteristic  scale  of  these  changes  is 
much  smaller  in  comparison  with  the  characteristic  length,  the  value  of  Wmax  is  again 
independent  of  the  shear  modulus  of  the  adhesive. 

5.3  The  effect  of  variable  adhesive  thickness  and  material  non-linearity 

Having  established  the  Glinka  hypothesis  to  be  applicable  to  bonded  joints  with 
global  bending  prohibited  let  us  now  numerically  consider  the  effect  of  both  variable 
adhesive  thickness  and  adhesive  plasticity.  We  have  seen  that  the  energy  density  is  a 
weak  function  of  the  adhesive  shear  modulus.  It  should  thus  follow  that  the  hypothesis 
will  also  hold  if  the  adhesive  thickness  is  variable. 

To  further  validate  this  hypothesis  let  us  consider  a  symmetric  double  lap  joint  with  a 
lower  adherend  (half)thickness  Ti  of  3mm,  an  upper  adherend  thickness  T?  of  2mm, 
and  Ei  =  E2  =  73,000  MPa.  The  adhesive,  which  was  assumed  to  be  FM73,  was 
allowed  to  be  either  one  or  two  layers  thick,  i.e.  0.1mm  or  0.2  mm  thick,  and  had  a 
Young’s  modulus  of  1890  MPa  and  a  Poisson’s  ratio  of  0.35.  Whilst  the  adhesive  was 
assumed  to  have  a  post  yield  slope  E’  of  250  MPa  and  the  yield  stress  of  the  adhesive 
was  allowed  to  vary,  see  Table  5.1.  Furthermore,  at  the  critical  end  of  the  joint  the 
thickness  ta'  was  assumed  to  vary  according  to  the  following  equation 

ta1  =  ta  +  A  cos(  n  x/2xi)  for  0  <  x  <  X]  (65) 

ta1  =  ta  for  X  >  X] 

where  ta,  A  ,  x,  and  Xi  are  the  nominal  adhesive  thickness,  the  amplitude  of  the 
perturbation  of  the  thickness,  the  distance  from  the  end  of  the  joint  and  a  characteristic 
length.  In  this  initial  study  we  chose,  on  the  basis  of  past  experience,  Xi  to  be  3  mm, 
and  varied  the  values  of  A.  At  this  stage  only  monotonic  loading  was  considered  and 
to  increase  the  amount  of  yielding  the  remote  stress,  in  the  lower  adherend,  was  taken 
to  be  144  MPa.  The  results  of  this  analysis  are  summarised  in  Table  5.1. 

To  further  illustrate  the  validity  of  Glinka’s  hypothesis  for  variable  adhesive  thickness 
let  us  consider  the  case  when  the  adhesive  is  two  layers  thick,  i.e.  ta  =  0.2  mm,  and  A 
=  0.05  mm.  This  corresponds  to  a  perturbation  in  the  adhesive  thickness  equal  to  25% 
of  the  base  line  thickness.  In  this  case  the  distance  over  which  the  perturbation 
occurred,  i.e.  xj,  was  varied  from  1.5  to  5mm  and  the  results  are  shown  in  Table  5.2. 
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Table  5.1  Percentage  difference  in  Vw  and  values  of  W,  g  =  360  MPa 


Yield  stress  and 
adhesive  thickness 

A  (mm),  the  local  increase  in  the  adhesive  thickness. 

0.025 

0.05 

0.075 

cfy  =  40  MPa 
ta  =  0.2mm 

18.1 

W  =13.27 

17.0 

W  =11.54 

15.5 

W  =10.22 

Gy  =  68  MPa 
ta  =  0.2  mm 

13.9 

W  =14.65 

11.5 

W=  13.12 

9.5 

W  =11.78 

Gy  =  40  MPa 
ta  =  0.1  mm 

25.0 

W  =16.95 

22.7 

W  =14.80 

20.5 

W  =12.87 

Gy  =  68  MPa 
ta  =  0.1  mm 

21.6 

W  =18.76 

17.7 

W  =16.75 

15.33 

W  =14.59 

From  this  work  we  see  that,  for  monotonic  loading,  although  variable  adhesive 
thickness  changes  the  stress  and  the  energy  fields  Glinka’s  hypothesis  is  still  valid. 
This  means  that,  for  the  present  class  of  problems,  even  if  there  is  variability  in  the 
thickness  of  the  adhesive  bond  the  energy  field  Wmax  can  be  computed  using  only  the 
elastic  solution  of  the  repair.  The  strength  of  the  repair,  provided  failure  is  due  to 
failure  of  the  adhesive,  can  be  then  estimated  from  a  purely  linear  elastic  analysis. 
This  will  occur  at  the  load  level  when  Wmax  =  Wcrjt  •  This  finding  significantly 
simplifies  the  design  process  for  repairs  with  a  variable  thickness  adhesive. 


Table  5.2  Percentage  difference  in  Vw  and  values  of  W  (MPa),  A  =  0.05  mm 


Yield  stress 
and  post 

yield  slope 

Remote  stress  g  (MPa) 

144 

200 

240 

288 

360 1 

Gy  =  68  MPa 
ta  =  0.2  mm, 
xi  =  1.5  mm 

0.0 

W=2.58 

0.4 

W=4.90  . 

2.5 

W  =6.81 

6.0 

W  =  9.11 

11.1 

W  =12.7 

Gy  =  40  MPa 
ta  =  0.2  mm, 
X]  =  1.5  mm 

8.7 

W=3.05 

5.5 

W=5.49 

2.2 

W  =6.85 

9.0 

W  =8.53 

16.7 

W=ll.l 

Gy  =  68  MPa 
ta  =  0.2  mm, 
xf=  3  mm 

0.0 

W=2.68 

0.5 

W=5.08 

2.8 

W  =7.03 

6.4 

W=9.39 

11.5 

W  =13.12 

Gy  =  40  MPa 
ta  =  0.2  mm, 
Xi  =  3  mm 

7.3 

W=3.09 

4.3 

W=5.59 

3.5 

W  =7.08 

9.4 

W=8.80 

17.0 

W  =11.54 

Gy  =  68  MPa 
ta  =  0.2  mm, 
xi  =  5  mm 

0.0 

W=2.80 

0.8 

W=5.28 

3.2 

W  =7.28 

6.9 

W=9.70 

12.1 

W  =13.53 

3.7 

W=5.75 

3.0 

W  =7.32 

9.9 

W=9.09 

17.5 

W=11.91 
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